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PREFACE 


This document contains the proceedings of the joint NASA/George Washington 
University Workshop on Failure Analysis and Mechanisms of Failure of Fibrous 
Composite Structures. The workshop was held March 23-25, 1982, at Langley Research 
Center and was attended by over 80 researchers from government agencies, aircraft 
manufacturers, and universities. The objectives of the workshop were to assess the 
present state of the art of failure analysis and current design practices, to iden- 
tify deficiencies in these technologies, and to identify directions for future 
research leading to a verified failure analysis capability for fibrous composite 
structures. 

The workshop was organized into the followin=g six sessions : 

1. Design Technology I 

2. Design Technology II 

3. Failure Mechanisms and Failure Criteria 

4. Delaminations and Stress Concentrations 

5. Compression Failure and Transverse Strength 

6. Effects of Defects and Experimental Methods 

Certain materials are identified in this publication in order to specify ade- 
quately which materials were Investigated in the research effort. In no case does 
such identification imply recommendation or endorsement of the product by NASA, nor 
does it imply that the materials are necessarily the only ones or the best ones 
available for the purpose. In many cases equivalent materials are available and 
would probably produce equivalent results. 


Ahmed K. Noor 
Mark J. Shuar t 
James H. Starnes, Jr. 
Jerry G. Williams 
Workshop Coordinators 
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INTRODUCTION 


Over the past several years, considerable effort has been made by the research, 
development, and design communities to understand the mechanisms of failure of 
fibrous composite materials and the effect of these mechanisms on the performance of 
structural components made from these materials. Experience has shown that the use 
of fibrous composite materials for new design applications requires a better under- 
standing of their failure mechanisms so that safe, reliable structural designs can be 
developed. The full weight— saving potential of composite structures can only be 
realized by understanding their failure mechanisms, developing the capability to 
predict accurately the effect of these mechanisms on structural performance, and 
transferring this knowledge to designers of composite structures. 

The joint NASA/George Washington University Workshop on Failure Analysis and 
Mechanisms of Failure of Fibrous Composite Structures was held at NASA Langley 
Research Center, March 23-25, 1982, to assess the current state of the art of failure 
analysis for composite structures and of current design practices, to identify 
deficiencies in these technologies, and to identify directions for future research 
leading to a verified failure analysis capability for composite structures. This 
document contains copies of the workshop presentations and Includes some of the 
subsequent discussions. An extensive bibliography on various aspects of the failure 
of composite materials and structures is also included. Some general comments on 
topics discussed during the workshop are presented in this introduction. 


Local Damage and Its Effect on Design 

The strength of composite structures has been shown to be reduced significantly 
by local damage and discontinuities. Traditionally, the concern for durability of 
metal structures has focused on tension-fatigue behavior. However, for composite 
structures, both tension- and compression-loading responses are important, and a 
large percentage of the reduction in strength due to damage can be measured in static 
tests. Experience has shown that defects and local damage can exist In a laminate 
prior to their visual detection. An assessment of the effect on strength of such 
damage is Important to the structural designer. 

Nondestructive damage assessment methods are commonly used to monitor damage 
growth. Surface inspection techniques include the use of edge replicates and moire 
fringe patterns. Ultrasonic attenuation, vibrothermography, and radiography also are 
used when damage measurements more definitive than surface inspection are desired. 
Test and analysis studies are used by structural designers to ensure that the struc- 
tural strengths and stiffnesses are not reduced to unacceptable levels due to damage 
propagation between inspection periods. 


Damage Growth and Failure Criteria 

Because of their importance to practical applications of composites, damage 
growth and failure criteria are current research topics. Damage growth in composite 
laminates can result from either matrix- or fiber -dominated failure mechanisms. 
Matrix-dominated failure mechanisms include ply delamination and matrix cracking. 
Failure mechanisms for fibers include fiber fracture and fiber microbuckling. Damage 


growth that results from matrix-dominated failures has received considerable 
attention. By contrast, limited work has been done on damage growth resulting from 
fiber failures. For matrix-cracking and fiber fracture mechanisms, fracture me- 
chanics has been used to establish a criterion for damage growth based on the strain 
energy release rate for a laminate. Using this criterion, the damage grows when the 
strain energy release rate reaches a critical. value calculated from the various 
fracture mode contributions. 

Failure criteria for composite structures are being studied widely. Failure 
criteria for composite materials subjected to multiaxial states of stress have 
traditionally been operational rather than mechanistic. These criteria are not 
derived from mechanics principles but are empirical statements of the material 
failure process. Nevertheless, composite failure criteria are useful tools for 
material characterization and for estimating the load-carrying capacity of composite 
structures. Initially, familiar and easy-to-use failure criteria for isotropic 
materials were generalized for anisotropic composite materials. Examples of these 
criteria are: maximum stress criteria, maximum strain criteria, and distortional 

energy criterion (e.g., Tsai-Hlll, Hoffman). Later, tensor polynomial criteria for 
composites were proposed. The maximum stress and strain criteria are popular because 
they are operationally simple. For multiaxial states of stress, these criteria are, 
in effect, curve-fitting schemes and have no analytic foundation. The distortional 
energy and tensor polynomial criteria have built-in generality and mathematical con- 
sistency. Most of the existing failure criteria agree with one another when used to 
predict uniaxial and pure shear strength data along material S3mimetry axes. However, 
the failure criteria do not always agree when applied to combined stress states. 
Better correlation between the analytical predictions using failure criteria and 
experiments have been observed when the failure criteria are applied to the entire 
laminate rather than to the unidirectional constituent layers. As the understanding 
of specific failure mechanisms for composite structures increases and as new mech- 
anisms are identified, new failure models and criteria may be required. 


Delamination 

Delamination is a common failure mode in composite structures because inter- 
laminar strength is the lowest strength in a laminated material. Interlaminar 
stresses can be caused by out-of-plane loads, surface impact, or interlaminar de- 
fects. These stresses are also induced near any free edge (including notches and 
through-the— thickness discontinuities) of a laminate subjected to inplane loading. 

Delamination is a design consideration for composites in the same manner as 
tensile strength, compressive strength, and fatigue life. Therefore, the analysis 
and control of delamination are being widely studied. Analytical and computational 
models have been developed for the accurate determination of interlaminar stresses 
near free edges prior to delamination. These models are based on the use of pertur- 
bation techniques, approximate three-dimensional elasticity solutions, finite dif- 
ferences, and more recently, finite elements. Classical techniques of linear elastic 
fracture mechanics have been used for predicting the onset and growth rates of 
delamination. 
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The number of techniques to control delamination is limited, but the following 
three procedures are being used: 

1. Proper choice of the laminate construction (e.g., modifying the stacking 

sequence) to reduce the interlaminar stresses 

2. Material selection (e.g., resins with improved fracture toughness) 

3. Reinforcement through the laminate thickness (e.g., stitching) 

Also, good quality assurance programs can eliminate large voids as sources of delami- 
nation. The use of the aforementioned suppression techniques can result in changing 
the mode of failure from delamination to another type. 


Global Instabilities and Component Failure 

The failure characteristics of structural components are generally more complex 
than those of the relatively simple coupon specimens often used for material charac- 
terization. Structural components can be subjected to combined loading conditions 
with both in-plane and out-of-plane loading components and can respond to the applied 
loads by structural response phenomena, such as component buckling. Extrapolating 
coupon data to large-scale structural components may be inadequate. Local damage or 
a local failure can cause the load distribution in a structural component to change 
significantly. The loading conditions and response phenomena of structural compo- 
nents can cause large interlaminar stresses in the component which cause failure to 
occur by mechanisms not observed in simpler coupon tests. 

The differences between the failure analyses for coupons and for structural 
components are also being examined. Heuristic failure analyses developed from coupon 
data may be useful for estimating failure of a structural component subjected to the 
same loading conditions as the coupon tests, but these analyses may need reevaluation 
when other loading conditions or failure mechanisms are likely to occur. Heuristic 
analyses based on uniaxial in-plane loading conditions may not be adequate for more 
complex loading conditions. A simple analysis for predicting delamlnat ion-related 
failures may not be accurate for other failure modes (e.g., shear crippling) if 
delamlnatlon is suppressed. No reliable failure analysis has yet been developed for 
combined in-plane and out-of-plane failures such as the skin-stiffener separation 
failure mode observed for stiffened panels with buckled skins. Research is currently 
under way to develop failure analyses suitable for large-scale structural components, 
but additional experience with large components and with more complex responses and 
loading conditions is needed before reliable, verified failure analyses are readily 
available. 


Analysis Methodology for Predicting Progressive Failure 

Although methods are currently available for predicting the strength of com- 
posite laminates subjected to unidirectional loading, analysis methodology is still 
evolving for predicting the progressive failure of composite structures subjected to 
multidirectional loading. Failure of a laminated composite structure may start at a 
local defect that is characterized on the mlcromechanlcs level and may grow to the 
sublaminate, laminate, and subcomponent levels. Analysis methodology for predicting 
the progressive failure of composite structures could start at the micromechanics 


3 


level and trace the failure through the aforementioned levels. Different analytical 
models would be needed at each level. At the first level the model might Include 
(a) distribution of local damage and defects, (b) consequent local states of stress 
and deformation and (c) nonlinear material response. At the sublaminate and laminate 
levels, the analytical model should be capable of predicting interlaminar stresses. 

At all levels, simple models are needed which capture the major response charac- 
teristics. However, a reliable verified failure analysis for a composite component 
can only be achieved if the sensitivity to damage in the various regions of the 
structure is known. 


Damage-Tolerant Design Technology 

Most high-strength graphite-epoxy composite laminates lack the ductile behavior 
at high strains which is characteristic of metals. Composite structures respond to 
local strain concentrations (such as those that exist in the vicinity of holes and 
impact damaged regions) by localized failure. This local failure reduces the local 
stiffnesses and transfers high strain concentrations into adjacent regions, which may 
cause the local damage to propagate. Recent studies indicate that residual strength 
of impact-damaged laminates can be Improved by Increasing the fracture toughness of 
the matrix material, and some investigators anticipate that a high-f allure-strain 
fiber could also improve damage tolerance characteristics. 

Specific matrix properties that constitute toughness are not well established; 
however, a high failure strain (greater than 4 percent) for the matrix with nonlinear 
stress-strain response at high strains may be a necessary but not a sufficient con- 
dition. Resin materials with high fracture toughness (referred to herein as tough 
resins) may Improve the impact damage tolerance characteristics of brittle composites 
by suppressing the delamination mode of failure. For brittle composites, delami- 
nation growth is the dominant failure mechanism observed in tests. The performance 
of tough-resin composites under cyclic loading is not yet well established. Local 
shear crippling appears to be the dominant failure mechanism for tough-resin lami- 
nates with damage or discontinuities such as holes. 

Achievement of damage tolerance for heavily loaded composite aircraft structures 
can be addressed at several levels of structural complexity. At the material level, 
tough resins and high-failure-strain fibers are being developed. Improved damage 
tolerance can also be achieved through the use of innovative structural configura- 
tions designed to arrest and contain the damage within a region of acceptable size. 
There have been recent advancements in damage containment for both tension- and 
compression-loaded applications. Finally, damage tolerance can be addressed at the 
global structural level by designing efficient sti^ctures capable of carrying the 
redistributed loads following local, regional failure. In summary, a combination of 
advancements including material improvements, damage containment configuration 
development, and global-local redistribution design methods are needed to meet the 
stringent damage tolerance requirements typical of heavily loaded aircraft components 
such as a wing. 


Future Directions of Research 

This workshop emphasized the need for a systematic study to understand the 
failure mechanisms of composite structures from the micromechanics to the component 
levels and to develop a reliable and verifiable capability for predicting progressive 
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failure. As new resins and fibers are developed, the effect of material nonlinearity 
on damage onset and propagation needs to be assessed. For design purposes, simple 
computational models that capture the multilevel aspects of failure are needed. 
Different models may be needed at each level. 
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COMMERCIAL TRANSPORT AIRCRAFT 
COMPOSITE STRUCTURES 


John E. McCarty 

Boeing Commercial Airplane Company 
Seattle, Washington 


This presentation will take a look at the role that analysis plays In the 
development, production, and substantiation of aircraft structures; the types, 
elements, and applications of failure that are used and needed; the current appli- 
cation of analysis methods to commercial aircraft advanced composite structures, 
along with a projection of future needs; and some personal thoughts on analysis 
development goals and the elements of an approach to analysis development. 


mTRQDUCTIQN 

ANALYSIS ROLE 
ANALYSIS ELEMENTS 
COMPONENT ANALYSIS 
ANALYSIS DEVELOPMENT 
SUMMARY 


8 



The analysis of structure Is the only truly feasible means of substantiating 
the strength, fatigue life, and damage tolerance of a majority of commercial aircraft 
structures. This is true simply because only a minimum number of critical locations 
on the aircraft can be validated by full-scale testing. Test data at the coupon, 
detail, structural element, and subcomponent levels provide half of the Information 
required to establish the critical margin of safety. The other half Is obtained by 
structural analysis. Therefore, the confidence and credibility of the analysis tools 
used are very critical to the acceptance of analysis as the prime means of structural 
substantiation of commercial aircraft structure. 


AREAS OF DISCUSSION 

FAILURE ANALYSIS 
FAILURE MECHANISMS 


ANALYSIS ROLE FOR COMMERCIAL AIRCRAFT: 

"THE PRIMARY MEANS OF STRUCTURAL 
SUBSTANTIATION FOR CERTIFICATION OF 
COMMERCIAL AIRCRAFT IS BY ANALYSIS . 

IT IS EXPECTED THAT THE ANALYSIS WILL 
BE SUPPORTED BY THE APPROPRIATE TEST 
EVIDENCE" 


This figure dramatizes the idea of why the analysis is really the only means of 
fully substantiating the structure of large aircraft. This in no way implies that 
the same is not true for small aircraft. The large number of structural details that 
must be reviewed in producing any aircraft requires that all Information required by 
both halves of the margins of safety equations be available. 


Transport Aircraft Composite 

Structures 


• Large primary structure composite development 
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The Federal Aviation Regulations (FAR 25) (ref. 1) established the requirements 
for structural certification of commercial aircraft. The list of numbered para- 
graphs shows those most pertinent to this discussion on advanced composite appli- 
cations. In addition to the FAR 25, the FAA has issued an advisory circular for 
composite structure (ref. 2). The advisory circular provides guidance to both 
Industry and the FAA as to the acceptable means of compliance with the FAR 25 for 
advanced composites. 


CURRENT REGULATIONS 
INVOLVED IN COMPOSITE 
CERTIFICATION 


MATERIAL AND f25.603 

PROCESS SPECIFICATIONS \25.605 


MATERIALS 

FABRICATION METHODS 


MATERIAL PROPERTIES 


'25.613 MATERIAL STRENGTH PROPERTIES AND 

DESIGN VALUES 

25.615 DESIGN PROPERTIES 

25.619 SPECIAL FACTORS 


PROOF STRENGTH 


f25.305 STRENGTH AND DEFORMATION 

125.307 PROOF OF STRUCTURE 


DAMAGE TOLERANCE 


25.571 FATIGUE EVALUATION OF FLIGHT STRUCTURE 


ADVISORY CIRCULAR AC NO. 20-107 DATED 7/10/78 
SUBJECT: COMPOSITE AIRCRAFT STRUCTURE 


11 



These lines selected from FAR 25 (ref. 1) and the advanced composite advisory 
circular (ref. 2) illustrate FAA's recognition of the major role that structural 
analysis plays in the substantiation of the strength and damage tolerance of a 
commercial aircraft structure. Both the regulations and the guidelines recognize 
the relationship of the types of material, experience with material, and structural 
configuration and concepts, as well as the supportive contribution of test data. 

FAR 25 & ANALYSIS 


FAR 25. 307... "STRUCTURAL ANALYSIS MAY BE USED ONLY IF 
THE STRUCTURE CONFORMS TO THOSE FOR WHICH EXPERIENCE 
HAS SHOWN THIS METHOD TO BE RELIABLE..." 


FAR 25.571(B) DAMAGE-TOLERANCE (FAIL-SAFE) EVALUATION... 
"THE DETERMINATION MUST BE BY ANALYSIS SUPPORTED BY TEST 
EVIDENCE AND (IF AVAILABLE) SERVICE EXPERIENCE..." 


ADVISORY CIRCULAR & ANALYSIS 
AC NO. 20-107 


3b... "THE EXTENT OF TESTING AND/OR ANALYSIS AND THE DEGREE 
OF ENVIRONMENTAL ACCOUNTABILITY REQUIRED WILL DIFFER FOR 
EACH STRUCTURE DEPENDING UPON THE EXPECTED SERVICE USAGE. 
THE MATERIAL SELECTED. THE DESIGN MARGINS. THE DATA BASE 
AND EXPERIENCE WITH SIMILAR STRUCTURE AND ON OTHER FACTORS 
AFFECTING A PARTICULAR STRUCTURE." 


5... "THE STATIC STRENGTH OF COMPOSITE DESIGN SHOULD BE 
DEMONSTRATED THROUGH A PROGRAM OF COMPONENT ULTIMATE LOAD 
TESTS IN THE APPROPRIATE ENVIRONMENT. UNLESS EXPERIENCE 
WITH SIMILAR DESIGNS. MATERIAL SYSTEMS AND LOADINGS IS 
AVAILABLE TO DEMONSTRATE THE ADEQUACY OF THE A?!ALYSIS 
SUPPORTED BY SUBCOMPONENT TEST." 


12 



This particular statement should be of interest to the workshop attendees since 
several of the presentations relate to the idea of material system allowables. 


ADVISORY CIRCULAR & ANALYSIS 


Ab. "THE MATERIAL SYSTEM ALLOWABLES SHOULD BE 
ESTABLISHED ON THE LAMINATE LEVEL BY EITHER 
TEST OF THE LAMINATE OR BY TEST OF THE LAMINA 
IN CONJUNCTION WITH A TEST VALIDATED ANALYTICAL 
METHOD" 


An understanding of which analysis tools relate to connnercial aircraft safety 
and which to lifecycle economics is required to form a proper perspective of the 
analysis substantiation requirements. This chart attempts to provide a view of this 
idea. The confidence required for safety is attained by conservative application of 
simple analysis methods supported by extensive testing, or realistic application of a 
variety of analysis tools of varying sophistication supported by appropriate levels 
of test evidence. 


DESIGN/ANALYSIS AND CERTIFICATION ELEMENTS 


SAFETY 


ECONOMICS 


STATIC 

STRENGTH 

EXTERNAL LOADS ANALYSIS 
MAT^L 8 SIR. ALLOWABLES 
FAILURE CRITERIA & ANALYSIS 

DAMAGE 

TOLERANCE 

FLAW GROWTH ANALYSIS 
RESIDUAL STRENGTH ANALYSIS 
INSPECTION REQUIREMENTS 

FLUTTER 

MARGIN 

STIFFNESS REQUIREMENTS 
FLUTTER ANALYSIS 

MAINTENANCE 

COST 

DURABILITY - FATIGUE ANALYSIS 
INSPECTION COST 
REPAIR COST 

PRODUCTION 

COST 

MAT'L COST 
DESIGNED'IN COST 
MANUFACTURING METHODS & COST 
INSPECTION METHODS & COST 
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To better describe the requirements of damage-tolerant concepts, two charts 
are presented with the following statements. (1) The anticipated damage that may 
arise from normal aircraft operations can include fatigue damage, manufacturing 
and/or maintenance flaws, or errors in undetected accidental damage. (2) The damage 
tolerance design structure must also provide a very substantial means of protecting 
against accidental damage sustained in flight as a result of such things as engine 
breakup, hall, and bird or other types of Impact damage. In these particular cases 
it is expected that the damage will be found upon completion of the flight. 


DAMAGE TOLERANCE 


"THE STRUCTURE MUST BE DESIGNED IN SUCH A WAY 
THAT ANY DAMAGE INCURRED FROM NORMAL OPERATION 
IS DETECTABLE BEFORE THE STRENGTH OR STIFFNESS 
OF THE STRUCTURE FALLS TO AN UNACCEPTABLE LEVEL." 
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In order to establish for this discussion the separation of the fatigue aspect 
of analysis from that of damage tolerance analysis, this chart gives a graphical 
representation of this separation. Simply stated, damage-tolerant considerations 
must include the idea of damage detectability. Therefore, the idea acts to provide 
a real design and application separation. Once detectable, the time to grow to 
critical is the area needing the damage tolerance growth analysis. 


DISCIPLINED ANALYSIS APPROACH TO DURABILITY/DAMAGE 
GROWTH DESIGN CONSIDERATIONS 



CYCLES - N 
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This chart presents a breakdown by category of some of the detailed aspects of 
the requirements of damage tolerance analysis. Categories 2 and 3 are the analysis 
areas of concern in this discussion. However, Category 1 was the means of meeting 
the damage tolerance requirements for the Boeing 727 advanced composite elevator. 


DAMAGE TOLERANCE ANALYSIS REQUIREMENTS 



STRUCTURAL CATEGORY 

TECHNIQUE OF 

TECHNOLOGY 




ASSURING SAFETY 

CONTROL METHOD 
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The ability to configure a new aircraft Is based on the state of the art of the 
available analysis tools. The structure configuration, structural concept, and 
material selection are all evaluated analytically. The payoff of this analytical 
effort Is the future successful production and operation of new aircraft. The 
ability to change and Improve the aircraft Is reduced with this advancement through 
the development, production, and operational cycle. The state of the art and the 
ease of application of the analysis of tools applied early In the design cycle are 
significant factors In developing a successful aircraft. 

DESIGN/ ANALYSIS ROLE IN STRUCTURAL 
AIRFRAME DEVELOPMENT 


• RESEARCH & DEVELOPMENT 

IDENTIFY NEW STRUCTURAL PAYOFFS 
EVALUATE MATERIALS S STRUCTURAL CONCEPTS 
GUIDE THE TEST PLANNING & DATA ANALYSIS 

• PRELIMINARY DESIGN & PRODUCT DEVELOPMENT 

TRADES - WEIGHT 8 COST 
CONFIGURATION 8 CONCEPTS 

• AIRCRAFT PROGRAM 

STRUCTURAL SIZING 
DRAWING RELEASE 

FORMAL ANALYSIS - FOR CERTIFICATION 
MARGINS OF SAFETY 
TEST/ANALYSIS CORRELATION 


ANALYSIS SUPPORTED DECISION EFFECT ON COST AND/OR WEIGHT 



AIRCRAFT DEVELOPMENT TIME 
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I An element In utilizing any analytical capability Is the fact that In a large 

design program there must be a discipline to the tools used. Each company publishes 
design manuals, stress manuals, and analysis programs In various forms to ensure 
some uniformity In analysis procedures. The key then Is the disciplined procedures 
that provide the benefits shown on this chart. 


ANALYSIS SUPPORT OF AIRCRAFT PROGRAMS 


KEY - DISCIPLINED PROCEDURES 


PROVIDES ACCEPTABLE TOOLS AT ALL LEVELS 

USABLE BY LARGE NUMBERS OF NON-SPECIALISTS 

CAUSES ATTENTION TO BE FOCUSED ON CRITICAL PARAMETE^^S 

PROVIDES COMMON QUANTIFIABLE BASE FOR DECISION MAKING 

CORRELATES FLEET AND TEST EXPERIENCE 

ESTABLISHES BASE FOR AIRFRAME SUSTAINING 
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This list will be familiar to all those in the structures technical community. 
The application of analysis of tools in each of these levels of structural signifi- 
cance will be discussed. These analysis classifications have been separated for 
discussion purposes into the areas of requirements analysis and capability analysis. 
The idea is that the structure is required to perform a set of structural functions 
or requirements, and the structure has been designed to provide a level of capability 
to meet those performance requirements. 


ANALYSIS PROCEDURE NEEDS 



20 



» 


A collection of analytical tools across all the analysis disciplines (strength, 
fatigue, and damage tolerance) provides a means of quantifying the structures 
capability margin (margins of safety) In each area. This and the following charts 
are aimed to Illustrate the disciplined approach In production analysis of structures. 


DESIGNED & REQUIRED STRUCTURAL PERFORMANCE 


PERFORMANCE EVALUATION: 


MARGIN OF SAFETY = 


CAPABILITY OF STRUCTURE 
STRUCTURAL REQUIREMENTS 


- 1 - 
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A discipline analysis procedure can be used not only to establish the margins of 
safety but also as a means of quantifying structural performance against weight and 
cost. The analytical tools should end up clearly defining both the requirements of 
the structure and the capability of the structure being designed. 


COMMON APPROACH TO STATIC, FATIGUE, DAMAGE GROWTH 
AND RESIDUAL STRENGTH DESIGN 


ANALYSIS 

ULTIMATE 

STRENGTH 

FATIGUE 

STRENGTH 

DAMAGE 

GROWTH 

RESIDUAL 

STRENGTH 

REQUIREMENT 


'"^REQUIRED 

nxp 

'^REQUIRED 

^RS 

CAPABILITY 


'"‘^DETAIL 

'^^’^DETAIL 

^RS 

MARGIN 

ftu 

‘'^DETAIL j 

''^REQUIRED 

“^•^DETAIL j 

‘^^'^REQUIRED 

•"rs j 

^RS 
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The idea of an analysis model being part of the requirements analysis does not 
imply a single model, but refers rather to all the analysis tools required to 
establish the requirements half of the margins of safety equation. The analysis 
model for structure can Include the following: finite element, interlaminar stres- 

ses, laminate stress analysis at the detailed stress analysis level, fatigue damage 
model, damage growth model, and fracture mechanics analysis. 
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The capabilities analysis recognizes that the structural form, the type of 
applied load, the material used, and the expected failure mechanism or mechanisms 
must be part of any evaluation. The material properties can be as simple as test 
data to develop the properties, l.e., tension, compression, or shear, or they can be 
established by analysis procedures, such as in a simple column analysis. Because of 
the intrinsic variability of material properties, a reliability analysis of the test 
data must be part of any procedure. These forms of capability analysis and the 
requirement analysis are only an example. 


ANALYSIS PROCEDURE 
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This section of the discussion will look at where we are in the development 
cycle of commercial aircraft advanced composite structure. I will try to express 
two ideas: (1) the application of the appropriate analysis tools to the component 
and recognition of the background that has been established to date, and (2) a 
projection of future needs. 


COMMERCIAL AIRCRAFT DEVELOPMENT OF COMPOSITE 
STRUCTURE 

CURRENT 

• CONTROL SURFACES 

• EMPENNAGES 

FUTURE 

WINGS 

BODY 
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This figure is shown here to establish the current production status of advanced 
composites on commercial aircraft. The primary control surfaces , of graphite/epoxy , 
have been designed and certified using analysis tools available to date . An example 
of the level of this analysis is illustrated in the discussion that follows on the 
727 advanced composite elevator. 
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The 727 elevator, which is Illustrated In the figure, shows an upper and lower 
skin of honeycomb design consisting of one ply of cloth at ±45° and a chordwlse 0° 
tape direction. The spars are laminate with ±45° dominant webs and with 0° dominant 
in the chords. The few ribs noted are primarily ±45° with some 0° in the chords and 
honeycomb webs. This design has been fully certified by the FAA. 


Elevator Structural Arrangement 
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The following lists of analysis considerations are to direct attention not to 
the detailed analysis but to those facts that were considered in the interplay 
between this structure and current analysis capability or the analysis capability 
available at the time of the development and certification of the 727 elevator. The 
analysis used was just at the level required to meet those certification and design 
requirements. Since, as noted on the chart, only ultimate strength certification 
requirements were needed, this is all that is addressed by the analysis tools 
available. 


727 ELEVATOR - ANALYSIS CONSIDERATIONS 
DESIGN CRITERIA 

SIMILAR STIFFNESS TO ALUMINUM ELEVATOR 
DAMAGE TOLERANCE CATEGORY n 
TEMPERATURE RANGE -75°F TO 180°F 
MOISTURE CONTENT 1% ± .1% BY WEIGHT 

RESULTS OF CRITERIA ON DESIGN 

LOW ULTIMATE DESIGN STRAIN 
NO FATIGUE CONSIDERATIONS 
NO DAMAGE TOLERANCE REQUIREMENTS 

CERTIFICATION CONSIDERATIONS 

1®^ OF A TYPE MUST TEST TO ULTIMATE 
WITH ACCOUNTABILITY FOR MOISTURE & TEMPERATURE 

ULTIMATE STRENGTH ONLY 
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The tools used In the analysis of the 727 elevator were those that have been 
used in the past on some large control surfaces, l.e., on the Boeing 747. Most 
Boeing control surfaces are currently analyzed by hand. The 727, being the first 
generic of its type at Boeing, incorporated the use of the finite element method. 

The Boeing 767 control surfaces were primarily analyzed by hand. For the reasons 
noted, other than the first generic, the finite element method was used on this 727 
elevator. The other analysis concerned the stability and the ultimate strength of 
both the surface panels and the web of the front spar. Simplified analysis was used 
to adequately establish the buckling capability of these honeycomb panels. The 
ultimate strains were simply the strains associated with the mav i miim strain in the 
most critical direction compared to point design allowables of the specific layup 
developed by test. 


727 ELEVATOR ANALYSIS REQUIREMENTS 

FINITE ELEMENT ANALYSIS 

NO RIBS BACKING UP HINGE FITTINGS 

METHOD OF ACCOUNTING FOR MOISTURE/TEMPERATURE 

CHECK STIFFNESS SIMILARITY 

DETAIL ELEMENT ANALYSIS (STRENGTH CHECK) 

SURFACE PANELS 

STABILITY - MODIFIED ALUMINUM METHOD 
ULTIMATE STRAIN - MAX. PRINCIPAL 

FRONT SPAR 

WEB SHEAR STABILITY - MODIFIED AL. METHOD 
ULTIMATE STRAIN IN SHEAR 
CHORD STRAINS - ULTIMATE STRAINS T & C 
JOINT STRENGTH BEARING, TENSION 
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This figure simply illustrates the finite element model used. Note the grid 
refinement aft of the hinge locations. The extra attention in this area was due to 
the fact that this was the first Boeing design that did not use a backup rib at each 
hinge location. The local load redistribution from the hinges was of particular 
interest. 


Elevator Finite Element Model 
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This figure and the next one Illustrate the simplicity of the data available for 
basic analysis. These strain cutoffs and modulus charts, along with some bearing 
allowables and the laminate testing that was part of the ancillary test program of 
the 727 elevator, are all that was provided for both the preliminary design as well 
as the formal analysis. 
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ELASTIC MODULUS 
AND ALLOWABLE STRAINS 

IN-PLANE SHEAR MODULUS 
AND ALLOWABLE STRAINS 

Tape 

Graphite/epoxy, BHS 8-212 Type III Class 1, per BAC 5562 
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One of the unique considerations of the analysis model was the use of the 
characteristics of the material, the finite element analysis, and the validation 
by full-scale ultimate load test of the strain distribution to establish a means of 
analytically accounting for the effects of moisture and temperature in the ultimate 
load test results. Since material acts In a linear manner and the response to both 
load-induced strains and moisture- and temperature- Induced strains Is also linear In 
response In fiber-dominated structure, the two can algebraically be added to estab- 
lish accurate strain In the material at all locations by analysis. The distribution 
for the loads analysis Is verified by the full-scale test. By comparing this anal- 
ysis strain level to those established by laminate level testing, a margin of safety 
can be adequately described. 


727 ELEVATOR ANALYSIS CONSIDERATION 
ULTIMATE LOAD TEST MOISTURE AND TEMPERATURE ACCOUNTABILITY 


• MATERIAL RESPONSE 

• FINITE ELEMENT ANALYSIS 

• STRAIN DISTRIBUTION VALIDATION 


• LINEAR STRAIN RESPONSE 

• FIBER DOMINATED AND MAJOR 
LOAD PATH DESIGN 
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This table illustrates the result of Including the effects of moisture and 
temperature in the requirements part of the margin of safety calculation. 


PRINCIPAL STRAINS, ALLOWABLES, AND MARGINS OF SAFETY - ATLAS PANEL 1135, LOAD CASE 125 


PANEL 

ENVIRONMENTAL 

1 PRINCIPAL STRAINS* 

ALLOWABLE PR INC. 

STRAINS'* 

MARGINS OF SAFETY | 

SURFACE 

CONDITION 

*-MAX 

^MIN 

^MAX 

^MAX 

^MIN 

Y 

'max 

^MAX 

^MIN 

Y 

'max 

EXTERIOR 

70° DRY 

1989 

-2404 

4393 

5310 

-4930 

8499 

1.67 

1.05 

.93 


70° WET 

1545 

-2294 

3840 

4260 

-4410 

6690 

1.76 

.92 

.74 


-75° DRY 

1754 

-2073 

3826 

7060 

-5670 

8156 

3.03 

1.74 

1.13 


-75° WET 

1463 

-2116 

3579 

5850 

-5790 

6422 

3.00 

1.74 

.79 


140° DRY 

2406 

-2833 

5239 

4450 

-4600 

8800 

.85 

.62 

.68 


140° WET 

1938 

-2700 

4638 

3600 

-3900 

6755 

.86 

.44 

.46 


180° DRY 

2693 

-3127 

5821 

3980 

-4400 

8608 

.48 

.41 

.48 


180° WET 

2221 

-2990 

5211 

3260 

-3650 

6770 

.47 

.22 

.30 

INTERIOR 

70° DRY 

1502 

-1629 

3132 

4930 

-5310 

8499 

2.28 

2.26 

1.71 


70° WET 

1006 

-1451 

2456 

4410 

-4260 

6690 

3.38 

1.94 

1.72 


-75° DRY 

560 

-1098 

1657 

5670 

-7060 

8156 

9.13 

5.43 

3.92 


-75° WET 

327 

-1183 

1509 

5790 

-5850 

6422 

16.71 

3.95 

3.26 


140° DRY 

2185 

-2087 

4273 

4600 

-4450 

8800 

l.ll 

1.13 

1.06 


140° WET 

1679 

-1898 

3577 

3900 

-3600 

6755 

1.32 

.90 

.89 


180° DRY 

2591 

-2364 

4956 

4400 

-3980 

8608 

.70 

.68 

.74 


180° WET 

2084 

-2175 

4259 

3650 

-3260 

6770 

.75 

.50 

.59 


* MICROSTRAIN 
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The following figures represent two of the test articles that were used (1) to 
demonstrate the analysis of the structure for the load strains and for the thermal 
and moisture strains and (2) to show that the analysis gave good correlation in 
fiber-dominated structure. The results of these two tests are illustrated on the 
three graphs. In addition to the 727 and the 737 tests, I have shown the similar 
compatibility of the linear strain effect previously discussed. These tests, along 
with the test data at the coupon structural element and subcomponent levels, which 
included the effects of moisture and temperature, formed the basis both for demon- 
strating the requirements of capability analysis validity and for providing the 
capability data to form the margins of safety used in the certification of both the 
727 and 737 advanced composite components. 






X DESIGN ULTIHATE 


Lower skin panel 















The failure of the front spar of the 727 elevator during test was for a Boeing 
design condition, not an FAA required condition. The failure occurred In a combined 
load area of the upper chord of the spar at the hinge fitting. The combined load In 
the fastener area was demonstrated to be the reason for this failure. A simple 
Interaction curve showed that the failure should have occurred. Good correlation 
between (1) the stress analysis at both the finite element level and the detailed 
hand analysis level and (2) the test data and Interaction curve shown gave a very 
good correlation with the failure noted. 


727 Elevator Failure Analysis 

Combined load on front spar chord 
hinge fitting attachment 


Failure of front spar in a Boeing 
fail-safe test 
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727 Elevator failure analysis 
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0.0046 nr)/nm,(in/in) , from room temperature, 
dry test specimens (fig 71) 


of 1.00 based on foo- = 709 MPa (102 800 Ibf/in ) from room temperature, 

dry test specimens (fig 68) 


Interaction Curve Bearing versus Tension Bypass Strain 
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A look at the 737 structure will give some insight into the levels of 
information and analysis capability that were required for this structure. Again, 
simple tools could be used since the structure is primarily designed by stiffness. 
This shows a breakdown of the structural elements of the 737 horizontal stabilizer, 
which consists of cover panels (co-cured I-stiffened panels), laminate front and 
rear spars, and honeycomb composite ribs. 


Stabilizer Inspar Structural Arrangement 
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- The 737 analysis considerations look at the design criteria Impact, and here 
again stiffness design dominated and therefore produced low strain levels. The 
damage tolerance requirements are shown to be category 2. This requires that the 
horizontal stabilizer be designed for large detectable damage with no growth. The 
same moisture/temperature characteristics and analysis procedure were established 
as were used for the elevator. Therefore, the 737 requires certification for both 
damage tolerance and ultimate strain. Again, a finite element analysis was applied 
to the horizontal stabilizer primarily because of the design, which carries only the 
two spars through the center section. Therefore, a very significant shear lag and 
load distribution problem needs careful analysis. Also, by using finite elements 
and cutting or removing structure to simulate damage tolerance requirements, the 
analysis procedure was easier to perform. 

737 STABILIZER - ANALYSIS CONSIDERATIONS 
DESIGN CRITERIA 

SIMILAR STIFFNESS TO AL, STABILIZER 
DAMAGE TOLERANCE CATEGORY #2 

CONTROL STRAIN TO ELIMINATE DAMAGE GROWTH 
(I.E., CATEGORY #3) 

TEMPERATURE/MOISTURE (SAME AS ELEVATOR) 

RESULTS OF CRITERIA ON DESIGN 

LON ULTIMATE DESIGN STRAINS 
NO FATIGUE CONSIDERATIONS 
NO FLAW GROWTH CONSIDERATIONS 
TOLERANT TO DISCRETE DAMAGE 

CERTIFICATION CONSIDERATION 

OF A TYPE MUST BE TESTED ULTIMATE 
WITH ACCOUNTABILITY FOR MOISTURE R TEMPERATURE 
ULTIMATE STRENGTH 
DAMAGE TOLERANCE (CATEGORY #2) 


737 HORIZONTAL ANALYSIS REQUIREMENTS 


FINITE ELEMENT ANALYSIS 

SPAR CARRY THROUGH TO CENTER SECTION ONLY 
METHOD OF ACCOUNTING FOR HOISTURE/TEMPERATURE 
CHECK STIFFNESS SIMILARITY 
DAMAGE TOLERANCE ANALYSIS 
MAJOR DAHAGE/FAIL SAFE 

DETAIL ELEMENT ANALYSIS (STRENGTH CHECK) 

SURFACE PANEL ANALYSIS 

STABILITY, STRENGTH S DAMAGE TOLERANCE 

SPARS 

WEB - STABILITY, STRENGTH S DAMAGE TOLERANCE 
CHORDS - STABILITY, STRENGTH S DAMAGE TOLERANCE 
JOINT - STRENGTH 

RIBS 

STRENGTH S SKIN ATTACHMENT 
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The breakdown of the finite element model is shown below. Notice the finer 
grid in the shear lag load distribution region toward the Inboard end. 


SUBSTRUCTURE 3 UPPER SURFACE 


SUBSTRUCTURE 4 TRAILING 
EDGE SUPPORT STRUCTURE 


SUBSTRUCTURE 1 INSPAR BOX 



SUBSTRUCTURE 5 CENTER SECTION 
SUPPORT STRUCTURE 


SUBSTRUCTURE 2 LOWER SURFACE 


737 Compotha Horizontal Stabilizer Finite Element Model Substrvcture Definition 
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Buckling Analysis 




BUCKLING INTERACTION- PANEL 69.9 cm (27.5 in) x 8 

VIPASA STRINGER BAYS 


Analysis Approach 


• Skin panels buckle at 
42% ultimate load 

• At ultimate 

• Bending carried by 
stringers and spar 
chords 

• Shear carried by 
skin 

• Stiffness matching 
constraints 

• Strain cutoffs 


STAB STAB 
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The pin removal in the lug in the rear spar simulated a fatigue failure in the 
aluminum inboard of the graphite structure. The center section Is an arrangement of 
a truss network which carries the torsional loads in the diagonals of the truss and 
the bending in the front and rear spars of the truss. A failure in this inboard 
section, particularly at the adjacent lug, would cause a significant load redistri- 
bution in the rear spar of the composite structure. This test, therefore, was felt 
to simulate adequately this type of possible failure mechanism. 


• Maximum positive bending 

• Rear spar lower pin removed 

• Front spar lower pin removed 

• Maximum negative bending 

• Front spar upper pin removed 

• Rear spar upper pin removed 


Damage Tolerance Test 
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The test set-up for the full-scale ground test is shown. Strain surveys to 
limit load were performed for four load cases. Thermal linkage functional tests 
at high and low temperatures were performed to verify the thermal compensating 
linkage. The lateral stiffness of the elevator attachments was determined. 

The test box was spectriun loaded for one-half lifetime to verify no intrinsic 
damage growth. Small cuts, impact damage, and damaged fastener holes were introduced 
Into the test box. The box was spectrum loaded for one full lifetime to demonstrate 
that visible damage will not propagate during one lifetime of spectrum loading. Four 
ultimate load conditions were applied to the test box. Following these tests, a 
sequence of tests were performed in which the lug pins were removed to simulate 
center section lug failures. The lightning protection system was verified by sub- 
jecting an outboard tip section to a lightning test. 


Full-Scale Ground Test 


PUMMY RIGHT-HAND 
STABILIZER 


-PRODUCTION CENTER SECTION 
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(NEW 737 COMPONENT) 
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SUPPORT 
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This photograph shows the full-scale test specimen mounted in the support jigs 
Loads were applied by a system of pads to simulate spanwise and streamwise load 
distribution. 







This photograph shows the stabilizer's center section interface. Attachment of 
the stabilizer is made with five bolts, three at the rear spar and two at the front 
spar. 



Test Configuration-Failure 


• Removed pin in upper lug rear spar 

• Simulated center section failure 

• Applied load case 401 0 (down bending) 

e 67% DDL required 

e Failure occurred at 61% DDL 
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It can be seen from these illustrations of the failure of the rear spar that the 
crack extension modes varied from a tension failure at the most inboard end, which 
went completely through the surface and through the thickness of the web, to the begin- 
ning of a shear failure, some of which initiated inter lamlnarly. The propagation 
of the flaws from the most inboard to the most outboard cracks shown was verified 
by the strain gage data, which showed the inboard gages going nonlinear earlier than 
the outboard gages. 


Failure Description 


STABILIZER STABILIZER STABILIZER STABILIZER 
STATION STATION STATION STATION 



• WEB CRACKS, 
FORWARD FACE 



• THROUGH THICKNESS 
CRACK 


• WEB CRACKS, 
AFT FACE 
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This figure indicates the loading which initiated the failure just below the 
fail-safe lug at the most inboard location. This failure propagated aft, creating a 
shear concentration in the web which then propagated in a shear mode, causing the 
damage seen in the previous photographs. This tension load between the fail-safe 
lug and the lower lug is caused by the dihedral change at the side of the body. 

There was an aluminum angle spanning the lug areas; however, during the design, some 
of the lug area of the angle was trimmed away for clearance purposes, allowing a 
greater imposed deformation to occur at this location. Since the web must perform 
with a compatible strain, this tension strain then initiated the failure. 


Damage Tolerance 

Failure Sequence 
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This photograph illustrates the detail of the tension failure. The lug areas 
are made of co-cured highly uniaxially oriented fiber slabs which are then wrapped in 
two channel halves with cloth to form the "I" section of the spar. The two halves 
and the cap strip are then bonded together in a single envelope bagging operation. 

The failure occurred where the wrapped plies of cloth turned the corner and there- 
fore took the tension load in the resin direction, as previously described. Only 
XX plies were continuous straight up the center of the web. Therefore, this 
compatible deformation initiated the failure. 




TENSILE FAILURE OF INBOARD W^B 


DELAMINATION OF FACE PLIES 
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It can be seen, as was previously mentioned, that the most Inboard gage was 
nonlinear at a load below which the outboard gage was still linear. Indicating that 
the failure was beginning to propagate from the Inboard to the outboard direction. 


DAMAGE TOLERANCE- 
TEST FAILURE STRAIN GAGE READINGS 
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The accuracy of the finite element analysis as compared to strain gage data is 
shown in this figure. The extremely high shear strains that occurred during this 
test are illustrated. The additional shear strain Induced by the tension failure 
would raise the local shear strain in the web above the capability established from 
test data. 


Damage Tolerance 

Finite Element and Strain Gage Comparison 
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After establishing the failure mechanism, a plan for repair and avoidance of 
this failure mode was established. The conclusions shown here played a critical 
role in establishing the repair process. 


Damage Tolerance 

Conclusions 


• Load level 

• 61% as tested 

• 67% required 

• Rear-spar web detail insufficient for fail-safe 
loading 

• This area only critical area of stabilizer for 
this condition 


52 



Only the areas of the horizontal stabilizer that are critical for this condition 
are shaded in this figure. 

Only Area of Structure Critical for 


Upper-Pin-Removed Fail-Safe Case 




REAR 

SPAR 


LOWER 

SURFACE 


Repair Plan 


• Design a reinforcement for the rear-spar web area 


• Verify adequacy of rear spar by analysis 


• Apply modification to all five shipsets 
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This photograph shows the steel repair plate on the rear spar. The steel was 
used to assure minimum tension elongation and to minimize the thickness in the fit-up 
areas . 
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In considering what Is needed In the future to support the analysis and design 
of composite wings and bodies for commercial transports, we must consider the pos- 
sible criteria requirements as well as the design goals. To Improve the weight 
savings and to sustain the kind of damage that can be expected In service. It Is 
believed that an Increased design strain capability Is needed, both In design and In 
material application. We will have to fulfill category #2 and possibly category #3 
of the damage tolerance requirements If the raising of the strain level causes flaw 
growth. Obviously, In moving to the wing, high end loads will be encountered. As an 
example, these high end loads will require careful design at the side of the body In 
order to remove the high concentrated loads from the stringers. This design analysis 
and test matrix will be particularly critical In the success of wing design. 
Similarly, diffusion Into large fuselage structural shells of concentrated loads 
(e.g., the keel beams) will be extremely Important, and careful analysis will be re- 
quired. The determination of buckling criteria Is essential to establishing a good 
weight savings capability for composites In fuselage shells. Good criteria and test 
data, as well as validating analysis, are particularly Important. We must also 
consider the possibility that new uses will require a wider range of environmental 
conditions. Finally, raising the strain level in these critical structure components 
may again require us to consider fatigue damage. We need to have a matrix fatigue 
damage rule. 


WING & BODY - FUTURE ANALYSIS CONSIDERATIONS 

POSSIBLE CRITERIA, REQUIREMENTS OR DESIGN GOALS THAT WILL 
INFLUENCE ANALYSIS NEEDS 

INCREASE IN DESIGN STRAINS 

DAMAGE TOLERANCE REQUIREMENTS 
CATEGORIES M2 & #3 

HIGH LOAD TRANSFER JOINTS 

SIDE OF BODY JOINTS 
KEEL BEAM 

BUCKLING CRITERIA (AT OR BELOW LIMIT?) 

HIGHER TEMPERATURE (/MOISTURE) REQUIREMENTS 
BODY (WITH INSULATION) 

FATIGUE CONSIDERATIONS 
MATRIX DAMAGE RULE 
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There are three or four areas that need to be considered in developing 
the analytical capabilities to be used in these major components of primary struc- 
ture. They sound fairly simple, particularly the ones associated with strength, but 
since we will be asked to design closer to the capability of the material, the 
ability to predict both the ultimate strength and the residual strength requirements 
accurately in the design is much more important for both the static and damage 
tolerance conditions. Therefore, high competence in both aspects of the analysis 
and in validation by test is of extreme importance. Interlaminar effects, working 
the material in its weak direction, are of extreme concern to everyone. We must be 
able to predict interlaminar flaw growth, first of all from a grow-no-grow criteria 
and then under cyclic load. Those loads which induce these normal load stresses or 
strains are often secondary kick loads induced by eccentricities or secondary load 
paths. In metal structures these loads are often not a major concern, except 
possibly in fatigue. The secondary load calculations must improve in accuracy 
since they will most likely be the initiating phenomena for interlaminar flaw 
growth and failure. 


WING S BODY - ANALYSIS DEVELOPMENT AREAS 

DEVELOPMENT MUST CONSIDER REQUIREMENTS OF BOTH THE 
REQUIREMENT ANALYSIS AND CAPABILITIES ANALYSIS 

1. STRENGTH PREDICTION 

ULTIMATE STRENGTH 
RESIDUAL STRENGTH 
HIGH CONFIDENCE 

2. INTERLAMINAR EFFECT 

STATIC STRENGTH 
GROWTH PREDICTION 

5. SECONDARY LOADING 
NORMAL TO LAMINATE 
BUCKLING INDUCED 
PRESSURE INDUCED 

ECCENTRICITY OR KICK LOADS INDUCED 
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It Is Important to remember that we must bring these analytical tools along as 
early as possible, not necessarily In the most perfected form, but In a form that 
can be used as we begin to make these early configuration trade-offs, material 
selections, and manufacturing decisions. Once a company has begun to change from 
those pieces of equipment currently available for machining and assembling aluminum 
to those required for fabricating and assembling advanced composite structure, a 
major commitment will have been made. Therefore, It Is Important that we consider 
the development of methods that will support us In the near-term development of these 
components prior to production commitment. Development and validation of these 
methods must go hand in hand, and test techniques and procedures must be able to 
produce consistent and valid results. In addition, those test results must be 
validated, along with the analysis on structures simulating real aircraft loading. 


WING & BODY - ANALYSIS DEVELOPMENT CONSIDERATIONS 


TIMING 


ROLE IN EARLY CONFIGURATION SELECTION 
IMPACT ON MATERIAL DEVELOPMENT & SELECTION 
PLANNING OF DEVELOPMENT PROGRAMS & TESTING 

VALIDATION 

SUPPORT METHOD WITH ADEQUATE TESTING 
TESTING INCLUDES REAL STRUCTURE APPLICATIONS 


Here some thoughts are presented on what structural analytical goals should be 
In terms of the type of tools needed. 


STRUCTURAL ANALYTICAL TOOLS GOALS 


• ESTABLISH A LEVEL OF CONFIDENCE IN THE ANALYTICAL PROCESS 

SUCH THAT IT IS THE PRIME TOOL FOR STRUCTURAL VALIDATION 

- QUALITY SUCH THAT ANY REQUIRED FULL-SCALE TESTING WOULD 
ONLY BE USED TO PROVIDE FOR AIRCRAFT GROWTH OR TO UNCOVER 
GROSS HUMAN ERROR IN APPLICATION OF THE ANALYSIS OR 
MANUFACTURING METHODS 

- CAPABILITY OF BEING IMPLEMENTED IN A SIMPLE AND ECONOMIC 
MANNER BY A LARGE GROUP OF ENGINEERS 

- USABLE WITHIN THE CONSTRAINTS IMPOSED BY AVAILABLE DATA 
AND FLEXIBLE ENOUGH TO IMPROVE WITH THE EXPANSION OF 
AVAILABLE DATA 

- CAPABLE OF IMPACTING THE DESIGN AS WELL AS BEING USED TO 
ANALYZE THE DESIGN 

- VALIDATABLE BY TEST EVIDENCE 
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Here are some thoughts on the analysis development aspects, whether related to 
advanced composites or to other types of structures. For successful development, I 
believe that those people Involved In research need a group working environment, and 
that working environment Includes enthusiastic management backing, adequate but not 
excessive budgets, facilities available for developing the analytical tools in terms 
of computer access, test facilities, and, of course, supporting organizations. But 
one of the key elements for success Is easy Interface with production analysts and 
design people to ensure that the methods are going to be truly usable by the produc- 
tion engineers. 

Making a good Initial choice of which development areas to start Is extremely 
difficult. In making prioritized lists we must make sure that we are looking for 
real needs and that we plan our program to meet those needs when they are required 
at the right depth. Let's not produce a program of perfection that is too late to 
support some of the earlier decisions in the development of advanced composite struc- 
ture. We must be aware that our development should show true progress, not rehash 
over and over again the methods already available. On the other hand, let's not be 
afraid to adapt current methodology and carry that information from past experience. 

Whether the method Is used In a timely manner will depend on how familiar It Is 
in form to something the current production stress analyst is familiar with. This 
takes an honest and realistic evaluation of what is available so as not to produce a 
replacement for something that Is adequate today. Are all the items that are needed 
to peripherally support the analytical method well defined? Finally, take the time 
to establish a well-developed method specification and review this specification with 
the potential users, so the end product will be adequately suited to the job. 

ANALYSIS DEVELOPMENT 

ITEMS NEEDED FOR SUCCESSFUL DEVELOPMENT 

• GOOD WORKING ENVIRONMENT 

ENTHUSIASTIC BACKING 
ADEQUATE BUDGET 
FACILITIES 

SUPPORTIVE ORGANIZATION 

EASY INTERFACE WITH METHODS USER 

• GOOD INITIAL CHOICE OF DEVELOPMENT AREA 

REAL NEED CLEARLY IDENTIFIED 
PLAN TO MEET NEED (APPROPRIATELY) 

WHEN NEEDED 
RIGHT DEPTH 

UTILIZES CARRY-OVER FROM PAST EXPERIENCE 
TRUE DEVELOPMENT PROGRESS 
OR 

MODIFICATION TO CURRENT METHODS 
ALTERNATE METHODS EVALUATED - (HONESTLY) 

REPLACEMENT 

INFORMATION AVAILABLE TO SUPPORT SOLUTION 
ESTABLISHED AND WELL DEVELOPED METHOD SPECIFICATION 

• END PRODUCT QUALITY SUITED TO JOB 
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In order to produce methods of high quality, keep the process simple. In 
concept and in method, this produces a tremendous payoff in terms of development 
time, cost, and ease of utilization, and certainly helps assure the early success of 
both development and application. Many who are in research for a long period of time 
tend to become perfectionists in method development. Let's be careful of this trap 
while making sure that each , item we add to the method is a true Improvement and 
expands the information produced by the analysis procedure. 

The key to a manager's good over-all development and management of method 
development is the balance of the budget he makes available relative to the quality 
level of the method developed. To expand on the idea of a method's usability and 
simplicity, it also should show adaptability to being part of other methods, in order 
to be available for future expansion. 


ANALYSIS DEVELOPMENT 


CONSIDERATION DEVELOPING QUALITY METHOD 

• SIMPLE IN CONCEPT & METHOD 

MINIMAL DEVELOPMENT TIME & COST 

USAGE & APPLICATION EASIER 

INHERENT SIMPLICITY - HELPS ASSURE EASY SUCCESS 

• REFINEMENT OF METHOD 

PERFECTIONISM 
TRUE IMPROVEMENT 
EXPANDS INFORMATION 
BUDGET VS. QUALITY LEVEL 

• ADAPTABILITY 

USABLE AS PART OF OTHER METHOD 
BALANCED FOR CURRENT USAGE 
& 

FUTURE EXPANDED APPLICATION 
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COMPROMISE BETWEEN SIMPLICITY AND THEORETICAL PERFECTION 
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METHOD 


SIMPLE 

METHOD 
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AS THE ENGINEER SEES HIS METHOD 


AS IT REALLY IS 
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Accuracy in method development is very important, and this accuracy should be 
validated by tests. To produce analytical tools for which there is no test valida- 
tion in the real world is useless, and one or two data points are not sufficient if 
you expect production stress analysis personnel to accept your methodology. Testing 
should be repeated again and again and should be validated with test data. In 
establishing a method, be sure that there are no additional unnecessary additions to 
the solution which really do not improve the accuracy or produce more Information, as 
previously noted. Don't join in analysis fads; make a good, unbiased judgement of 
the need for the method. Consider design handbook methods for end products to be as 
important as a computer method. The majority of detail stress analysis is still 
performed by hand today. Finally, keep in mind that the production analyst is your 
final user. If analysis methods of similar nature can be joined to form a single 
standard procedure, do it. Provide the flexibility for all uses and users. Can it 
be used in a disciplined manner by large organizations? This is critical. 


ANALYSIS DEVELOPMENT 


• ACCURACY 

IMPROVED SOLUTION (OVER CURRENT) 

NO UNNECESSARY ADDITION 

• FADS 

METHOD OF IMPLEMENTATION 

COMPUTER VS. DESIGN MANUAL (FORMATS) 
NEW LOOK TO OLD SOLUTION 

PREJUDGEMENT OR PREJUDICED APPROACH 

POOR PROGRESSIVE JUDGEMENT 
EMOTION & SNAP JUDGEMENT 
(PRESSURE - TIMING) 

0 PRODUCTION ANALYSIS & STANDARDIZATION 

FLEXIBILITY FOR ALL USES & USERS 
AVAILABLE TO ALL PROJECTS 
DISCIPLINED USE BY LARGE ORGANIZATION 
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Proper and timely application requires a follow-up; this means provide support 
when it is needed, keep the user up to date on any modifications, listen to his 
Inputs, don't be defensive. Think about expanding the method's usage and adaptability 
to new applications, and continue to evaluate it against new methods and recognize 
when it must be replaced. Then, when it must be replaced, proceed with a new develop- 
ment as needed. 


ANALYSIS DEVELOPMENT 

PROVIDE ADEQUATE FOLLOW-UP 

CONTINUE INTERFACE WITH USER 

PROVIDE SUPPORT WHEN NEEDED 

KEEP USER UP TO DATE ON ANY MODIFICATION 
USER MANUAL (REVISIONS) 

LISTEN TO USER INPUTS (ACCEPT INPUTS) 

EXPANDED USAGE (SCOPE & TIME) 

CONTINUOUS REVIEW FOR ADAPTABILITY 
OR 

NEW APPLICATION 

CONTINUE EVALUATION AGAINST NEW METHOD 
(RECOGNIZE REPLACEMENT) 

PROCEED WITH NEW DEVELOPMENT WHEN NEW REAL NEED IDENTIFIED 
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Finally, here are some rules for success in the analysis development: 


(1) Coordinate carefully between the researcher and user; be sure that there 

are no other alternate methods before you start the development, and 
have alternate methods available and assessed should your development not 
produce results . 

(2) Allow time for corrections and updates, particularly after you expose the 

method to the user; be sure, as you plan your development, that you plan 
the means of validation, since you want to have high confidence in its 
results . 

(3) It is a good idea not only to have those specifications developed and 

reviewed by the user, but to have him on board in frequent discussions 
during the development . 

This will give the program timely evaluation, build confidence in its usage, and 
also ensure that it is used correctly so that you are not blamed for its error or 
inaccuracy. Remember, many methods operate in a production environment for a long 
period of time. They must be effective over this total period of time. 


ANALYSIS DEVELOPMENT 


DEVELOPMENT RULES FOR SUCCESS 

• COORDINATION BETWEEN 

DEVELOPER & USER 

t BACKUP OR ALTERNATE METHOD 
AVAILABILITY ASSESSED 

t TIME ALLOWED FOR CORRECTION & UPDATE 
AFTER USER EXPOSURE 

• ESTABLISH OR ESTABLISHED MEANS OF VALIDATION 

INSURANCE OF HIGH CONFIDENCE IN METHOD 

• DEVELOPER/USER BOTH RETAIN PESSIMISM 

• USER ON BOARD FROM START OF DEVELOPMENT 

TIMELY EVALUATION 
BUILDS CONFIDENCE IN USAGE 
ASSURE CORRECT USAGE 

METHOD OPERATES EFFECTIVELY IN PRODUCTION ENVIRONMENT 
FOR A LONG PERIOD OF TIME 
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SUMMARY 


f ANALYSIS PLAYS A MAJOR ROLE 

IN CERTIFICATION & DEVELOPMENT 

t ANALYSIS APPLICATION NEED TO BE 
APPROPRIATE TO JOB 

• ANALYSIS METHODS NEED TO BE 

TEST VALIDATED 

• ANALYSIS FUTURE NEEDS ARE: 

PREDICTION ACCURACY 
INTERLAMINAR STRESSES 
SECONDARY LOADS 

• ANALYSIS DEVELOPMENT MUST RECOGNIZE: 

REAL NEEDS 

ADVANTAGES OF SIMPLICITY 


65 


REFERENCES 


1. Airworthiness Standards: Transport Category, Airplanes. FAR Ft. 25, FAA, June 

1974. 

2. Composite Aircraft Structures. Report No. AC 20-107, FAA, July 1978. 


66 



A PRELIMINARY DAMAGE TOLERANCE METHODOLOGY 
FOR COMPOSITE STRUCTURES 


D. J. Wilkins 
General Dynamics 
Fort Worth, Texas 


PREVIEW 


The work described In this presentation was supported by the F-16 Airplane 
System Program Office, Aeronautical Systems Division, United States Air Force, as 
part of the F-16 Production Fleet Kanagement Program. As indicated In Figure 1, 
the presentation will briefly explore General Dynamics' certification experience 
for the primary, safety-of-f light composite structure applications on the F-16, 
will describe the rationale for the selection of delamlnatlon as the major issue 
for damage tolerance, and will outline the modeling approach selected. The 
development of the necessary coupon-level data base will be referred to and briefly 
summarized. The major emphasis will be the description of a full-scale test where 
delamlnatlon growth was obtained with which to demonstrate the validity of the 
selected approach. A summary will be used to review the generic features of the 
methodology. 


• CERTIFICATION EXPERIENCE 

• DELAMINATION RATIONALE 

• MODELING APPROACH 

• DATA BASE DEVELOPMENT 

• FULL-SCALE VERIFICATION 

• SUMMARY 

Figure 1 
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DEFINE THE WORDS 


The Air Force, in MIL-STD-1530A, Airplane Structural Integrity Program, 
sets forth specific definitions of the words "durability" and "damage tolerance. 
For the purpose of this presentation, the definitions are given in Figure 2. 


DURABILITY. THE ABILITY OF THE AIRFRAME TO RESIST 

CRACKING, CORROSION, THERMAL DEGRADA- 
TION, DELAMINATION, WEAR, AND THE EF- 
FECTS OF FOREIGN OBJECT DAMAGE FOR 
A SPECIFIED PERIOD OF TIME 


DA.MAGt: TOLERASCE THE ABILITY OF THE AIRFRAME TO RESIST 

FA ILURE DUE TO THE PRESENCE OF FLAWS, 
CRACKS, OR OTHER DAMAGE FOR A SPECIFIED 
PERIOD OF TIME 


Figure 2 



CERTIFY COMPOSITES BY TEST 


In contrast to the combined analysis/test procedures used to certify metallic 
parts on the F-16, the composite primary structure parts had to be certified 
completely by test (Figure 3) . Full-scale room temperature, dry static and fatigue 
tests of the vertical tail and the horizontal tail were supplemented with beam 
tests representing critical load paths. Beams were tested statically at environ- 
mental conditions ranging from cold, dry to hot, wet. Coupled load/environment 
fatigue tests were also performed. All of these tests were necessarily defined 
only in terms of the contractually specified "design usage." 

For reference, the vertical tail employs thick graphite-epoxy skins bolted to 
aluminum spars and ribs, and the horizontal tail uses graphite-epoxy skins bonded 
to full-depth aluminum honeycomb core and bolted to a titanium hub/spar sub- 
structure. 



Figure 3 
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COMPOSITES ANALYSIS IS NECESSARY 
FOR TRACKING USAGE CHANGES 


As we moved into Fleet Management activities, it became apparent that the 
design-usage-related test program was less than adequate for meeting the needs of 
the Air Force Logistics Command in terms of making maintenance decisions. And 
tests of each individual airplane /damage situation were out of the question. It 
was quickly concluded that composites analysis was necessary for tracking fleet 
and individual aircraft usage changes and to make repair decisions (Figure 4). 
Because few "zero-defect" parts are really built, the analysis must be able to 
relate defect size to time in service. The analysis should address the typical 
flaws occurring in manufacturing (durability flaws), those resulting from unusually 
bad manufacturing errors (damage tolerance flaws) , and large damages occurring as 
a result of service usage. The issue is how these various flaws behave as a 
function of severe, normal, or mild usage of the airplane. The indicated 
relationship of flaw size to a buckling limit is illustrated in the following 
figure. 


DEFECT 

SIZE 



Figure 4 


71 





DELAMINATION IS MOST PREVALENT GROWTH MODE 


In research work previously performed for the Air Force (Reference 1) , we 
discovered that if any life-limiting, growth-related failure mode was likely to 
happen in a composite structure, it is probably delamination. Figure 5 shows about 
15 inches of the overall 88-inch length of a box beam tested with a severe combina- 
tion of bending and torsion. A delamination propagated in the compression surface 
skin until the surface finally buckled and caused the beam to fail. Some insight 
into the nature of delamination in graphite-epoxy can be obtained by exploring the 
toughness of some typical laminates in each of the three directions, as is outlined 
on the next chart. 


Figure 5 



GRAPHITE-EPOXY FRACTURE TOUGHNESS 


Figure 6, adapted from the results of Reference 2, plots fracture toughness 
(Kjq as measured by edge-notched beams or center-notched tensile coupons) as a 
function of the percentage of 0° plies in a balanced, symmetric laminate containing 
only 0° and ±45° plies of unidirectional tape. The trends are typical of the 
Thornel 300/Narmco 5208 class of graphite-epoxy, but the absolute numbers are not 
strictly valid because no ASTM standards have been developed for fracture toughness 
of composites. Nevertheless, we would agree that a ±45° laminate has a reasonable 
toughness slightly less than alumlniom. As 45° plies are replaced by 0° plies, the 
toughness would be expected to increase until it was about doubled. On the other 
hand, as 45° plies are replaced by 90° plies, the toughness is expected to decrease 
until it reaches a minimum represented by a crack running completely in the resin 
between the fibers. The real point is that in the Interlaminar mode, the toughness 
is always at the minimum. 

These low values of interlaminar toughness are only significant as they relate 
to the out-of-plane loads discussed in the next chart. 



Figure 6 
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SOURCES OF OUT-OF-PLANE LOADS 


Shown schematically in Figure 7 are five of the most common design details 
from which real hardware is constructed. Each of these details, even under in- 
plane loading, gives rise to interlaminar normal and shear forces. The interlami- 
nar forces act directly on the plane of minimum toughness. As a result, 
delaminatlons may form, and if the in-plane forces are compressive, buckling 
becomes a pertinent issue. 


FREE NOTCH PLY BOND BOLTED 

EDGE (HOLE) DROP JOINT JOINT 
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TRANSFORM THE PROBLEM 


Creation of a compression coupon data base Is Impractical because Interlaminar 
tension and shear interact with local and overall buckling so that specimen size, 
end conditions, and lateral support make scaleup difficult. A simpler, more direct 
approach (Figure 8) Is to use an analysis to define the transfer function between 
applied loads and local rupture forces. Then the experimental effort can focus on 
tension and :shear failure mechanisms. 



ANALYSIS 

P 

•^(TENSION) 

■H’K 

^ISHEAR) 

~r~ 

o;r= F(P) 


TENSION TEST 



SHEAR TEST 



Figure 8 
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DELAMINATION ANALYSIS USES FRACTURE MECHANICS APPROACH 


The capability is available from adhesive fracture technology (Reference 3) 
for computing the profiles of strain-energy release rate around a delamination in 
a laminated plate under general loading. As schematically depicted by Figure 9, 
the profiles can be broken into three components for use in a fracture-mechanics- 
based failure theory (such as Reference 4). The calculations can be performed with 
finite element methods coupled with the virtual crack closure technique (Reference 
5). 



Figure 9 
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COMPOSITES DAMAGE TOLERANCE APPROACH 


With this background, an approach (Figure 10) was taken for damage tolerance 
of composites on the F-16. Achievable goals were carefully set for this study by 
emphasizing manufacturing defects (such as planar voids) rather than service 
defects, concentrating on delamination as the major mechanism, developing a fracture- 
mechanics-based model, and doing it in five increasingly complex phases. The first 
three phases defined acceptable coupon test specimens, checked the influence of 
environmental effects, and sought a method for predicting unidirectional delamina- 
tion growth under spectrum loads. The fourth phase was intended to extend the test 
base to limited two-dimensional growth in beam specimens representing thin and 
thick areas of the horizontal tall. Its results are not included in this presenta- 
tion. Correlating delamination growth in a full-scale horizontal tail to that 
predicted by the coupons was the goal of the final phase. 


• EMPHASIZE MANUFACTURING DEFECTS 

• CONCENTRATE ON DELAMINATION 

• DEVELOP A FRACTURE MECHANICS-BASED MODEL 

• COMPLETE A 5-PHASE PROGRAM 


• Coupon Screening 

• Environmental Effects 


• Beam Model 
Verification 


• Tail Model 
Verification 



Figure 10 
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COUPON TESTS 


The coupon test phases (Figure 11) were performed with a Mode I coupon and a 
mixed mode coupon with about 75% Mode II response. As reported in Reference 6, 
static tests produced the critical strain-energy release rates and their environ- 
mental behavior, and fatigue tests measured the delamination growth rates for 
constant amplitude and spectrum loading. 


MODE I 


MODE II 




Figure 11 
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COUPON SUMMARY 


Figure 12 outlines some of the important results of the coupon phases, 
results have since been supplemented by the research reported in Reference 


• DEVELOPED UNIQUE TEST SPECIMENS AND METHODS 

• MEASURED CRITICAL STRAIN-ENERGY RELEASE RATES 

• FOUND GROWTH BEHAVIOR TO OBEY SIMPLE POWER LAWS 

• CHECKED ENVIRONMENTAL EFFECTS 

• PREDICTED SPECTRUM GROWTH ACCURATELY 

Figure 12 


These 
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TAIL MODEL VERIFICATION 


An introduction for the final demonstration phase of the program is presented 
in Figure 13. 


• OBTAIN A TAIL WITH KNOWN DEFECTS 

• PERFORM FATIGUE TEST USING TRUNCATED SPECTRUM 
•DO FRACTURE ANALYSIS 

• CORRELATE MODEL TO GROWTH OBSERVATIONS 



Figure 13 



DEFECT LOCATIONS 


Horizontal tail unit number 110 was selected for test and was brought back 
after 144 flight hours of normal service. One reason for its selection was that 
it contained a "comer defect," a 2.1- inch long by 0. 3-inch wide delamination 
between the 51st and 52nd plies of the 52-ply graphite-epoxy skin, that had 
been repaired by injection of resin. Figure 14 shows the locations of the corner 
defect, the titanium hub-spar substructure, and the boundary of a ply termination 
separating the 12- and lO-ply thicknesses in the skin. In an area of high stress 
coinciding with the ply termination, two additional defects were introduced. One 
was an impact damage produced by dropping a 3.25-pound steel dart with a 1-lnch- 
diameter hemispherical nose from a height of 2 feet. The other was a routed 
3- inch hole in the skin down to the aluminum core, representing the clean-up hole 
around a typical skin repair that never got a patch bonded over it. 


-5 SKIN SIDE 


-6 SKIN SIDE 




Figure 14 
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TAIL TEST HISTORY 


The tail was tested In an air-conditioned lab environment. Figure 15 outlines 
the testing sequence. 


• 1 LIFETIME OF DESIGN USAGE SPECTRUM LOADING 

• 2 SPECTRUM REPETITIONS WITH LOADS MAGNIFIED BY 1.3 

• 4162 FULLY REVERSED CYCLES OF 1 .3 TIMES MAXIMUM 
SPECTRUM LOAD 

• RESIDUAL STRENGTH TEST TO 246% DESIGN LIMIT LOAD 

Figure 15 
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TAIL FLAW GROWTH 


The overall views of both tail surfaces (Figure 16) show the outlines of 
delaminations created by the testing. The impact damage grew only slightly 
during the first few loads and exhibited no further growth. Delaminations grew 
from the 3- inch hole during the constant amplitude cycling. Large delaminations 
were produced in the root area of both surfaces. The first large increment of 
delamination at the 3-inch hole and at the root occurred during the first 2235 
cycles of the constant amplitude cycling during which no inspections were made. 



CORNER DEFECT 

• SLIGHT GROWTH DURING THIRD LIFETIME 

• MASSIVE GROWTH DURING CONSTANT AMPLITUDE CYCLING 
•STABLE PROPAGATION 


Figure 16 
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ROOT AREA MATH MODELS 


The finite-element based fracture analysis used the "zoom-in" technique 
shown in Figure 17. Each succeeding detailed model was loaded with deflections 
from the previous model using MSG NASTRAN. 



THREE DIMENSIONAL 
FINE GRID 


Figure 17 
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PLY CRACKING APPARENTLY TRIGGERED DELAMINATION 


Figure 18 shows the 3-D brick model of the skin, substructure, and core in 
the local area of the "corner defect," as well as its top view and an internal 
view of the local geometry. The model showed a strong through- the- thickness 
stress concentration associated with the bonding of the skin to the titanium 
spar. Transverse stresses in the bottom +45° ply were much higher than the 
ply could carry, so its stiffnesses were modified to represent a cracked ply. 
Introduction of a delamination by releasing nodes as indicated gave a strain- 
energy release rate that was 90% of the critical static strain-energy release 
rate and thus was near the upper end of the delamination growth curve referred 
to in Figure 11. 
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DELAMINATION GREW FROM 3- INCH HOLE 


A closer view of the impact damage and 3-inch hole growth indications is 
presented in Figure 19, along with the number of fully reversed constant 
amplitude cycles associated with them. 



Figure 19 
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3-INCH HOLE MATH MODELS 


Figure 20 depicts the "zoom-in" technique used for the 3-lnch hole math 
models. 



THREE DIMENSIONAL 
FINE GRID 


Figure 20 



SEQUENCE OF CRACK ANALYSIS 


The edge of the 3- inch hole is shown in Figure 21, along with the stacking 
sequence of the plies , all of which were modeled individually in the 3-D brick 
model. The model delamlnation was located between a -45° ply and a 0° ply. A 
node was released from the edge, then in all three possible directions, to 
indicate the preferred growth direction, which was along the ply termination. 
Succeeding nodes were released to "grow" the delamination. 



Figure 21 
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FRACTURE ANALYSIS AT P^x LOAD 

Results of the models, scaled to the peak load In the original design usage 
spectrum, are plotted in Figure 22. The analysis showed that a delamination of 
0.085 inches should have resulted from that load. The plot of strain-energy 
release rate as a function of crack length was scaled to the square of the applied 
load and combined with the coupon static fracture and delamination growth rate 
data to produce the results shown in Figure 23. 

DAMAGE LOCATIONS EDGE OF HOLE LOOKING DOWN 



Figure 22 
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ANALYSIS AGREES WITH OBSERVED GROWTH 


Figure 23 sunmarlzes and correlates the growth calculations and observations. 
The resolution of the hand-held ultrasonic Inspection technique at the edge of the 
hole was not good enough to locate the small delamlnatlons predicted by the 
analysis. However, edge cracks were observed by surface replication of the cut 
surface of the hole. The agreement between calculations and measurements for the 
larger delamlnatlons Is quite satisfactory. 



CALCULATED MEASURED CALCULATED MEASURED 

CRACK CRACK CRACK CRACK 

FIRST LIFE SECOND & THIRD REPETITION 



Figure 23 



SUMMARY 


The highlights of the program are reviewed and summarized in Figure 2A. 


COUPONS 

• DEVELOPED UNIQUE TEST METHODS 

• FOUND SIGNIFICANT DIFFERENCES BETWEEN PEEL & SHEAR 

• PREDICTED SPECTRUM GROWTH 


TAIL 

• OBTAINED GROWTH ONLY UNDER EXTREME CONDITIONS 

• DEMONSTRATED HIGH DAMAGE TOLERANCE 

• SHOWED STABLE GROWTH IN REAL HARDWARE 

• DEVELOPED MODELING PROCEDURE 


Figure 24 


COMPOSITES DURABILITY AND DAMAGE TOLERANCE ANALYSIS 


The methodology, which General Dynamics believes to be equally applicable to 
both dnrability and damage tolerance analyses, is outlined in Figure 25. Its 
input is a flight -by-flight load history that is coupled to the fracture analysis 
to produce a sequence to strain-energy release rate/delamination size relationships. 
The coupon data base can then be entered to judge the potential for static or 
cyclic growth including the Influence of moisture and temperature. A simple 
load-by-load growth model produces predictions of delamlnatlon size as a function 
of flight hours. This methodology has many attractive features, but the most 
important one is that it should be applicable to virtually any laminated composite 
structure. 



Figure 25 
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DESIGN AND ANALYSIS OF COMPOSITE STRUCTURES 
WITH STRESS CONCENTRATIONS 


S. P. Garbo 

McDonnell Aircraft Co., McDonnell Douglas Corporation 
St. Louis, Missouri 



INTRODUCTION 


This presentation provides an overview of an analytic procedure which can be 
used to provide comprehensive stress and strength analysis of composite structures 
with stress concentrations. The methodology provides designer/ analysts with a user- 
oriented procedure which, within acceptable engineering accuracy, accounts for the 
effects of a wide range of application design variables. The procedure permits the 
strength of arbitrary laminate constructions under general bearing/bypass load 
conditions to be predicted with only unnotched unidirectional strength and stiffness 
input data required. 

In particular, the presentation includes: 

(1) A brief discussion of the relevancy of this analysis to the design of 
primary aircraft structure 

(2) An overview of the analytic procedure with theory/test correlations 

(3) An example of the use and interaction of this strength analysis relative 
to the design of high-load transfer bolted composite joints 

Many of the results presented were taken from two recently completed MCAIR 
programs funded by the U.S. Air Force and the U.S. Navy (refs. 1 and 2). 
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DESIGN/ANALYSIS ITEBATIONS 


The usual project sequence of design/analysis Iterations is Illustrated. The 
emphasis of methodology discussed in this presentation is in the areas of (1) 
detailed stress analysis about a structural discontinuity or stress riser which 
results in local stress concentrations and (2) strength analysis performed on the 
stress (strain) solutions to determine load capability. 


^ INTERNAL LOADS > 
AND/OR 

BOLT LOAD DISTRIBUTION 
V ANALYSIS y 


DESIGN MODIFICATIONS 

\ • GEOMETRY j 

\ • MATERIALS / 
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F/A-18A WING STRUCTURE 


The planview of the F/A-18A wing structure is presented to illustrate the 
prevalence of stress risers in typical primary aircraft structure. This production 
structure is con 5 >rised of monolithic composite skins mechanically fastened to 
aluminum spar and rib substructure. Evident are numerous stress risers due to 
required access doors, cutouts, and pylon post holes required for attachment of 
external equipment. Additionally, fastener holes are required along all spars and 
ribs to mechanically attach wlngsklns to substructure, and local load introductions 
(e.g., leading- and tralling-edge flap attachments) and major structural splices 
(e.g., wingfold) result in many required high-load transfer joints. The point to 
be made is that for current and near-term aircraft, accurate analysis of the effects 
of stress risers on structural efficiency is of prime importance. 
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STIFFNESS AND STRENGTH DESIGN DATA 

The items listed under lamina (unidirectional ply) and laminate Indicate 
variables which affect stiffness and strength. Fundamental differences between the 
mechanical behavior of metals and composites are due to the composite stiffness and 
strength anisotropy and its material Inhomogeneity. These factors complicate 
considerably the complete characterization of lamina mechanical properties relative 
to metals and introduce the necessity of additional strength characterization 
requirements at the laminate level for each layup variation. The objective of our 
analysis development has been to provide a comprehensive and general laminate 
strength analysis based primarily on lamina stiffness and strength data, thus 
minimizing the level of laminate testing currently required for applications. 


LAMINA LAMINATE 



• IN-PLANE 
LOADINGS . . 


TENSION 
^ COMPRESSION 
SHEAR 




IN-PLANE 
LOADINGS . . 


BIAXIAL 1 

off-axis; ■ ■ ■ 


TENSION 
COMPRESSION 
' SHEAR 
BEARING 


• INTERLAMINAR 
LOADINGS 

• ENVIRONMENT 

• NONLINEARITIES 


• INTERLAMINAR 

LOADINGS 

• MATERIAL TAILORING 


/pressure 

(BENDING 

LAYUP 

^ STACKING SEQUENCE 
ANISOTROPY 


• GEOMETRY... 

• ENVIRONMENT 


* 

HOLE SIZE 
- HOLE TOLERANCE 
FASTENER TYPE 

I. 


• NONLINEARITIES 
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AS/3501-6 GRAPHITE-EPOXY LAMINA PROPERTIES 


This table summarizes stiffness and strength data along the principal material 
axes of the AS/3501-6 graphite-epoxy material system. These data are usually 
generated by testing unnotched specimens and most data are routinely obtained during 
material quality assurance testing performed by companies at the time material 
shipments are received. These data are the basis upon which general laminate 
stiffness, stress, and strength analyses are performed using the procedures 
discussed In this presentation. 


MECHANICAL 

PROPERTY 

rtdA 

RTW^ 

ETW^ 

E-|t (MSI) 

18.85 

18.85 

18.54 

Ei= (MSI) 

18.20 

18.20 

17.80 

E2 (MSI) 

1.90 

1.63 

0.74 

Gi2 (MSI) 

0.85 

0.85 

0.38 


0.30 

0.30 

0.30 

Fi<“ (KSI) 

230 

230 

236 

F-|CU (KSI) 

321 

179 

111 

F 2 "' (KSI) 

9.5 

5.6 

4.3 

F2'“ (KSI) 

38.9 

33.4 

12.4 

Fi2 (KSI) 

17.3 

19.8 

11.9 


/K Room temperature dry (as manufactured) 

Room temperature with 1.05% moisture 
Elevated temperature ~ 250°F and 1 .05% moisture 



INTERNAL LOADS ANALYSIS 


While this is not the main subject of this presentation, it is important to 
realize that the needed detailed stress and strength analyses, especially for 
bolted composite joints, are directly related to Internal load distribution 
analysis. Should local laminate failures occur due to fastener loadings, joint- 
member and fastener-to-joint-member flexibilities (spring rates) will be altered 
and internal load paths will change. The variables listed are important to accurate 
internal load analysis and have direct effects on laminate point strength analysis . 


• GEOMETRY. . jIN-PLANE 

Ithru-the-thickness 

• MEMBER STIFFNESSES 

• FASTENER STIFFNESSES (SPRING RATES) 

• MANUFACTURING TOLERANCES 

• MULTIPLE FASTENERS 

• LOADING COMPLEXITIES 

TENSION 

COMPRESSION 

SHEAR 

VARIATIONS 

• MATERIAL NONLINEARITIES 


102 



DETAILED STRESS DISTRIBUTION ANALYSIS 


This list represents some of the more important variables typical of general 
eiircraft structure. The initial thrtast of our strength analysis development was to 
provide a user-oriented procedure which would account for the effects of fastener 
holes on laminate strength for arbitrary laminates under general biaxial and 
fastener bearing loads. Currently we are extending these procedures to account for 
the general class of stress risers, through-the-thlckness variables, and finite- 
geometry effects . 


• GENERAL LOADINGS 


• GEOMETRY 


• THROUGH-THE-THICKNESS 


• MATERIAL NONLINEARITIES 

• ENVIRONMENT 


BIAXIAL 

< BEARING/BYPASS 
BENDING/PRESSURE 


HOLE SIZE 
CUTOUT SHAPE 

i WIDTH (EXTERNAL BOUNDARIES) 
MULTIPLE STRESS RISERS 
.HOLE TOLERANCES 


r 




STACKING SEQUENCES 
CLAMP-UP (TORQUE-UP) 
JOINT ECCENTRICITY 
FASTENER TYPES (CSK) 
FOUNDATION STIFFNESS 


I^ORTHOTROPY 
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FAILURE PREDICTION OF BOLTED COMPOSITE JOINTS 


An overview of the MCAIR analytic approach Is Idealized In this Illustration. 
The principle of superposition Is used to obtain laminate stress and strain 
distributions due to combined bearing and bypass loads. Once laminate stresses and 
strains are known, failure analysis Is performed at a characteristic dimension (Rq) 
away from the hole. This bolted joint stress field model (BJSFM) generalizes 
Whitney and Nulsmer's postulated characteristic dimension concept to apply to 
anisotropic laminates under general biaxial loadings (ref. 3) . Failure predictions 
are made on a ply-by-ply basis by comparing material lamina allowables with elastic 
point stresses calculated for each ply at a characteristic dimension from the edge 
of the hole. Various material failure criteria (Tsai -Hill, maximum stress, etc.) 
available for user selection can then be used to compare stresses with lamina 
allowables . 


This program Is available to the Industry through the Air Force Flight Dynamics 
Laboratory and Is detailed In reference 1. This user-oriented program can be used 
to generate complete stress (strain) solutions at or about the hole, or just the 
failure load, critical ply, and failure location. 


( II t I I I I I I I I I I 



II ( I I I I I I I I I I I 
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CIRCUMFERENTIAL STRESS SOLUTIONS UNDER PURE BEARING LOAD 


The need for detailed distributions In orthotropic laminates Is Illustrated by 
using the BJSFM program to determine circumferential stress distributions about a 
fastener hole loaded by pure bearing. A 50/40/10 layup and a 70/20/10 layup from 
the 0°, ±45°, and 90° family of ply orientations were analyzed. Solutions for both 
of these layups are compared with solutions obtained for an Isotropic material. In 
the first two comparisons, bolt bearing was directed along the 0° ply orientation 
(principal material axis) . In comparison with the dashed line representing Isotropic 
material stress solutions, orthotropy In the two composite layups causes a marked 
change in maxlmiom stress and locations. This directional dependence Increases with 
the degree of orthotropy (higher percentage of 0° plies relative to load direction) , 
as can be observed by conq)arlng the solutions of the 70/20/10 layup with the 50/40/10 
layup. Additionally, If bolt bearing Is directed off the principal material axis, 
the circumferential stresses change markedly. The dashed line illustrates the fact 
that solutions for an Isotropic (metal) material are unaffected by change in bolt 
bearing direction. 



— Isotropic 
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EFFECTS OF FINITE GEOMETRY 


To evaluate the effects of finite geometry on stress solutions and establish 
limits on the accuracy of infinite plate solutions, a comparison of BJSFM solutions 
with finite element solutions was performed. In this illustration, circumferential 
stresses normalized to the average bolt bearing stress are plotted about the half- 
circle from directly in front of the neat-fit bolt (0°) to directly behind the bolt 
(180°). The dashed line indicates the infinite plate BJSFM solutions, the solid 
line represents the BJSFM solution corrected using the DeJong finite-width approxi- 
mation method, and the triangular symbols represent the finite element solution. 
Results indicated that for an e/D of 9, the BJSFM solutions, corrected for finite 
width, correlate extremely well with finite element solutions where W/D ratios were 
greater than 4. However, as indicated in this graph, at an e/D value of 3, BJSFM 
infinite-plate solutions that were approximately corrected for finite widths (long- 
dash line) , while significantly improved over the original BJSFM infinite-plate 
solution (short-dash line), still differ from the correct finite geometry solution 
(triangle symbol line) . Full correction for finite geometry is anticipated using 
boundary collocation procedures currently under development at MCAIR. 


SMALL EDGE DISTANCE 
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THROUGH-THE-THICKNESS ANALYSIS 


The previous results Ignore thickness-direction variables. These variables 
can be Important and additional procedures are required to account for their effects. 
Some of the related through-the-thickness variables are listed here. New procedures 
using beamr-on-elastic foundation theory are currently being developed at MCAIR to 
examine the effect of these design variables on fastener- to- joint-member spring rates 
and on bearing stress distributions in the thickness direction. Preliminary results 
of one parametric study are shown that relate the peaking of bearing stresses 
(relative to an average bearing distribution) to fastener head fixity. 



• JOINT ECCENTRICITY 

• INTERLAMINAR STRESSES 

• STACKING EFFECTS 

• THICKNESS EFFECTS 

• ORTHOTROPY EFFECTS 

• FASTENER FIT 

• INTER-PLY REINFORCEMENTS 


PEAKING 

FACTOR 



HEAD FIXITY 



STRENGTH ANALYSIS 


/ 

/ 
i 

/ 

Given a detailed ply-by-ply stress analysis, strength predictions using the 
characteristic dimension failure hypothesis can be made. Point material failure 
criteria (Tsai-Hill, maximum strain, etc.) are used to predict ply matrix, fiber, or 
shear failure about the fastener hole at a distance Rq away from the hole boundary. 

This is automatically performed by the BJSFM program by comparing the complete stress 
(strain) solution about the fastener hole boundary with user-selected material 
failure criteria. Critical plies, failure loads, and locations are identified. 

However, in addition to failure at a point, the "patterns" of failure must be 
evaluated to predict overall joint modes of failure such as bearing, shearout, or 
net section. Examples of these predictions and correlations with data are presented 
in the following illustrations. 


• LOCAL AND GLOBAL 


"load 

LOCATION 
CRITICAL PLIES^ 


. . FAILURE MODES 



MATRIX CRACKING 
FIBER RUPTURE 
SHEAR 

COMPRESSION 

INTERLAMINAR 


FIRST PLY 
PROGRESSIVE 


• GENERAL LOAD CONDITIONS 

• ENVIRONMENTAL FACTORS 

• MATERIAL NONLINEARITIES 
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BEARING/BYPASS STRENGTH PREDICTIONS 


Examples of laminate strength predictions under general bearing and bypass 
loadings are presented for a wide range of layups. The higher stress predictions 
of each of the solid and dashed line sets represent predictions at the characteristic 
dimension away from the hole boundary. The solid line represents predicted ply fiber 
or shear failure and correlates within ±10 percent with the test data for this wide 
range of layups. The dashed lines represent predicted matrix failure (i.e., matrix 
cracking due to tensile load transverse to fibers) but no correlation of this 
predicted failure with changes in mechanical properties was observed, and this type 
of ply failure was ignored in all subsequent analyses. The lower strength predic- 
tions of each of the solid and dashed line sets represent application of failure 
criteria on the hole boundary. The difference between predictions on the hole 
boundary (R^ = 0.0 inch) or away from the hole boundary (Rj. = 0.02 inch) correlates 
with the conservatism inherent in strictly using elastic stress concentration 
factors to predict strength. 


PREDICTED PLY MATRIX, FIBER, 
AND SHEAR FAILURES 





BYPASS 



■“ Ply fiber or shear failure 

“ ““ ““ Ply matrix failure 



0 20 40 60 80 100 
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HOLE SIZE AND LAYUP EFFECTS 


Hole size effects are accounted for analytically. Stress gradients at the edge 
of the hole vary with hole size and laminate orientation. Smaller diameter holes 
produce steeper stress gradients which decay rapidly as the distance from the hole Is 
Increased. This confines maximum stress to a smaller area In the laminate; thus a 
laminate has greater load-carrying capability with a small hole than with a larger 
hole. Analytical correlation is shown with strength data for various hole sizes and 
laminates loaded In tension and con^resslon. 


TENSION ^ 
Rc = 0.02 IN. 


COMPRESSION 
Rg = 0.025 



HOLE DIAMETER - IN. HOLE DIAMETER - IN. HOLE DIAMETER - IN. 
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OFF-AXIS TESTING 

Test data are compared with predictions for a laminate with a countersunk hole 
that Is loaded to failure unlaxlally at various angles relative to the principal 
material axis. Countersunk fastener holes were analyzed as being equivalent In 
diameter to a constant-diameter hole which displaces the same volume of material as 
the countersunk hole. The correlation shows that the BJSFM procedure can account 
for countersunk hole configurations, strength anisotropy, and off-axis loadings while 
using a constant characteristic dimension (Re) • 




REPRESENTATIVE LOAD-DEFLECTION DATA FOR 50/40/10 LAYUP 


The previous theory-test correlations presented strength data for laminates 
with an unloaded hole. As bearing stresses are applied to the hole boundary, 
bypass strength decreases and noncontlnuous or nonlinear load-deflection behavior 
Is observed prior to specimen ultimate failure. At low bearing stress conditions, 
strength predictions based on first-ply fiber failure correlate very closely with 
the ultimate specimen strength. However, as bearing loads Increase, first-ply fiber 
and shear failures correlate with the occurrence of Initial nonlinear or discon- 
tinuous load-deflection behavior. 
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BEARING-ONLY LOAD CONDITION 



The local nature of failure due to pure bearing loads Is examined analytically 
using the BJSFM procedure. The results shown indicate the critical plies and 
contours of equal strain levels for an 80-ksi bearing load when first ply failure Is 
predicted. Ply failures occur at 0° and 45° due to excessive shear. The contours 
indicate where failure patterns in specimens would likely occur. The C-scan 
documents the occurrence of local damage at this bearing stress level, and 
qualitatively agrees with predicted failure contours. 


SHEAR-CRITICAL 0 “ PLY SHEAR-CRITICAL - 45 “ PLY 



fgr = 80 KSI 
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BEAKING-PLUS-BYPASS LOAD CONDITION 

The analysis of the prevloxis pure bearing case was extended to include the 
added effect of a bypass load equivalent to a stress of 30 ksi. The critical plies 
and contours are indicated. The patterns have shifted to the net section and now 
the 0° ply has become fiber critical (before it was shear critical) ; the 45° ply 
is still shear critical, but this is now at the net-section region as well as in 
front of the bolt. The same specimen previously loaded to 80 ksi of bearing only 
was loaded to that same bearing load and then further loaded to a bypass stress of 
30 ksi. The C-scan again qualitatively verified predicted failure patterns. 


FIBER CRITICAL 0^ PLY SHEAR-CRITICAL -45^ PLY 
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OFF-AXIS BEARING TEST DATA 


In the of f-the-principal-material-ax±s study, the BJSFM procedure was used to 
examine the directional dependence of pure bearing strength in a 50/40/10 laminate . 
Test data with scatter bands are compared with predictions. Onset of joint 
nonlinear load deflection behavior is indicated by open triangle symbols. Specimen 
ultimate bearing strengths are indicated by open circles and occurred at approxi- 
mately 50 percent higher loads. Specimens failed predominantly in bearing shearout 
at off-axis load directions ranging from 0° to 22-1/2°, tension-cleavage at the 45° 
off-axis load direction, and net-section failures at the 67-1/2° to-90° off-axis 
load directions. Solid symbols indicate specimens tested with no fastener torque-up; 
all others were torqued to 50 in. -lb. Strengths were predicted based on ply fiber 
or shear failure using the 0.02-in. characteristic dimension and the maximum strain 
failure criterion. Compared to off-axis loading of unloaded hole specimens, two 
important characteristics are predicted: (1) laminate strength is relatively 

constant over the entire range of off-axis load directions, and (2) ±45° and 90° 
plies are equally critical within 10 percent. This is in contrast to the unloaded 
hole case in which a pronounced sensitivity to off-axis bypass loading was observed 
and singular ply orientations were usually critical. 

Unlike loads at onset of nonlinear joint load deflection behavior, test results 
indicate that ultimate bearing strengths are Insensitive to off-axis loading effects 
at all but the 0° and 90° loading directions. In all cases, specimen failures 
occurred after permanent hole elongations had developed. The onset of the 
associated transition from linear to nonlinear failure modes was predictable using 
the linear-elastic BJSFM procedure. 

(50/40/10 LAYUP) 


BEARING 

STRESS 

KSI 
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HIGH-LOAD LUG 


A further verification of the BJSFM strength procedure was performed using a 
boundary collocation procedure to account for finite geometry effects. The 1.50- 
in. -thick specimen was designed to contain two different lug sizes that would be 
loaded in pure bearing. The lower left photograph indicates the failed left end, 
and the lower right C-scan documents the close-to-failure damage of the right end 
of the specimen. Failure modes of both ends were as predicted. 


BEARING 


Note: Dimensions in inches 


PREDICTED 

SHEAR FAILURE 
154 kips 


PREDICTED 

NET-SECTION 
159 kips 


STRENGTH ENVELOPES FOR APPLICATIONS 


The utility of the previous analytic procedures was recently demonstrated in a 
MCAIR/NADC-funded program (ref. 2) to design high-load- transfer metal-to-composite 
bolted joints applicable to a fighter wing root area. This graph indicates predicted 
strength envelopes used to design the required joint for either a monolithic 55/44/4 
laminate (lower curve) or a locally softened 52/44/4 laminate (horizontal, upper 
curve) where 0° plies were replaced by ±45° plies local to fastener hole areas. For 
the monolithic case, the limit cutoff represents the predicted first-ply shear 
failures under pure bearing loads . 



BEARING • KSI 
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JOINT THICKNESS VERSUS NUMBER OF BOLTS 


Based on a design ultimate load requirement of 30,000 Ib/ln., the previous 
strength envelope for the monolithic laminate was used to predict the number of 
In-line bolts required versus overall joint thickness. Cutoffs based on fastener 
shear strength are also Indicated. 




LAYUP VERSUS OVERALL JOINT THICKNESS 


Based on a seven-bo It- In- tandem joint configuration, the BJSFM laminate strength 
procedure was used to efficiently evaluate the effect of layup variations on joint 
thickness. The results of this study are plotted for layups having 4 to 25 percent 
of their 90° plies aligned with the chord direction of the baseline wing; the 
remaining plies were oriented in the spanwlse (0^^) direction. 
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SCARFED MONOLITHIC-LAMINATE TEST SPECIMEN 


Based on the results of the previous two analytic studies, scarfed monolithic- 
laminate and scarfed sof tened-laminate specimens were fabricated to verify predic- 
tions. The graphite-epoxy and titanium members are shown with the 0 . 50- in. -diameter 
fasteners used in the monolithic-laminate design. The specimen is representative of 
a single column of a wing root splice that would be spaced every 2.50 in. across the 
chord of the wing root. 
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SCARFED SOFTENED-LAMINATE TEST SPECIMEN 


The scarfed sof tened-laminate splice is shown assembled and strain gaged prior 
to test. The sof tened-laminate splice test specimen represents a single spanwise 
column of the wing root design which would be repeated approximately every 10 in. 
across the chord of the baseline wing root area. The taper of the planform was 
designed to simulate the gradual increase of spanwise loading due to increased wing 
bending moments encountered closer to the root area. 
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REPRESENTATIVE STATIC FAILURES 


The scarfed monolithic-laminate splice (upper photo) and sof tened-laminate 
splice (lower photo) are pictured after static failure. As predicted, both specimens 
failed in the net-section area of that fastener location having the highest bypass 
load. The static joint failure load for the monolithic-laminate was 68 kips; the 
sof tened-laminate failed at 125 kips. Both loads were within 10 percent of 
prediction. 





EFFECT OF HOLE CLEARANCES ON BOLT LOADS 


To investigate the cause of scatter in bolt load distributions, the effect of a 
0.003-in. and 0.006-in. clearance at each of the fasteners was analyzed. The results 
illustrated in this figure offer a possible explanation of the scatter and indicate 
the importance of this variable. 




FASTENER NUMBER, N 
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EFFECT OF HOLE CLEARANCE ON LOCAL MEMBER LOADS 

The occurrence of a U. 003- In. clearance on each of the fasteners was evaluated 
relative to its effect on joint member loads. Unlike its effect on bolt loadings, 
when the area of the joint members is taken into account, the effect of clearance 
on joint member loading is relatively small. 


HOLE CLEARANCE = 0.003 IN. 
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CONCLUSIONS 


As a backdrop for conclusions, this sketch idealizes the evolution of structural 
testing typical of project hardware development programs. Usually testing proceeds 
from simple coupon specimens to more structurally detailed element specimens, then to 
subcomponent box beams, and finally to full-scale static and fatigue-test articles. 
For the application of composites the testing requirements are extensive and costly, 
especially for establishing material allowables, because of loading possibilities and 
geometry effects which change for each laminate variation. 

The BJSFM analysis procedure should considerably reduce the need for coupon 
testing relative to laminate variations, hole size, general bearing/bypass load 
interactions, and environmental test conditions. 

Analysis and experimental data were compared to verify that the BJSFM procedure 
accurately predicts laminate strength under general tension and compression bearing/ 
bypass loadings. Theoretical studies were reported which demonstrated method utility 
and the need for detailed stress solutions. Correlations between the closed-form 
BJSFM stress solutions and finite element solutions were reported which defined 
current program finite geometry limitations. 

Three benefits are provided by using the BJSFM procedure. (1) Test data 
requirements are minimized, requiring only knowledge of mechanical properties for 
the unidirectional-ply and one-laminate configuration. (2) Design conservatisms and 
unconservatisms can be reduced through this accurate, comprehensive, and inexpensive 
analysis, (3) parametric evaluations of failure mechanisms can easily be con- 
ducted to provide guidance for the design of efficient bolted composite joints. 



FULL SCALE STATIC AND FATIGUE 
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INTRODUCTION 


The safe-life and damage-tolerant design approaches presented herein 
apply to both metallic and fibrous composite helicopter structures. 
However, to maintain the theme of this workshop, this presentation will 
emphasize the application of these design approaches to fibrous compos- 
ite structures. 

In the current generation of Army helicopters, such as the UH-60 
Black Hawk shown in Figure 1, composite materials make up as high as 17 
percent of the airframe and rotor weight. With the advent of major 
helicopter composite structures R&D projects, such as the Advanced 
Composite Airframe Program (ACAP) , and Manufacturing Methods and Tech- 
nology (MM&T) projects, such as UH-60 Low Cost Composite Blade Program, 
it is estimated that within a few years composite materials could be 
applied to as much as 80% of the airframe and rotor weight of a heli- 
copter in a production program. Along with this application is the 
essential requirement that sound, definitive design criteria be developed 
in order that the composite structures possess high fatigue lives for 
economy of ownership and good damage tolerance for flight safety. 

Safe-life and damage-tolerant criteria are applied to all helicopter 
flight critical components, which are generally categorized as follows: 

Dynamic components - main and tail rotor system, which includes 
blades, hub and rotating controls, and drive train which includes trans- 
mission, and main and interconnecting rotor shafts. 

Airframe - fuselage, aerodynamic surfaces and landing gear. 



Figure 1 
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TYPICAL STRESS SPECTRA FOR A HELICOPTER BLADE 
ROOT END AND TRANSPORT AIRPLANE ROOT END 

The dynamic components previously described tend to be fatigue 
critical. Although fatigue problems have been known to occur within 
the helicopter airframe, it for the most part is designed by strength 
and minimum gage. 

Figure 2 shows a typical comparison of the stress history for a tirans 
port helicopter blade root end and that for a typical commercial trans- 
port wing root. The helicopter dynamic components tend to be loaded at 
frequencies that are a multiple of rotor speed - typically 2 to 20 Hz. 
This loading rate represents 7 million to 70 million cycles per 1000 
operating hours. Thus these components experience a much higher fre- 
quency of cyclic loading but with a narrower spread of the loading than 
the fixed wing aircraft component. 
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SAFE EIPB PHILOSOPHY 


The safe-life design approach is stated in Figure 3. It applies to 
statistically predictable failures which are those occurring from the 
random combination of fatigue strength and applied fatigue loads. 


ALL FATIGUE CRITIUL COMPONENTS OF 
THE HELICOPTER ARE DESIGNED TO A 
SPECIFIC SERVICE UFE, IN OPERATIONAL 
HOURS, AND ARE REMOVED FROM SERVICE 
AT OR BEFORE THIS ELAPSED TIME SO 
THAT THE PROBABILITY OF A FATIGUE 
FAILURE IS REMOTE 


Figure 3 



SAFE -LIFE METHODOLOGY 


The methodology for the fatigue design of primary structural com- 
ponents is schematically outlined in Figure 4. The three basic elements 
of safe-life design include; 

- Working Fatigue S-N Curve - developed based upon coupon tests to 
determine curve shapes, full-scale specimen tests to provide data points 
for constructing a mean curve, and then a statistical reduction of the 
mean curve to a "working" curve upon which the damage calculations are 
based. A minimum of six full-scale test specimens are generally 
required to establish the mean curve. The statistical reduction is made 
on the stress (and not cycles) for it is the high cycles region of the 
S-N curve at which fatigue loading occurs. 

- Design Loading Spectra - developed from the mission profile 
which is a prediction of how the aircraft shall be used, and the flight 
loads which are developed for the various maneuvers and flight condi- 
tions of the helicopter's mission profile. Obviously in the development 
of a new helicopter system, these loads must initially be either pre- 
dicted or empirically derived from measured loads on similar helicopters. 

- Cumulative Damage Theory - the combination of loads and fatigue 
strength by means of a cumulative damage theory results in the component 
fatigue life. 

The safe-life methodology outlined herein typifies that used by most 
US and European helicopter manufacturers (references 1 through 3) 
although specific approaches for analysis of the loads and the material 
fatigue data vary with the manufacturers. 

For the UH-60 Black Hawk and the YAH-64 Apache, the Army's more 
recent helicopter development programs, the minimum required fatigue 
life for all dynamic components has been 5,000 and 4,500 flight hours, 
respectively. Many of the components on these aircraft have exceeded 
their minimum fatigue life requirement. 



Figure 4 
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FAILURE RATE VS FLIGHT HOURS 

The safe -life design approach, when properly executed, places 
statistically predictable failures at an extremely remote probability 
of occurrence. This approach does, however, assume the structure to 
be flaw free. Over the years operational experience has shown (refer- 
ences 4 through 6) that the safe -life approach by itself is not adequate 
because of unpredictable failures resulting from such sources as manu- 
facturing defects and service -induced damage (see Figure 5) . Further- 
more, review of this operational experience has shown these unpredictable 
failures to far outnumber the predictable ones. 



HOURS 


Figure 5 
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DAMAGE SOURCES 


The statistically nonpredictable failures discussed in the previous 
figure result from one of three categories of damage - initial, as 
originating from material defects or manufacturing-induced anomalies, 
service- induced or ballistic. Examples for each of the damages are 
presented in Figure 6. Ballistic damage is certainly an in-service 
occurrence; however, because of the extreme nature of this damage and 
its associated short residual life requirement (30 minutes to roughly 
5 hours) , is generally addressed separately from the other damage types, 
and shall not be discussed herein. 


• DESIGN ERRORS 

• MATERIAL-AND MANUFACTURING-INDUCED DEFECTS 

- IMPROPER LAYUP 

- INCORRECT/INCOMPLETE CURING 

- MACHINING ERRORS 

- HANOLING/ASSEMBLING AREAS 

• SERVICE INDUCED DAMAGE 

- FOREIGN OBJECT DAMAGE |F00| 

- MAINTENANCE ERRORS 

- GROUND HANDLING 

• BALLISTIC DAMAGE 

- SMALL ARMS PROJECTILES 17.62mm fc 12.7mm| 

- HIGH EXPLOSIVE INCINDIARY |23mm HEI-T| 


Figure 6 
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EARLY DAMAGE TOLERANCE REQUIREMENTS 


Although damage tolerance as we know it today is a relatively new 
technology, its importance was recognized as long as four centuries ago. 
Near the end of the fifteenth century, technical notes were written on 
what must have been one of the first requirements for damage -tolerant 
design. They were in the notebook of Leonardo da Vinci in which he 
discussed the physics of flight and the design of "flying machines," and 
are shown in Figure 7 (reference 7) . 

In the early 60 's a fail-safe design philosophy began to evolve 
within the rotory wing industry and has been used in numerous applica- 
tions for both metallic and fibrous composite structures (references 8 
through 10). The safeguard in the fail-safe approach is that damage, 
induced or fatigue, will be detected by inspection procedures before it 
grows to such an extent that the residual strength of the structure is 
reduced below a safe level. 

In the early 70 's a third design philosophy - damage tolerance - 
began to evolve. It is similar' to the fail-safe approach but further 
considers the growth of damage resulting during manufacture or service 
usage. Therefore, fundamental to the damage-tolerant approach is an 
understanding of structural performance in the presence of defects or 
damage. 


"IN CONSTRUCTING WINGS ONE SHOULD MAKE ONE CORD TO BEAR THE 
STRAIN AND A LOOSER ONE IN THE SAME POSITION SO THAT IF ONE 
BREAKS UNDER STRAIN THE OTHER IS IN POSITION TO SERVE THE SAME 
FUNCTION".' 


FROM UONAIDO DA VINCI'S NOTtlOOK 
ON DfSIGN OF "FLYING MACHINES' 

15TH CEHTUIY 


Figure 7 
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DAMAGE TOLERANCE PHILOSOPHY 

Damage tolerance is a safety requirement and equates to the fail- 
safe design approach employed for flight critical components by the 
rotary wing industry since the early 60 's. Damage tolerance is the 
ability of the structure to resist failure due to the presence of 
defects, cracks, or other dcimage for a time period sufficient to enable 
their detection. 

The essential elements of dcimage tolerance are presented in 
Figure 8. In designing to this philosophy, the criteria must address 
both the static residual strength and damage propagation for the 
structure under consideration. 

The fail-safe and damage tolerance design philosophies each have 
the common objective of providing structural integrity at a reasonable 
level of assurance for all safety-of-f light structures, that is struc- 
ture whose failure could cause direct loss of the aircraft. Reference 
5 points out an important distinction between fail safety and damage 
tolerance. Fail safety as it was defined prior to the evolution of 
damage tolerance is based upon the premise that one cannot be certain 
that cracks (or damage) will not initiate at some time during the air- 
craft life and these cracks must be detected before the strength drops 
below a certain level. Damage tolerance, on the other hand, assumes 
the existence of initial flaws in the structure and the structure is 
designed to retain adequate residual strength until damage is detected 
and corrective actions taken. 

Although the damage-tolerant approach has begun to see a more wide- 
spread usage within the Army, there is still considerable work to 
be accomplished in developing specific requirements and procedures for 
this design approach. 


The basic philosophy of damage tolerant design is 
based upon: 

• THE ACCEPTANCE THAT DAMAGE WILL OCCUR 

• AN ADEQUATE SYSTEM OF INSPECTION SO THE 

DAMAGE MAY BE DETECTED 

• AN ADEQUATE STRENGTH MAINTAINED IN THE 

DAMAGED STRUCTURE 


Figure 8 
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RESIDUAL STRENGTH AND DAMAGE GROWTH 

The dcimage- tolerant requirement, which is illustrated schematically 
in Figure 9 , addresses both the residual strength and the damage 
propagation for the structure under consideration. 

The residual strength is the amount of static strength available at 
any time during the service exposure with damage present. Safety is 
ensured by designing to requirements wherein damage is never allowed to 
grow and reduce the residual static strength of the structure below a 
specified value, such as that corresponding to the maximum load to be 
experienced in service. 

The damage growth curve presents damage size as a function of time. 
The period of time required for damage to grow from a minimum detect- 
able size, such as by in-service inspection, to a critical size is the 
detection period. The critical size generally corresponds to the maxi- 
mum service load residual strength, although for composites this value 
could be dictated by stiffness loss to a critical value. 




KESIDUAL STtENGTH 


DAMAGE GROWTH 


Figure 9 



DAMAGE TOLERANCE APPROACHES 


Damage tolerance may be achieved in one of two different methods as 
shown in Figure 10. In fail-safe structures, damage tolerance (and 
their safety) is assured by the allowance of partial structural failure, 
the ability to detect this failure prior to total loss of the structure, 
the ability to operate safely with the partial failure prior to inspec- 
tion and the maintenance of specified static strength throughout the 
period. Fail-safe structures are usually comprised of multiple elements 
or load paths such that damage can be safely contained by failing a load 
path or by arresting a rapidly running crack. In safe crack growth 
structures, damage tolerance is achieved by sizing the structure, using 
fracture mechanics or empirical data, such that initial damage will grow 
at a stable, slow rate and not achieve a size large enough to fail the 
structure prior to detection. A special case within this category is 
nonpropagating- defect structures wherein the structure is designed to 
sufficiently low stress levels for virtually no propagation of the 
largest defect that would not be detected during a production inspection 
or that could be incurred in service. 



Figure 10 
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SAFE- LIFE AND DAMAGE TOLERANCE REQUIREMENTS FOR ARMY HELICOPTERC 

The structural design specification for helicopter MIL-S-8698 
imposes requirements for safe-life fatigue design of flight critical 
structures but not for fail-safe or damage tolerance design (Figure 11). 
In the more recent major helicopter system procurements, the Army has 
required fail safety or damage tolerance for the helicopter primary 
structure. These requirements have been stated in general terms as 
follows ; 

"The primary structure shall incorporate materials, stress levels, 
and structural configurations that will minimize the probability of loss 
of the aircraft due to damage of a single structural element (including 
control system or dynamic components) , or due to propagation of unde- 
tected flaws, cracks, or other damage. Slow crack growth, crack 
arrestment, alternate load paths and systems, and other available design 
principles shall be used to achieve this capability. " 

Because of the general nature of these requirements, they were 
successfully applied to both metallic and fibrous composite structures 
(reference 11). It has been left up to the helicopter manufacturer to 
define the specifics such as damage or flaw sizes, inspection intervals, 
etc., subject, of course, to Government approval. 


• MIL-S-8698, STRUCTURAL DESIGN REQUIREMENTS, HELICOPTERS 

- SAFE LIFE DESIGN REQUIREMENTS 

- NO FAIL SAFE/DAMAGE TOLERANCE DESIGN REQUIREMENTS 

• SYSTEM DEVELOPMENT SPECIFICATIONS (UH-60, YAH-64) 

- FAIL SAFE/DAMAGE TOLERANCE DESIGN REQTS 
■ REQUIREMENTS ARE GENERAL 

- TAILORED TO METALLIC STRUCTURES 

Figure 11 
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COMPOSITE MAIN ROTOR BLADES 


The Army's work in composite main rotor blades has spanned well 
over a decade. Composite blades are either in production or at various 
stages of planned implementation for nearly all inventory helicopters. 

Four composite blade designs are shown in Figures 12 through 15, 
and although they were each developed by a different helicopter manu- 
facturer, there is a certain degree of commonality in the fail-safe and 
damage-tolerant structural design approaches. Two blades are in pro- 
duction - the CH-47D blade (Figure 12) manufactured by Boeing Vertol and 
the K-747 blade (Figure 13) manufactured by Kaman Aerospace. The Hughes 
Helicopters AH-64 composite blade (Figure 14) and the Sikorsky UH-60 
composite blade (Figure 15) are currently under manufacturing develop- 
ment in MM&T programs and are planned for future production implementa- 
tion. 

All blades are of damage-tolerant, high fracture toughness mate- 
rials - either fiberglass or Kevlar. The UH-60 blade spar will be 
approximately 50% graphite; this is necessitated by the desire to main- 
tain the same 9^5% blade airfoil thickness as the titanium spar blade 
that it is to replace, and at the same time to match flap, torsion, and 
chord stiffnesses. To accomplish this with fiberglass would necessi- 
tate increasing the inboard airfoil thickness to 14.5%. 



Figure 12 
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K-747 COMPOSITE BLADE 



$PAI ATTACHMIHT DIIAIl fUU SFAN S CLASS 
STIAr WOUND AlOUND AlUMtNUIII flTTIHG 


CAST iroxv WEDGiS 



IIP WEIGHT • IIASS 
LEADING EDGE ElOSION GUAND* 
THEIMOPLASTtC POIYUIETHANI 
TIP CAP 



NOMEX HONEYCOMI 
jCOIE 

MUIT1>CEU mAMENT WOUND 
S GIASS/EPOXY SPAt 

fILAMENT WOUND XmAt/EPOXY 
TIAIUNG EDGE MCMIEt OIIGiN 


EllAMENT WOUND ^ AS” 
GLASS/ EPOXY SKINS 
(TOP A lOTTOM) 


Figure 13 


YAH-64 COMPOSITE MAIN ROTOR BLADE 


AHPOIl SICTIOW 



Figure 14 



COMPOSITE MAIN ROTOR BLADES (CONTINUED) 

All blades have a redundant load path root end except for the CH-47D 
blade; its single lug root end configuration was dictated by hub 
interface requirements. Hence fail safety in the root end has generally 
been achieved through redundancy. 

For the constant airfoil section, damage tolerance is achieved with 
the material selection and relatively low design operating strain 
levels, typically 1500 yin. /in. to 1800 yin. /in. 

Two of the blade concepts - the K-747 and AH-64 - also contain 
structural redundancy of the spar in the airfoil region through a 
multi-tubular (or cell) concept; this is to a large degree for ballistic 
tolerance. 

In general, the damage tolerance of these blades was developed and 
demonstrated through subcomponent and full-scale testing. Analysis was 
generally by classical techniques, with fail safety or damage tolerance 
assessed by assiiming loss of a certain portion of the structure. 


UH-60 COMPOSITE MAIN ROTOR BLADE 


AKFOIl SICIIOM 



SiMPtriEDSHEATH 


A 


ROOT END 


ROTOR DIAMETER = SS’I” 

RUDE CHORD = 21" 

RUDE WEI6HT = 211 LRS 
MANUFACTURINC RROCESS: WET FIIAMENT 
WINDIN6/C0CURIN6 


BUNDCD HBfcBGlASS lAMINAIE 



Figure 15 
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COMPOSITE ROTOR BLADE DAMAGE TOLERANCE DEMONSTRATED BY HLH ROOT END TEST 


Some of the early R&D emphasis on composite blades centered on the 
root end, for the success of composite applications to the rotor blade 
hinged to a large degree on the development of a root end that was 
structurally efficient, fail safe, and at the same time producible. 
Results from damage tolerance testing on the root end of the HLH blade 
which was fiberglass construction and weighed 780 lbs are summarized 
in Figure 16. After testing this specimen for 490 hours to demonstrate 
safe life, the specimen was damaged to varying degrees, from severing 
one of the four root end lugs to removing a 6" x 12" section, and again 
tested as shown in the figure to demonstrate its fail-safe properties. 
At the cutout section the axial, flap, and chord stiffness were reduced 
by approximately 50 percent, and torsion was reduced by 85 percent. 



• UNDAMAGED AFTER 

240 HOURS AT 100% OF DESIGN LOADS 
250 HOURS AT 300% OF DESIGN LOADS 


• SUCCESSFULLY WITHST.OOO 100% OF DESIGN LOADS 
40 HOURS WITH LUG CUT 

130 HOURS WITH LUG CUT AND FORWARD UPPER STRAPS SLOTTED (1/8 IN. SPANWISE BY 6 IN. CHORDWISE) AND THEN 
SUSTAINED LIMIT LOAD WITH NO DAMAGE PROPAGATION 
30 HOURS WITH LUG CUT AND FORWARD UPPER AND LOWER STRAPS CUT (12 IN. SPANWISE BY 6 IN. CHORDWISE) 


Figure 16 
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COMPOSITE ROTOR BLADE DAMAGE TOLERANCE 
DEMONSTRATED BY UTTAS ROOT END TEST 

Investigations have shown that initial deterioration of composites 
under repeated loading is manifested in a gradual loss of stiffness 
which is attributable initially to failure or breakdown within the 
resin. This factor, combined with another important feature that this 
loss in stiffness can occur, depending upon the stress level, many 
thousand cycles before static failure, yields the possiblity of inherent 
damage tolerance for many composite structures. 

This detection method of stiffness loss has been demonstrated in 
full-scale fatigue testing of fiberglass rotor blades. Reference 8 
reports that during fail- safe testing of the Boeing Vertol UTTAS rotor 
blade root end, which was fiberglass/epoxy construction, there was a 
noticeable decrease in torsional stiffness after approximately 2 x 10® 
cycles at elevated test loads. This testing is summarized in Figure 17. 
Although the blade still retained its full structural integrity, the 
change in stiffness was deemed a failure. The loads were dropped to high- 
speed level-flight values and the equivalent of 30 flight hours was run, 
with no further measured degradation. Test loads were then increased 
to the original high level and testing continued until blade deflec- 
tions exceeded machine capability. Even here there was no structural 
failure. Dynamic analyses, based upon reduced stiffness values, showed 
that the vibrational levels would have risen to a level noticeable to 
the pilot, but would still have been within a safe operating range. 
Additionally, regions of surface delamination became visible along with 
the stiffness loss. A key issue that must be addressed with this detec- 
tion procedure is tying the indication, such as increased vibration 
level, to the correct source. 



RATIO: 

DEFLECTION TO 
INITIAL DEFLECTION 


1 8g LOAD (277% DESIGN TORQUE) 


ONE MILLION CYCLES AT - 
100% DESIGN LOAD 


Figure 17 
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UH-60 COMPOSITE TAIL ROTOR SPAR 


The UH-60 composite tail rotor spar provides a good example of the 
need for damage tolerance criteria for composites . This component is 
in production on the UH-60. It is graphite/ epoxy and is highly complex 
in configuration, being made up of more than 300 individual detailed 
plies in order to satisfy strength and stiffness distribution require- 
ments. Because of this complexity and the manufacturing process used, 
the rejection rate of this component during early production was quite 
high. Because the effect of manufacturing defects (See Figure 18) had 
not been quantified for this component, the acceptance criteria was 
stringent. 

Three steps were taken to reduce this unacceptably high rejection 
rate: 

1. Alternate layup 

2. Press grade graphite 

3. Inspection and acceptance criteria 

Out of the three corrective steps, implementation of improved 
inspection techniques and an improved acceptance criterion, based upon 
actual testing of specimens with known defects (voids, porosity and 
inclusions) , had the most significant impact in significantly lowering 
the rejection rate. 


QUAL TEST SPEC #3 



Figure 18 
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UH-60 COMPOSITE TAIL ROTOR SPAR FATIGUE TEST 


Fatigue testing of a tail rotor spar is shown in Figure 20. 



Figure 20 
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UH-60 COMPOSITE TAIL ROTOR SPAR FATIGUE TEST RESULTS 

Figure 21 shows that spars with various defects fell within 
the scatter of spars that were void-free; this enabled establishment 
of a new acceptance criteria. 



Figure 21 


SUMMARY/CONCI.US IONS 


Helicopter flight- critical structures, although limited primarily 
to rotor blades, have been successfully designed for damage tolerance. 
This has been accomplished through redundancy, use of damage “tolerant 
materials and sizing to low operating stress levels. Design and 
demonstration have been based on laboratory testing. 

The increased application of composites to more complexly loaded 
structures, such as bearingless rotor hubs which inherently will operate 
at higher strain levels, further amplifies the need for a sound damage - 
tolerant criterion. Development of the criterion must address those 
items listed in the last bullet of Figure 22. 


0 LABORATORY TESTING THE BASIS FOR DESIGN AND DEVELOPHENT OF DAMAGE 
TOLERANT COMPOSITE STRUQURES 

0 DAMAGE TOLERANT CRITERIA FOR COMPOSITES NEEDED, ESPECIALLY IN 
LIGHT OF 

- COMPOSITES APPLICATION TO MORE COMPLEXLY LOADED COMPONENTS 

- INCREASED STRAIN LEVELS 

0 DAMAGE TOLERANT CRITERIA DEVELOPHENT MUST INCLUDE 

- PREDiaiVE CAPABILITY 

- DAMAGE SIZE REQUIREMENTS 

- NDT CAPABILITY/DAMAGE DETECTABILITY 

- DESIGN REQUIREMENTS 

- ENVIRONMENTAL EFFEQS 

- TEST/QUAL REQUIREMENTS 


Figure 22 
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UTS = 1102 MPa 
(160 ksi) 

... UTS/Vf = 2250 MPa 
(327 ksi) 


jjlj = 1.00 - .091 log N 


X 1 .00 - .096 log N 


20 40 60 80 

SLOPE OF S-N CURVE, B, (MPa/ decade) 


5 6 0 1 2 3 4 

LOG (CYCLES TO FAIL, N) 


Single cycle strength versus slope of 
S-N curve, nonwoven glass fiber 
composites, tension-tension fatigue 
at R = 0 to 0. 1 


Normalized tensile fatigue life curves 
for E-glass strand without matrix and 
with several matrices (5 Hz, R = 0.10). 
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FATIGUE SHAPE PARAMETER 














Characteristic lifetime (Cycles) 


Vieibull parameters for fatigue life distributions 







TRANSVERSE TENSION AND LONGITUDINAL SHEAR 


HIGHLY FATIGUE SENSITIVE 

MORE THAN 60% REDUCTION IN STRENGTH 
AT 10® CYCLES 

FAILURE RESULTING FROM FAST CRACK 
PROPAGATION 

NO INDICATION OF GRADUAL CHANGES IN 
MODULUS AND STRENGTH 



Transverse S-N data for B/Ep (Ultimate tensile 
strength = 60.9 MPa) 



CYCLES TO FAILURE 

Shear S-N data for B/Ep- (Longitudinal shear strength 
= 66.7 MPa, interlaminar shear strength = 81.4 MPa) 
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LAMINATE FATIGUE 


FAILURE PROCESSES 
FIRST PLY-FAILURE 

ASYMPTOTIC INCREASE IN CRACK DENSITY IN PLIES 

DELAMINATION 

LAMINATE FRACTURE 

















STRENGTH THEORIES OF COMPOSITES 


STATUS AND ISSUES 


Edward M. Wu 

Lawrence Livermore National Laboratory 
University of California 
Livermore, California 



Central Issues 


XLntL^ Uniaaaial iftMMt 


Fiber controlled? or Fiber dominated? ^ 
Matrix controlled? or Matrrix dominated? 



Should tensile & compressive strength be related? 



I 

^ h\ 

1 


X, 


XLmLuf- '€om/bimAd 


Are longitudinal & transverse strength 
independent? or coupled? a 



How to measure strength coupling? 


ROLE OF MATRIX BINDER IN LONGITUDINAL STRENGTH (REFS.1 AND 2) 


MATRIX BINDER PROVIDES LOCAL REDUNDANCY 



^6^ 


6= INEFFECTIVE LENGTH 
~ lOd 


FIBER BREAK WITHOUT MATRIX 

NO. OF LOAD CARRYING 
FIBERS 


WITH MATRIX 
NO. OF LOAD CARRYING 
FIBERS 


0 3 P/3 3 P/3 

1 2 P/2 3-8 ~ P/3- 


LONGITUDINAL STRENGTH DEPENDS ON MATRIX EFFECTIVENESS 
PARAMETER 6 ^E^^/Ej. 0^^^^ ®|NTERFACe| 




LONGITUDINAL COMPOSITEFAILS SEQUENTIALLY (REF. 3) 


n 





GLOBAL LOCAL 





LONG. COMPOSITE 
FAILURE IS SEQUENTIAL 



SMALL 6 MORE 
UNSTABLE 
BUT STRONGER 


LARGE a MORE 
UNSTABLE 
WEAKER 


MATRIX 

INCREASES 

a 


• LONG COMPOSITE FAILURE IS SEQUENTIAL NO WELL-DEFINED PLANE OF FAILURE 


LONGITUDINAL STRENGTH: FIBER-DOMINATED OR FIBER-CONTROLLED 

HIGH LONGITUDINAL STRENGTH OF FIBER IS DUE TO: 

DEFECTS ARE MINIMIZED BY SMALL FIBER DIAMETER 

• a 600 ksi (STEEL - 200 ksi) 

• • ult 


BUT EXISTENCE OF A SINGLE FLAW LEADS TO FAILURE 
(WEAKEST LINK OF CHAIN) 


/. LARGE SCATTER 


COEFF. OF VAR. SHAPE 
FIBER 20% 6-8 

STEEL 3-5% 25-50 


WITH WITHOUT 

MATRIX MATRIX 

BUNDLE STRENGTH 550 ksi 350 ksi 
BUNDLE SCAHER 4-5% 20-25% 


• LONGITUDINAL COMPOSITE STRENGTH IS FIBER-DOMINATED 
WITH SUBSTANTIAL MATRIX INFLUENCE 











Ln< STRENGTH > 


7*8 8 8.2 8.4 8.6 



STRENGTH 

KEWLRR 43 PILAMENT <5CM> 
INTRINSIC STRENGTH <.02CH'MIN,23C),N=54 


Ln< STRENGTH > 

7.8 8 8.2 8.4 8.6 



STRENGTH 

KELVAR 49-332/T403 STRAND 
INTRINSIC STRENGTH (lCn/NIN>iN«ie8 

17 «; 


CO 





TRANSVERSE STRENGTH: MATRIX-DOMINATED OR MATRIX-CONTROLLED (REF. 4) 


t t t t 









• 

• 


• 



• • 


• — ■ 






. ^ 


• 


• 






• 
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Oj Cj O3 

NUCLEATION GROWTH COALESCENCE FAILURE 

INFLUENCED BY FIBER PROPERTIES 

TRANSVERSE STRENGTH IS MATRIX-DOMINATED, 

INFLUENCED BY FIBER AND PACKING, NOT WELL-DEFINED PLANE OF FAILURE 


STRENGTH COUPLING UNDER COMBINED STRESS 



0 





• BROKEN FIBER INITIATES TRANSVERSE CRACK. 
THEREFORE REDUCES TRANSVERSE STRENGTH 



• STRENGTH COUPLING 
IS EXPECTED FROM 
PHYSICS OF FAILURE 


0 





TRANSVERSE CRACK REDUCES LOAD TRANSFER. 
INCREASES 6. REDUCES LONGITUDINAL STRENGTH 




FAILURE SURFACE FOR STRENGTH COUPLING (REFS. 5 AND 6) 


WEAKENED FROM 
REDUCTION OF 
I LOAD SHARING 





a 


1 


f(ai) = o,Fj.Fj.aaj+F„^a,.o.o,+ ...<1 
TENSOR POLYNOMIAL EXPANSION OF SURFACE 


NONCOUPLING STRENGTH CRITERIA 


“2 



MAX. STRESS 
(RANKIN) 



MAX. STRAIN 
(ST. VENANT) 


BOTH CAN BE EXPRESSED AS TENSOR POLYNOMINALS (REF.5) 



TRANSVERSE/TRANSVERSE COUPLING 
NO TRANSVERSE/LCNGITUDINAL COUPLING 
(REF.7) 






FAILURE TESTS REQUIRES WELL-DEFINED 
FAILURE MODES 




if a/ < 2ag* 

Tensile failure 
mode 


if a/ >2ag* 

Combine failure 
mode 



# ^ * 
<"6 

Tensile failure 
mode 


Plane of Failure does not necessarilj correspond to 
applied stress 

Plane of Failure of composites are seldom well defined 

Failure Criterion based on applied stress 

is NOT MechanisHc; it is Operational 


CONSISTENCY IN STRENGTH THEORY FORMULATION 


Output 


Mtl constant 

Input 


- 

®ijkl 

®ki 


Deformation | 

•1 

e.. 

>1 

- 

“ii 

AT 

Strength 

• 


( ) 

".1 


Ojj<aJ uniaxial 

o., < f (ctT ) combined 
•I H 

f {a..) - F,j o.j + F,.^, o.. 0 ^, + a., + . . . 


F . etc material constants dimensions 

1 


' 1 


stress 

» 

3 

Xt 


• Independent of mtls coordinates 

• Allow mathematical operations 






FAILURE STRENGTH MEASUREMENTS 


Design loading conditions (direction of stress vector) which is 
capable of discriminating failure criteria 


Stimuli 



One dimensional stress give rise to 3D strain 
One dimensional strain give rise to 3D stress 
Consistent characterization desirable 


THREE-DIMENSIONAL FAILURE CRITERION 

Tensor-polynominal Failure Criterion: 

FjOj + F.OjO. + . . . = 1 i.j = 1 ,2. 3,4,5, 6 



2-0 a failure surface 3-D a hyper-surface (6) 

• How many independent tests? (12) 

• What are the tests? 
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NO. OF INDEPENDENT TESTS 
Failure tensor F^, F^ follow tensor transformation rules 



• 3-D failure criterion for lamina (a) only 
Symmetry condition of orthotropic lamina 2 = 3 
Tensor Notation Contracted Notation 


22 

^33 

f=2 

= F 3 


1122 

~ ^1133 

F,2 

= «=13 

> 12-4 = 8 

2222 

~ ^3333 

*^22 

= >=33 

1 

1313 

~ ^1212 

*=55 

“ ^66 



Component not associate with planer properties 

^2233 ^^23^ Transverse strength coupling 
^2323 <^ 4 * 


INDEPENDENT EXPERIMENTS 
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EXPERIMENTS ASSOCIATE WITH THICKNESS DIRECTION 




Experiment No. 7 Experiment No. 8 

Biaxiel compression Shear test 



NOW. INTERPRET AS INTERNAL STRESS (A MATERIAL CONSTANT). 
STRESS ANALYSIS MUST OPERATE ON (o^ - 
WHICH IS AN OPERATIONAL INCONVENIENCE 
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Off-axis sample to measure shear strength 
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DETERMINATION OF COUPLING COEFFICIENTS 


• EXPERIMENTALMEASURMENTATOPTIMIZED STRESS RATIO ®1 (REF.5) 


SIMPLIFIED STRESS RATIO; 


FOR DUCTILE MATRIX 


MEASURE AT — ~ 

®2 ^2 



FOR BRIHLE MATRIX 



HEURISTIC ESTIMATION - BOUND ^ F^^ 


IF Xj HAS SMALL SCAHER THEN F^^ ^ SMALL 

IF X^ HAS LARGE SCAHER THEN ►LARGE 


WITHIN BOUND 


COMBINED STRESS, ANISOTROPIC STRENGTH 





SUMMARY 


MODELING OF SEQUENTIAL FAI LURE - RECENT PROGRESS 

SHAPE OF COMBINED STRESS FAILURE SURFACE - LONGITUDINAL 
AND TRANSVERSE COUPLING EXPECTED AND OBSERVED 

FOR STRENGTH-COUPLED FAILURE SURFACE, TENSOR POLYNOMIAL 
IS OPERATIONALLY AHRACTIVE: 

COUPLING COEFFICIENTS (F,.)CAN BE MEASURED OR HEURISTICALLY 
ESTIMATED 

READILY EXTENDABLE TO 3-D AND HIGHER ORDER 
MOST ISSUES IN PROPER ORDER 

PROBABLISTIC REPRESENTATION OF FAILURE SURFACE NOT YET RESOLVED 
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ANALYSIS OF DELAMINATION 


Frank W. Crossman 

Lockheed Palo Alto Research Laboratory 
Lockheed Missiles and Space Company, Inc. 
Palo Alto, California 



INTRODUCTION 


This paper reviews our understanding of the influence of delamination on the 
fracture sequence and final fracture of graphite-epoxy laminates. The emphasis here 
is on tensile fracture. 

The paper is organized as follows: 

(1) Review of the common fracture modes in laminates and their interaction 

(based on Lockheed and NASA studies) 

(2) Review of free-edge stress analysis and prediction of delamination sites 

(based on studies by Lockheed Palo Alto Research Laboratory and 
Drexel University) 

(3) Consideration of the problem of predicting delamination fracture in 

regions containing a high stress gradient (based on studies by Lockheed 
for NASA) 

(4) Analysis of delamination fracture by calculation of strain energy release 

rate (based on studies by Drexel University and Lockheed) 

(5) Consideration of the Interaction of interlaminar and intralaminar fracture 

(based on studies by Drexel University and Lockheed for the Air Force 
Office of Scientific Research) 

(6) Analysis of the interrelationship between strength and stiffness (based on 

studies by Lockheed for NASA) 

(7) Consideration of the correlation between strain energy analysis and 

statistical analysis approaches in predicting the criticality of a 
laminate crack (based on studies by Lockheed, Cornell University, and 
Drexel University) 

The following people contributed to the observations, analysis methods, and 
methodology of delamination fracture analysis described in this paper: J. G. 

Bjeletlch and W. J. Warren (Lockheed Palo Alto), A. S. D. Wang and G. E. Law 
(Drexel University), T. K. O'Brien (NASA Langley), and J. T. Ryder (Lockheed Rye 
Canyon) . Finally, many researchers have contributed to the understanding of 
delamination and its influence on laminate fracture. These people are acknowledged 
in the partial bibliography at the end of this paper. 
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Five cases that involve delamination fracture are presented. These are: 


(1) Free-edge delamination 

(2) Compression delamination 

(3) Delamination interaction with intralaminar cracking 

(a) Ply peeling 

(b) Delamination after transverse cracking 

(c) Simultaneous delamination/ transverse cracking 
These cases are illustrated in the five figures which follow. 


INTERLAMINAR FRACTURE - CASE 1 : FREE EDGE PLY DELAMINATION 




SINGULAR INTERLAMINAR STRESSES AT CRACK-TIP 
OPENING MODE SHEARING MODE 

ACTUAL CRACK IS IN A COMBINED MODE 
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INTERLAMINAR FRACTURE - CASE 2: COMPRESSIVE LAMINATE FAILURE MECHANISMS 



A 



CRACK TIP 


SINGULAR INTERLAMINAR CRACK TIP 
NORMAL STRESS 

CRACK GROWTH IS PRECIPITATED BY 
STRUCTURAL BUCKLING ACTION 

UNSTABLE CRACK GROWTH DETERMINES 
THE COMPRESIVE STRENGTH OF LAMINATE 


INTERLAMINAR FRACTURE - CASE 3 (a): PLY-PEELING (ALSO A CASE 

OF PLY DELAMINATION) 
(SURFACE DELAMINATION) 

A«x 



V 
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INTERLAMINAR FRACTURE - CASE 3(b): 

*FREE EDGE DELAMINATION (ALONG 25/90 PLANE. MIXED MODE) 



INTERLAMINAR FRACTURE - CASE 3(c): 

"DELAMINATION EMANATING FROM THE ROOT OF A TRANSVERSE CRACK (MIXED MODE) 
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ULTIMATE TENSILE STRENGTH VERSUS STACKING SEQUENCE 


Studies were conducted In the mid-1970's to examine the Influence of laminate 
stacking sequence on the tensile properties. The results shown below and in the next 
figure compare the ultimate strengths of a number of quasi-lsotropic laminates 
constructed of T300/934 tape and fabric (ref. 1). A large variation in strength was 
found. This variation was not predicted from laminated plate analysis where stacking 
sequence does not alter the in-plane ply stresses when the specimen is loaded under 
tension. 


TENSILE STRENGTH OF QUASI-ISOTROPIC T300/ 934 LAMINATES AS A FUNCTION 
OF PLY STACKING SEQUENCE (LAMINATES DESIGNATED HMF WERE FABRICATED 
FROM WOVEN GRAPHITE FABRIC) 

{From ref. 1} 


MATERIAL 

THICKNESS 

(IN.) 

PANEL 

NO. 

LAYUP SEQUENCE 

+ STANDARD DEVIATION 

60 70 80 90 

1 1 ' ■ ' 1 ( 

T300/934 

0.042 

13 

[0/+45/-45/90) 

8 

1 ' I ' ^ ' 1 

T300/934 

0.042 

15 

■ 1 ■ 1 ■ 1 ■ 1 

(90/45/-45/0J 

^ • 1 . 1 . 1 • 1 

T300/934 

0.042 

14 

' 1 ■ 1 ■ 1 ' 1 

(0/90/45/-45) 

® 1 1 • 1 . 1 1 J 

T300/934 

0.042 

17 

' 1 • I ■ 1 ' 1 

[90/0/45/-45] ^ 

^ 1 1 1 1 1 1 1 1 

HMF330C/34 

0.054 

2 

' 1 I 1 ' i i 1 

[(0/90)/ (+45)] 

“ 1 1 1 1 t- ' ■ ' 

HHF330C/34 

0.108 

4 

1 1 1 1 1— 1 1 1 

((0/90)/(90/0)/(+45)/(+45)] O 

1 1- - 1 1 1 1 1 1 

T300/934 

0.042 

18 

' r ' 1 ' 1 ' 1 

[45/-45/90/0J ^04 

^ 1 1 1 i 1 i 1 1 

T300/934 

0.042 

16 

1 1 1 1 I 1 1 1 

[45/-45/0/90I 

1 1 1 j 1 i 1 i 

HMF330C/34 

0.054 

1 

' 1 ' 1 ' 1 ' 1 

[(+45)/ (0/90))^ Q 

. i • i 1 i 1 i 

HMF330C/34 

0.108 

3 

• 1 • 1 1 1 • 1 

[(+45)/(+45)/(0/90)/(90/0)l 

Li 1 1 1 ! ! ! ! 
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DELAMINATION AND TENSILE STRENGTH 


An examination of several independent studies confirms the observation that 
stacking sequences that are prone to delamination exhibit generally lower ultimate 
strengths. 


TENSILE STRENGTH VARIATION WITH STACKING SEQUENCE 

EIGHT-PLY, GRAPHITE-EPOXY. QUASI-ISOTROPIC LAMINATES 



90 h 


TENSILE 80- 
STRENGTH 

°ult- 

ksi 70 - 


60 


50 


O LOCKHEED 
• NASA 
V BOEING 
O NASA 
O NASA/VP I 
O AFFDL 
A GWU 




[+45/0/-45/90]s [+45/0/90]s [0/+45/90]s [+45/90/-45/0]s [0/90/t45]s [90/t45/0]s 

V V ' 


DELAMINATION 


NO DELAMINATION 
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WHAT CAUSES DELAMINATION 

Strain gages placed on tensile specimen edges to measure the "through-the- 
thlckness" normal strain (z-strain) as a function of tensile strain (x-straln) 
showed that the delamlnat Ion-prone laminates developed a negative Poisson ratio 
(e.g., the z-strain was positive during tensile loading). Thus interlaminar 
interaction that was unaccounted for in laminated plate analysis was providing a 
driving force for delamination. 


THROUGH-THICKNESS NORMAL STRAIN VERSUS APPLIED TENSILE STRAIN 
£x AS A FUNCTION OF T300/934 PANEL TYPE DESCRIBED IN PREVIOUS FIGURE 

iFrom ref. 1 } 
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WHERE DOES DELAMINATION OCCUR 


Examination of the edges of tensile specimens before and after failure revealed 
that the location of delamination depended on the particular stacking sequence and 
in some cases involved multiple delamination at several locations through the ply 
thickness. Transverse cracks were seen to initiate both before and after delami- 
nation, depending on stacking sequence. The extent of delamination was found to be 
greater at the edges than in the interior of the specimen. 


x-z EDGE VIEW OF (0/45/-45/90)^ LAMINATE AFTER LOADING TO 9S% UTS 

(From ref. 1 ) 
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EDGE VIEW OF (45/-45/0/90) LAMINATE AFTER LOADING TO 95% UTS 

(From ref. 1 ) 



LONGITUDINAL x-z SECTION OF (45/-45/0/90)^ TAKEN NEAR MIDPLANE 
BY SECTIONING AND POLISHING 

(From ref. 1 ) 
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EDGE VIEW OF (90/45/-45/0)^ LAMINATE AFTER LOADING TO 95% UTS 

(From ref. 1 ) 
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ANALYSIS OF INTERLAMINAR STRESSES 

Consideration of stress equilibrium and conservation of angular momentum led to 
an understanding of how interlaminar stresses arise and the development of simple 
models to approximate the level and sign of these stresses. 

ANISOTROPY IN STIFFNESS OR EXPANSION COEFFICIENTS 
CAUSES PLY OR LAYER STRESSES 


aL. 


V| E| 4L| + 




E 


0-| ■ -j- (4L^ - 4L|) 

0 


■ r 

0 






1 

r 

2 i 

H H- 


i 

? 



1 

•^1 

2 


K- 


IN-PLANE PLY STRESSES RESULT IN FREE-EDGE STRESSES 
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NUMERICAL ANALYSIS OF INTERLAMINAR STRESSES 



Finite-element and finite-difference analyses have been applied to models of an 
y— z section of a laminate under uniaxial loading in the x direction. The through- 
the-thickness Poisson ratio predicted from the analyses agreed well with experi- 
mentally measured values. 



5.0H 


CO 

M 

N 



2.5 5.0 7.5 10.0 12.5 15.0 


Y AXIS 

STACKING SEQUENCE DEPENDENCE OF THROUGH THICKNESS 
CONTRACTION IN T300/934 QUASI -ISOTROPIC 


LAMINATES (ref. 1) 




Panel No. 

Stacking 

Sequence 

V 

xz 

Experimental Calculated 

Ultimate 

Tensile 

Strength 

(ksi) 

IS 

l90/45/-45/0|g 

0.80 

0.72 

87 

13 

(0/45/-45/90) 

-0.05 

o 

o 

K) 

72 

16 

145/-45/0/90) 

s 

-0.30 

-0.34 

63 
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COMPARISON OF NUMERICAL METHODS 

The detailed stress analysis was found to depend somewhat on the fineness of 
the model discretization and the order of the assumed variation of displacement or 
stress field in the case of finite-element analysis. Stress gradients In both the 
y and z directions were found in these analyses, as shown In the figures that follow 


THIS STUDY 

PIPES AND PAGANO (ref. 2) 


FREE EDGE STRESSES IN [o, 90. 90. o]^ T 300/934 







PREDICTION OF DELAMINATION LOCATION 

As shown In the figures, the location of peak tensile stress In the analysis 
often agreed with the experimentally observed delamlnatlon location. 


DISTRIBUTION OF Oz NEAR FREE EDGE UNDER TENSILE LOAD 

IN X DIRECTION 


EIQ 





THE PROBLEM OF STRESS GRADIENTS 


A problem arose, however, when the numerical analysis of Interlaminar stresses 
associated with moisture gradients was considered. The figures show the moisture 
concentration as a function of y,z position at two different times during the 
moisture absorption process. At a time that was small with respect to the equilib- 
rium time, the gradient of concentration is very steep near the free edge. In the 
limit of very small times, this gradient is infinitely steep. Since absorbed 
moisture is known to swell the composite locally, a gradient of hygroscopically 
induced stresses is introduced. 

AVERAGE WEIGHT GAIN (LOSS) DURING ABSORPTION (DESORPTION) 

RELATIVE TO EQUILIBRIUM MOISTURE CONTENT M^ 

{From ref. 3} 



MOISTURE CONCENTRATION C/Cg AFTER ABSORPTION TIME 
(M/Ms«4%) ^ 


{From ref. 3f 



12 


y/h 


16 
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PREDICTION OF DELAMINATION UNDER A STRESS GRADIENT 


The results of stress analysis In a unidirectional laminate during desorption of 
moisture from an initially saturated condition show that a transverse normal stress 
greater than the transverse tensile strength of the material is present during the 
first part of the desorption period, but gradually reduces to zero upon drying out. 
Although strength theories based on these local stresses predict delamination during 
desorption, no experimental evidence of delamination was found in a series of 
experiments conducted under a NASA Ames contract (ref. 3). It was recognized during 
the analysis of these results that stresses calculated at a distance from the free 
edge or interlaminar region which is on the order of the fiber diameter cannot be 
physically real. The analysis is based on the assumption that fiber and matrix 
phases are smeared together at the lamina level. That assumption is not appropriate 
on the microstructural scale. At this level, the interaction of individual fibers 
and intervening matrix must be considered. Thus the analysis of interlaminar 
stresses is inherently limited to a cross-sectional area which contains on the order 
of 10 fibers. In graphite-epoxy laminates this results in a minimum finite-element 
size of approximately 0.001 x 0.001 in. 

a, THROUGH-THICKNESS DISTRIBUTION NEAR 
FREE EDGE OF UNIDIRECTIONAL 8-PLY 
LAMINATE DURING DESORPTION FROM 
EQUILIBRIUM MOISTURE CONTENT 

{From ref. 

(MPa/M%) 


0 30 60 
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PLY THICKNESS AFFECTS DELAMINATION ONSET STRESS 


Rodlnl and Elsenmann (ref. 4) found that the onset stress In a series of 
quasl-lsotroplc laminates varied with ply thickness as approximately 1 over the 
square root of laminate thickness. Delamlnatlon stress analysis based on continuum 
theory Is Independent of the scale of the structure and was unable to predict the 
experimental results. Rodlnl et al. proposed that Weibull strength theory which Is 
dependent on specimen volume could be used to explain the results. 

A second approach, based on fracture mechanics, was taken at Drexel and Lockheed 
to account for the observed size dependence of delamlnatlon and transverse cracking 
onset in a series of tensile-loaded (25/-25/90n)s laminates where n ranged from 
1/2 to 8 (ref. 5). 


PLY THICKNESS EFFECTS 


T 300/5208 LAMINATES 

AT ONSET OF DELAMINATION 


EXPERIMENT 

PREDICTION 

(±45/0/90); 

45.3 ksl 



45.3 

49.0 ksi 


45.5 


(± 452 / 02 / 902 ) 

38.3 


s 

38.6 

36.0 


37.8 


(± 453 / 03 / 903 ) 

28.7 


s 

30.7 

29.0 
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CASE STUDY USING T300/934 GRAPHITE-EPOXY LAMINATES 


Laminate configuration: 

(±25/90n)s n = 4, 1. 2, 3. 4. 6. 8 

Loading condition: 

Uniaxial laminate tension .5^^ or 

Uniform curing temperature AT = T^ - T^ ^ 

o ■ » ^ 

Types of cracking studied: 

1. Transverse cracks in 90° layer 

2. Edge delamination along 90/90 interface (mode I) and 25/90 

interface (mixed mode), and transverse-crack-induced delamination 

3. Edge-delamination-induced transverse cracks 
Objectives of study: 

1. Onset and growth behavior of various modes of damage 

2. Their interrelations 

3. Physical and/or geometrical parameters influencing these behaviors 

4. Whether minimechanics approach is viable tool in studying damages in 

composite laminates 

5. How to bridge gap between micromechanics and macromechanics 









OBSERVED TENSILE STRAIN Ex FOR TRANSVERSE CRACKING, DELAMINATION, 
AND ULTIMATE FAILURE VERSUS NUMBER OF 90° PLIES IN T300/934 



12 3 4 


n 



ANALYSIS OF DELAMINATION STRAIN ENERGY RELEASE BY CRACK CLOSURE 


The crack closure method pioneered for delamination analysis by Rybicki et al. 
(ref. 6) was applied to the analysis of strain energy release rate associated with 
both delamination and transverse cracking. The analysis accounted for the con- 
tribution of thermal stresses to the total energy release rate. The crack closure 
method also allows the determination of the G values associated with each crack 
opening mode. 


Y.v 



Y,v 



■ a + Aa ' ■■ 

The work required to close the crack extension Is approximated by 

AW«»[F (u,-u ) + F (v -V ) + F (w -w )]/2 
xfg yfg zfg 

where F^,F^,F^ are the components of the nodal forces required to close nodes 
f and g together. Thus, the energy release rates for the three crack extension 
modes are approximated by: 

G_ F (v,-v )/2Aa 

I yfg 

hi - P^(Uf-“g>^2Aa 

F (w,-w )/2Aa 
III zfg 
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ENERGY ANALYSIS 


G = - ^ 

da Aa 

Energy released in crack opening equals work required to close crack: 
-Au = %F • D 

Let Fg, Dg be the elastic solution for a loading of = 10"® and let F 
be the elastic solution for a loading of AT = -1°F. 

Then 

® ' ilr + D^(iT)) 

This may be rewritten as 



THE CRITICAL FLAW CONCEPT 

The next figure below Introduces the concept of an effective flaw which 
preexists in a material prior to loading or one which grows subcritically in the 
sense of satisfying the Griffith critical stress criterion based on independently 
measured stiffness and toughness. 

If a critical flaw is not assumed to exist, then the analysis of transverse 
cracking of progressively thicker 90° layers in the (25/-25/90n)s laminate predicts 
that the energy release rate goes to infinity as the ply thickness is increased 
without limit. The consequence of this is that the predicted stress level for 
transverse cracking approaches zero for very thick 90° layers. This is inconsistent 
with the finite strength of transversely loaded unidirectional composites. 


METHOD OF LINEAR FRACTURE MECHANICS 

• CONCEPT OF ENERGY RELEASE RATE 

• G-CURVES AS FUNCTIONS OF CRACK EXTENSION 



FOR MODE-I COLLINEAR CRACKING IN ISOTROPIC ELASTIC SOLIDS. 
Gj = (tt/E) ao^ 

CRITERION FOR CRACK GROWTH: 



HENCE 

»c= y(<5cE)/K) 
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DETERMINATION OF THE CRITICAL FLAW SIZE 


To determine the critical flaw size a^, we postulate that It can be calculated 
by knowledge of the Independently measured elastic modulus E the fracture toughness 
6];^ and the unnotched tensile strength In the longitudinal and transverse direc- 
tions. For T300/934 composites, the substitution of these quantities Into the 
Griffith equation 


a 


c 



gives Sj, (transverse) = 0.010 In. and a^, (longitudinal) = 0.020 In. Since these 
are both much larger than the fiber diameter, our assumption of homogeneous ply 
properties In deriving the effective flaw size Is consistent with the absolute size 
of the flaw. The transverse critical flaw size was applied to the analysis of 
transverse and delamlnatlon crack Initiation in the analysis. The relative sizes 
of ac and am are important In assessing the stability of delamlnatlon growth, 
which is found to depend on the 90° ply thickness in the (25/-25/90n)s experiments. 

BEHAVIOR OF G-CURVES IN FREE-EDGE DELAMINATION 





CONSIDER l±em. 


EDGE DELAMINATION MAY OCCUR 


1. ALONG 90/90 INTERFACE 
(MODE I, OPENING) 

2. ALONG -0/90 INTERFACE 
(MODES I, m MIXED 
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FRACTURE MODE AND SEQUENCE PREDICTION 


Energy release rate analysis is qualitatively and quantitatively successful 
in predicting experimental behavior. 

OBSERVED TENSIII STRAIN Ex FOR TRANSVERSE CRACKING, DELAMINATION, 

AND ULTIMATE FAILURE VERSUS NUMBER OF 90O PLIES IN T30O/934 


z 

o 



,2L 

0 


_L 

.5 


J_ 

4 

n 


J 1 1 I 

5 6 7 8 


COMPARISON BETWEEN THEORY (DOTTED LINES) AND EXPERIMENT 
(PREDICTION BASED ON - 1.4; = 0.008" AND AT - 225°F) 



n 
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IMPROVING THE RESISTANCE TO DELAMINATION 


From the fracture mechanics analysis of delamination onset, it is clear that the 
onset stress for delamination can be raised if the critical energy release rate (or 
delamination toughness) of the material can be improved. The table given below 
shows (1) that order-of-magnltude Improvements in matrix toughness do not translate 
into large Improvements in the composite prepregged from tape, and (2) that 
composites fabricated from cloth prepreg show greater improvement in toughness 
than do the tape materials. The references presented emphasize the Importance of 
microstructural factors in controlling (and limiting) improvements in delamination 
toughness. 


FRACTURE ENERGIES OF GRAPH ITE -EPOXY MATRIX COMPOSITES 


Material 

Ref. 

Resin or 
Interlaminar 
Fracture 
Energy 
(KJ/m^) 

Comments 

5208 

7 

0.076 

TGMA/DDS Matrix 

Code 69/HTS Graphite Tape 

8 

0.076 

M 

5206/Morg II Graphite Tape 

9 

0.17 

II 

F263/T200 Graphite Cloth 

7 

0.36 

II 

5206/Morg 1 1 Tape 




UD Cross Fiber 

9 

20.4 


(+ 45)5 Cross Fiber 

9 

7.8 


Ql Cross Fiber 

9 

9.6 


A 

10 

0.16 

MY750 / NMA / BDMA 

A/HMS Tape 

10 

0.24 


B 

10 

0.33 

MY750 / Piperidine 

B/HMS Tape 

10 

0.28 


C 

10 

1.40 

MY750 / 3.2% CTBN 

C/HMS Tape 

10 

0.37 


D 

10 

2.2 

MY750 / 6.2% CTBN 

D/HMS Tape 

10 

0.36 


E 

10 

3.20 

MY750 / 9% CTBN 

E/HMS Tape 

10 

0.49 


205 

7 

0.27 

Diglycldyl ether bisphenol A/ 

205/T300 Graphite Cloth 

7 

0.60 

Novelac/Dicyandiamlne Catalyst 

205+CTBN 

7 

5.1 

13x12 3000 tow BOO cloth 

205+CTBN/T300 Graphite Cloth 

7 

4.6 

CTBN + HYCAR additives 
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STIFFNESS/STRENGTH RELATIONSHIPS 


A current Lockheed/NASA study is emphasizing the development of a structural 
analysis approach to assess the Influence of transverse cracking and delamlnatlon 
the resulting laminate stiffness and strength properties. 


OBJECTIVE 

• DEVELOP " STRUCTURAL” ANALYSIS FOR LAMINATES WITH 
ARBITRARY LAYUPS 

• PREDICT RELATIONSHIP OF MATRIX CRACKING AND DELAMINATION 
TO LAMINATE STIFFNESS AND TENSILE STRENGTH 

PREMISE 

• A LAMINATE IS A STRUCTURE 

• MATRIX CRACKS AND DELAMINATIONS MAY OR MAY NOT BE 
DAMAGE 

• A "STRUCTURAL” ANALYSIS OF ARBITRARY LAYUPS IS NEEDED 


SHEAR STIFFNESS REDUCTION 


Finite-element analysis was used to determine the effective in-plane transverse 
and shear moduli as a function of transverse crack density. The results indicate 
that for the experimentally determined crack densities, the loss of in-plane stiff- 
ness due to transverse cracking is very small. 


SHEAR STIFFNESS REDUCTION 
AS A FUNCTION OF CRACK SPACING 



S 
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LAMINA STIFFNESS LOSS AS A FUNCTION OF CRACK SPACING 

MASTER CURVES 



STIFFNESS LOSS DUE TO MATRIX PLY CRACKING 

T300 - 5208 GRAPHITE EPOXY 


LAYUP 

90° PLIES 
S^ ^22 

f22)o 

+ 45 O PLIES 

1 ®12 

' («12)o 

-45^' 

S 

h 

PLIES 

«12 

(<^ 12)0 

^ ^xx, % 

[0/90]j 

3.33 

65 

— 

— 

— 

— 

1.2 

[»± 45 ] j 

— 

— 

14.6 

95 

14.2 

95 

0.9 

[0/90/±45]^ 

2.37 

52 

13.6 

95 

— 

— 

2.2 
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STIFFNESS CHANGE WITH DELAMINATION GROWTH 

Correlation of delamlnatlon size and stiffness shows a much more significant 
effect on stiffness associated with delamlnatlon, which can be rather simply 
calculated from a modified laminate analysis. 

MEASUREMENT OF STIFFNESS CHANGE 
WITH DELAMINATION GROWTH 





RULE OF MIXTURES ANALYSIS OF STIFFNESS LOSS 
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STIFFNESS LOSS AND DELAMINATION SIZE 


A comparison between model prediction of stiffness loss and experimental data 
is given below. 


STIFFNESS LOSS AS A FUNCTION OF DELAMINATION SIZE 

[+45/0/90] T300-5208 GRAPHITE EPOXY 
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TRANSVERSE-CRACK EXTENSION 


Transverse cracking of a ply typically terminates at the interface of a ply of 
differing orientation. Microscopic and X-ray examination of the crack tip by 
researchers at VPI&SU and Lockheed, among others, has failed to reveal significant 
damage to the fibers in the neighboring ply, at least when the specimen has-been 
subjected to only static loading (private communication, K. L. Reifsnider). Recent 
analysis of the stress concentration and energy release rate associated with the tip 
of the transverse crack has confirmed the lack of damage potential associated with a 
transverse ply crack. 

The figure below shows the finite-element model discretization of a transverse- 
ply crack. The ply arrangement for the quasi-isotroplc configuration is shown, but 
three other layups were examined with the same element model. The stress concen- 
tration associated with a transverse 90°-ply crack is much less in the quasl- 
isotropic case than in the transverse (90) g laminate. This is a result of the 
relatively small load carried by the cracked ply in the quasi-isotropic laminate. 

The minimum element size used in this analysis contains approximately 10 fibers. 
Analysis of stresses with elements on a fiber scale would lose physical meaning, 
given the homogeneous-ply-stiffness assumptions of the analysis. 


LOCAL STRESS CONCENTRATION DUE TO MATRIX CRACK 

[0/90/ + 45 ]^ 




LAYUP 

SCF 

0/90/±45‘ 

s 

1.1 

0 / 0 / ± 45 ' 

s 

1.6 

[Oj, 

8 


1.7 

[ 90 ] 

8 


1.5 
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WEIBULL STRENGTH ANALYSIS 


Phoenix (ref. 11) has derived a model which can be used to obtain the mean 
strength of a unidirectional composite as a function of its voliune. The volume is 
given by the product of the number of fibers in cross section (n) and the length 
divided by the critical length (m) . For a typical unidirectional tensile specimen, 
nm is on the order of 10^. The figure below shows the characteristic strength 
distribution function of a single critical length of fiber and the shift of that 
function associated with a volume of nm = 10^. Taking W(x) = 0.5 as the mean 
strength of a given volume, one can construct a table of strength normalized to a 
single fiber strength or that of a typical tensile specimen. 


CHARACTERISTIC DISTRIBUTION FUNCTION Wlx) FOR BUNDLE STRENGTH 
AND ASSOCIATED WEIBULL APPROXIMATIONS (ref. 11). (ALSO SHOWN IS 
STRENGTH DISTRIBUTION FOR A BUNDLE WITH # FIBER ELEMENTS) 



0.2 0.3 0.4 0.5 0.6 0.7 


DIMENSIONLESS LOAD. x/Xg 
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TENSILE STRENGTH VERSUS VOLUME 


Given that the stress concentration factor for a transverse crack in a quasi- 
isotropic laminate is approximately 1,1 and the voliime associated with the finite 
element nearest the transverse crack tip is approximately 5 x 10^ in nm units, we 
can see from the volume dependence of strength derived in the previous table that 
"failure" of that element is not expected on the basis of the Phoenix theory 
(ref. 11). 


UNDIRECTIONAL TENSILE STRENGTH VS SPECIMEN VOLUME** 


CASE 


Ln (nxn) X/X(CASE A) X/X(CASE B) 


A, A single fiber of the 

critical length 0,0 

B, A tensile specimen of 

,080 X ,5 X 10, in, 17,9 

C, A single ply ,005 in, 
thick of length ,04 in, 
near tip of transverse 

crack 9,65 

D, A single layer of fibers 
of length ,04 in, at the 
tip of the transverse 

crack 7,13 


1,0 


,25 


,35 


.41 


4.0 


1.0 


1.4 


1.6 


** A fiber spacing of .0004 in. and a critical length of .04 in. 
were assumed. 
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STRAIN ENERGY RELEASE RATE ANALYSIS 



The figure presented below shows the results of strain energy release rate 
analysis at the tip of a transverse crack for the geometry shown in the inset. 
Several types of laminates were analyzed. The energy release rate is plotted as a 
function of tensile strain. The longitudinal fracture toughness of graphite-epoxy 
is on the order of 80 to 90 ft-lb/ln^. A transverse crack in a (02/902)s laminate 
does not release enough energy to allow the transverse crack to propagate across the 
0° ply. Thus the energy release rate analysis approach provides the same prediction 
of transverse crack stability as the Welbull strength analysis. 


STRAIN ENERGY RELEASE RATE ANALYSIS OF MATRIX CRACKS 
IN EIGHT PLY GRAPHITE EPOXY LAMINATES 
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FAILURE BY STRUCTURAL LOAD TRANSFER 

The analysis of the final fracture event for tens lie- loaded laminates does not 
appear to depend on the details of ply microcracking. As experimental observations^^ 
of laminate damage growth under tension have accumulated, the Importance of local 
delamlnatlons in transferring a local load increase to the major load-carrying plies 
had been postulated by Reif snider and colleagues (ref. 12). The schematic edge 
section of laminate damage under static (top) and fatigue (bottom) loading shows 
regions of the fatigued specimen where essentially all of the applied load must be 
routed into the 0° plies around regions of extensively delaminated and tranverse- 
cracked 90° and 45° plies. Prediction of these loads will require more detailed 
structural modeling. 


PREFRACTURE PAHERNS NEAR FRACTURE SITE (ref. 13) 
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"STRUCTURAL" FAILURE ANALYSIS 


The degree to which the details of crack interaction must be modeled to predict 
tensile failure are being investigated in several ongoing research programs. 
Structural analysis of composite response can be carried out routinely with numerical 
methods lite finite-element analysis at several different structural levels; 
mlcrostructural, laminate structural, and macrostructural. However, this does not 
provide the total solution ta fracture analysis of composites. Our understanding of 
the structural implications of cracking at each level must be Improved to provide 
the proper failure criteria for the structural analysis. 

FRACTURE OF STRUCTURAL COMPOS ITES 




Structural 

Mechanics 

Geometric 

imperfections 

Post-buckling 

Damage 


Macro- 

Mechanics 

Local Damage 
at Stress 
Concentrations 


Minl- 

Mechanlcs 

Ply Fracture 
Delamination 


Micro - 
Mechanics 

Fiber Fracture 
Matrix " 
Interface " 


FractureO==ONDE 


What is the significance of damage or defect on the structural integrity? 
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Him III iiiiiiii II I I 


I 

CONCLUSION AND DISCUSSION / 

1. Methods of classical fracture mechanics can be viable crack growth criteria 
for matrix-dominated sublamlnate cracks, such as transverse cracking and 
ply delamlnatlon. 

2. The energy release rate Increases (generally) with the size of crack, which Is 
constrained by the physical dimension of the laminate's ply structure, such as 
the thickness of the piles. Thus, ply thickness has a profound Influence on the 
growth process of sublamlnate cracking, and even on the final failure of the 
laminate. 

3. The f Inlte-element/crack-closure technique Is a useful computational tool to 
simulate sublamlnate crack growth, provided the crack paths are relatively 
simple. 

4. The DIB (dllodobutane) enhanced X-radlography Is a useful NDE tool for monitoring 
the growth of sublamlnate cracks. 

5. For T300/934, the value of Is In the range 1.2 to 1.5 In-lb/ln^, but other 

experiments found It In the range 0.5 to 0.8 (refs. 3 and 6). (See sketch.) 


~ 1 . 2 - 1. 5 


This difference should be verified by additional experimentation. 

6. For the best fit of data, the critical material crack size a^ is taken as 
0.008 in. This should be verified by micromechanical investigation. Generally, 
a^ is a function of the size of inhomogeneity (fiber), the stiffness property 
surrounding the crack, and other processing variables. However, it is not known 
whether a^ might also depend on the ply structure of the laminate. 

7. In the study of delamlnatlon, all predictions were based on the total energy 
release rate G = Gj + Gxi + Gju for both mode-I and mixed -mode cracks. 
Experimental correlation indicates a need for a mixed-mode crack criterion, 
such as 



^l/^lc 


^II^^IIc 


^IIl'^^IIIc 


= 1 


There Is no firm evidence that Gjjj. and actually exist, and experi- 

mental means to measure them Independently has not been developed. 
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Some of the observed cracks were actually *planar cracks'. In the case of a 
crack front with a complicated contour, the degree to which a 2-D crack growth 
criterion is readily applicable is not known. 
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NOTCH STRENGTH OF COMPOSITES 


James M, Whitney 

Air Force Wright Aeronautical Laboratories 
Wright Patterson Air Force Base, Ohio 45433 



(0/± 45/90)25 T300/5208 


05 = 71.7 KSI 



^a = 0.12 IN 






STATISTICAL BEHAVIOR OP COMPOSITES 


WEI BULL STATISTICAL STRENGTH THEORY 

R (S)-EXP |-B(S)j 
BIS) -RISK-TO -BREAK 
TWO PARAMETER WEI BULL DISTRIBUTION 



TENSION: 
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R(S) 


BENDING 



PT; 


(ot + 2 ) \ 

Z|(«+ 1)2 ^$0^ 


H (CX + 1)2 

(CX+ 2 ) Vj^ 



3 PT: 



r('-T 

2 (o<+D^ \SqJ 


S. = S„ T2 (0< +1)2 


V< 


’ey 


BENDING OPT): 




1.52.®^= 15 
1.33, «= 25 


BENDING ( 4 PT): 



(o< + 1)2 
(o<+ 2) 



ex = 15 
oc= 25 
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THICKNESS EFFECTS 


TENSION: 



BENDING; ^ L (o^+ i)2 

(3 PT) " ° Tbh^ 

_ 




h2>hl 


END 

TAB 
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TAB 
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CIRCULAR HOLE APPROXniATE STRESSES 


Oy(X,0)=*^ 


CO 


{ I 3(t) -< Kt- 3)[5(-|)*-7(-|)*] j 


X >R 



CIRCULAR hole: 

POINT STRESS CRITCRION 



2 

[2 + «,*+3e,«-(K®-3)(5£,*-7€,«)] 



_R 

R +do 


AVERAGE STRESS CRITERION 

2 ('-g 2 > 


[ 2 - + ( K? - 3 )( • ) ] 



_R 

R + Oq 
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FINITE WIDTH CORRECTION - 


CIRCULAR HOLE 


=Y3 (2R/W) CTn 


Y3(2R/W)= 


2+(l-2R/wf 

3M-2R/W) 
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CENTER CRACK 


CTy(X.O)« 



X 


3 


K, X 



X>c 


POINT STRESS CRITERION 




c 

C +(1q 


AVERAGE STRESS CRITERION 



FINITE WIDTH CORRECTION - CENTER CRACK 

S “ Y^(2c/w)Ojj 

Yj^(2c/w) » 1 + 0.128(2c/w) - 0.288(2c/w)^ 
+ 1.52(2c/w)^ 

2c/w < 0.7, ERROR < 0.5% 
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[0/+45/90] SCOTCHPLY I002 


Vf = 0.50 
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MODIFIED THEORY 


R = HOLE RADIUS 

Rq = NORMALIZING FACTOR = UNIT LENGTH 
m = EXPONENTIAL PARAMETER 
C = NOTCH SENSITIVITY FACTOR 

^ ^ 2 

[(2 + )^ + 3\S - (K^ - 3) - 7A^i] 
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INFLUENCE OF STACKING SEQUENCE ON NOTCHED STRENGTH 



Log^O^ (inch««) 


SUPERHIPOSING NOTCHED STRENGTH 
MATERIAL 1 AND MATERIAL 2 



R ■ ARBITRARY RADIUS FOR MATERIAL 1 

R ■ RADIUS IN MATERIAL 2 PRODUCING SAME 

PERCENT STRENGTH REDUCTION 
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DEFINE EXPONENTIAL SHIFT PARAMETER 



EXPONENTIAL SHIFT PARAMETER 






COMPARISON BETWEEN WEIBULL THEORY AND MODIFIED THEORY 



- 6/2 < X £ 6/2 
QS.y 

-h/2 < z < h/2 



L 

[‘■m y)]“dY 


t{0, y) = 1 + l/2p* + 3/2p* 

- (K? - 3)(5/2p® - 7/2p«) 

where p = R/y 


o-A'/o-„ = [1 + 1/2^* + 3/2|« 

- (K? - 3)(5/2f« - 7/2D]-’ 


where f = R/(R + d„) and do = R"’/K 


Loading geometry - volume of numerical 
integration. 


In do = In R — In K 
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Shap* 

paramatar 

a 

Exponantlal 

paramatar 

m 

Notch 

aanaltivity 

K 

Corralation 

coafficlant* 

r 

3 

0.31 

3.9 

0.985 

5 

0.39 

5.5 

0.987 

10 

0.52 

9.7 

0.990 

12 

0.57 

10.8 

0.986 

15 

0.56 

14.0 

0.983 

20 

0.61 

16.8 

0.996 

30 

0.57 

27.8 

0.977 


‘ Corrvlatlan coafllclant from Mn lu«t aquara* fM ot linaar aquatian, 
la <(, - III la A - In K; minimum thfaa palMa (rndH) uaad tor aach givan a. 



WEIBULL SHAPE PARAMETER a 

Exponential parameter m as a function 
of shape parameter a. 



Notch sensitivity K as a function 
of shape parameter a. 
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COMMENT 


An interesting observation has been made (ref . 1) concerning the Whitney "point 
stress" or "average stress" criterion. In particular, the center-notched unidirec- 
tional laminate with no notch tip damage is shown to have a square root type stress 
distribution with an equivalent notch length which differs from the actual number 
of broken fibers by a small but constant amount, independent of the number of broken 
fibers. This is consistent with the assumptions leading to the "point stress" and 
"average stress" criterion. This is not the case if notch tip damage (in the foirm 
of matrix yielding or transverse broken fibers) is present, and it is shown that 
the damage produces stresses in the unbroken fibers which are less severe than a 
square root behavior, and also that an analogous equivalent notch length does not 
exist . 


REFERENCE 

1. Goree, J. G.; Dharani, L. R. ; and Jones, W. F.: Mathematical Modeling of 

Damage in Unidirectional Laminates. NASA CR-3453, 1981. 


James G. Goree 
Clemson University 


COMMENT 

The Weibull distribution function is based on a weakest link concept of failure, 
which in general is not applicable to composites. This concept underlies the size 
effect which would account for the difference between bending and uniform tensile 
strength in purely brittle materials. The difference between tensile and bending 
Strengths in composites is more likely related to the "gradient" effect, which 
allows localized failure and stress redistribution, as in the case of elastic- 
plastic materials. 

I. M. Daniel 

Illinois Institute of Technology 
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COMMENT 


The present point stress and average stress criteria relate the notched strength 
of a laminate to the average strength of a relatively long tensile coupon. Tests 
of notched specimens in which microstrain gages have been placed at or near the edges 
of the holes have measured strains much larger that those measured in an unnotched 
tensile coupon. Furthermore, orthotropic stress concentration analyses of failed 
notched laminates have also indicated that failure occurred at strains much larger 
than those experienced on tensile coupons with normal gage lengths. Earlier, both 
Hahn (ref. 1) and Wu (ref. 2) presented data that related fiber strength to gage 
length. This suggests that the high strains at the edge of a hole can be related 
to the very short length of fiber subjected to these strains. Since the length of 
fiber that is highly strained is proportional to the hole size, this would explain 
the higher strains to failure measured at the edge of smaller holes. 

Lockheed has attempted to correlate a series of tests of several laminates with 
holes ranging from 0.19 to 0.50 in. Although the average stress criterion correlated 
well with test results for hole sizes equal to or greater than 0.50 in., it over- 
estimated the laminate strength in the range of hole sizes from 0.19 to 0.38 in. It 
thus appears that we need a theory that is based on the mechanics of failure and is 
more generally applicable to the range of hole sizes and the varieties of laminates 
found in aircraft construction. 


REFERENCES 


1. Hahn, H. T. : Failure Mechanisms. Failure Analysis and Mechanisms of Failure of 

Fibrous Composite Structures, NASA CP-2278, August 1983, pp. 153-171. 
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COMPRESSION FAILURE OF 


COMPOSITE LAMINATES 


R. B. Pipes 

Center for Composite Materials 
College of Engineering 
University of Delaware 
Newark, Delaware 



This presentation will attempt to characterize the compressive behavior of 
Hercules AS-1/3501-6 graphite-epoxy composite. This involves: 

Studying the effect of varying specimen geometry on test results 

Determining end condition factors for IITRI fixture 

Determining transition region between buckling and compressive failure 

Defining failure modes and developing analytical models to describe 
and/or predict these modes 



Modified IITRI compression fixture. 
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angulor across width brooming 



ongulor ocross thickness biongular across thickness 
Failure modes 


dimension units - inches 

(centimeters) 



four specimens tested at each L/r ratio 



L/r 

Test coupons 
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lao 


Ultiinate stress (cr^) versus 
slenderness ratio (L/r) 
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^ 2 
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L/r 


End condition factor (1/k ) versus 
slenderness ratio (L/r) 


Maximum compressive strength occurred at 25 < L/r < 30 
End condition factor increased as L/r increased 
Failure modes change with varying L/r 


Conclusions 



Induce fatigue damage in the vicinity of a circular hole in a composite 
laminate through compressive loading 

Observe and characterize the nature and extent of damage 

Determine the influence of Induced damage upon residual compressive 
failure 


Objectives 


(R=0.1, frequency of 10 Hertz, 
ciinbient conditions) 


Specimen 

S-Level 

Cycles 

1 

(Static) 

— 

2 

(Static) 

— 

3 

(Static) 

— 

4 

0.50 

IxlO® 

5 

0.50 

1x10® 

6 

0.50 

IxlO® 

7 

0.60 

IxlO® 

8 

0.60 

IxlO® 

9 

0.60 

1x10® 

10 

0.60 

5x10® 

11 

0.60 

5x10® 

12 

0.60 

5x10® 


Fatigue program 



<=> 


A: [02/±45]53 

[0/0/45/-45/0/0/45/-45/0/0/45/-45/0/0/45/-45/0/0/45/-45] g 

B: [0/45/0/-45]3g <=> 

[0/45/0/-45/0/45/0/-45/0/45/0/-45/0/45/0/-45/0/45/0/-45]^ 

C: [0/±45/90]5g <=> 

[0/45/-45/90/0/45/-45/90/0/45/-45/90/0/45/-45/90/0/45/-45/90] 

D: [90/0/+45153 <=> 

[90/0/45/-45/90/0/45/-45/90/0/45/-45/90/0/45/-45/90/0/45/-45] 

Laminate configurations 



Modified IITRI compression fixture 



r 

MICRO-MEASUREMENTS 

STRAIN OAOE A< EA-OC-I2SBZ-3SO 
STRAIN «A8E S' EA-06-I2SAC-3S0 

Notched fatigue specimen 


(Lood Direction) 



AXIAL SECTION 
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Cracks 


- 45 ® 

Cracks 
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t 


45 ® 
Cracks 


Fatigue crack pattern 



Cracked laminae by quadrant 
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Q'f f f 


Antisymmetric distritution of 
shearing stress T 
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Laminate "D"; [ 90/0/±45]5g; axial section (50x) 



Near-surface delamination growth - laminate "A" 
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Near-surface delanlnation growth - laminate "C" 


1 

n 

Lood 4 
Direction 

n 

n 

Q 

Near-surface delamination growth - 
laminate "D" 
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Two modes of compressive failure observed: 

Diagonal shear - predominant mode in fiber-dominated laminates 
Net compression - predominant mode in quasi-isotropic laminates 

Both failure modes are characterized by local instability of individual 
lamina or small lamina subgroupings 

Mode of failure is related to the nature of the specimen delamination and 
the fatigue- Induced cracking 

Angle plies crack in transverse tension during fatigue cycling with the 
cracked layers determined by quadrant location about the hole 

Primary direction of near-surface delamlnatlon progression is dependent 
upon the fiber direction of the laminate surface ply 

Laminate stacking sequence affects the nature of cracking as well as the 
predominant failure mode 

Failure mechanisms are not significantly different for the two resin systems 


Conclusions 


• DEVELOPMENT OF AN EXPRESSION FOR THE BUCKLING OF THE 
DEBONDED REGION ABOVE AN IMPLANTED FLAW. 


• VERIFICATION OF THE BUCKLING ANALYSIS THROUGH THE USE 
OF EXPERIMENTAL TESTING. 


• DEVELOPMENT OF AN UNDERSTANDING, THROUGH EXPERIMENTAL 
TESTING, OF THE INFLUENCE OF IMPLANTED FLAWS ON THE 
COMPRESSIVE PERFORMANCE OF SEVERAL GRAPHITE/EPOXY 
LAMINATES. 


Objectives 
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• DEBOND DEFECTS MAY ALLOW FOR PREMATURE COMPRESSIVE 
FAILURE OF LAMINATES. 

• INSTABILITY ANALYSIS MUST TREAT ASYMMETRIC LAMINATES 
AND CONSIDER TRANSVERSE SHEAR DEFORMATIONS. 

• BOUNDARY GEOMETRY INFLUENCES BUCKLING LOAD MORE THAN 
LAMINATE ASYMMETRY. 

• LAMINATE GEOMETRY CFIBER ORIENTATION AND STACKING 
SEQUENCY) STRONGLY INFLUENCES DEFECT CRITICALITY. 

• DEFECT INSTABILITY PRECEDES INTERLAMINAR CRACK PROPA- 
GATION AND ULTIMATE FAILURE. 

• SIMPLIFIED BUCKLING ANALYSES MODEL INSTABILITY 
INITIATION. 

• STRENGTH LOSSES OF 50 PERCENT WERE OBSERVED. 


Debond defect studies. 


FOR AN ISOTROPIC MATERIAL 

D = 


Eh- 


12(1-v2) 


THE POTENTIAL ENERGY FOR A PLATE WITH A SINGLE IN-PLANE 
LOAD IS 

Vz] ^ f ^ ” 2(l-v)(w w ~ ) ] 

J \ XX yy XX yy xy' * 


J' 


+ tN^w^] JdS = Constant 


INSERTING OUR ASSUMED DEFLECTION SHAPE 


w(r) = a[i- (i)^ 


PERFORMING THE INTEGRATION AND THE VARIATION YIELDS THE 
BUCKLING CONDITION FOR THE ISOTROPIC CASE 


N 


-32D 


crit 
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INSERTING THE DEFLECTION SHAPE AND PERFORMING THE INTEGRATION 
YIELDS 


N 


-24 

*crit 




BY COMPARING ELEMENTS OF THE BENDING STIFFNESS MATRIX WE SEE 
THAT FOR AN ISOTROPIC MATERIAL 


“ “22 - “ 


“l2 = 

“66 ■ “ 


INSERTING THESE INTO THE ORTHOTROPIC RESULT PRODUCES THE 
ISOTROPIC CASE. 


N 


- 24 

"crit 


D + D ^ -oD + 2 (1 -v)/2d '| _ -^D 


2 


J 



Sandwich construction 



















LOAD (KQ) 



1000 2000 3000 4000 5000 6000 7000 


STRAIN (/tin/in) 

Load/strain data for sample Fg-1-2 
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Failure of a sandwich beam composite 
with a 1.5~in. circular defect located 
at the near surface 


LOAD (LBS) 














DEFECT DIAMETER (in.) 


concwitratMl kMd 
ooinpr«uion tollur«'*~ 


I.O(r 

.8 ■ 

.6 

.6 

.4 


Defect influence 
on (tflSjzs 
Lominote 


^eonioreuion 
buckling failure! 


I no defect 
I Teflon Film 
Kopton Bog 


10 


20 


30 


40 


50 


DEFECT DIAMETER (mm) 


40 in 
>c 




to 


Influence of Kapton bag and Teflon film 
defects on residual strength of the 
[±45]2s laminate 


DEFECT DIAMETER (in.) 



Effect of laminate thickness on residual 
strength of the [0/±45]ns laminate 




1.0-in. circular near-surface 
delamination prior to testing 


1.0-in. circular near-surface 
delamination after buckling 



1.0-in. circular near-surface 
delamination after 2000-kg 
applied load 






Load compressive failure in graphite/epoxy (200x) 


• STUDY VARIOUS DEFECT CONFIGURATIONS TO EVALUATE THE 
SENSITIVITY OF THE LAMINATE TO DEFECT MATERIALS AND 
CONSTRUCTION. 


• EXTEND THE 
THROUGHOUT 


WORK TO DEFECTS LOCATED AT OTHER POSITIONS 
THE LAMINATE THICKNESS. 


• CONDUCT TESTING DESIGNED TO PRODUCE BUCKLING IN AN EFFORT 
TO VERIFY THE BUCKLING RESULT. 


• PRODUCE A SOLUTION TO THE TWO-DIMENSIONAL FRACTURE PROBLEM 
WHICH STEMS FROM THE BUCKLED GEOMETRY. 


Recommendations for future work 


LAMINATE RESIDUAL STRENGTH IS A FUNCTION OF FAILURE MODE 
FOR THE HIGH STRENGTH CONFIGURATIONS BUT NO INFLUENCE ON 
STRENGTH WAS SEEN FOR THE LOW STRENGTH LAMINATES. 


CERTAIN LAMINATES FAIL BY A CHARACTERISTIC MECHANISM AND 
THIS MECHANISM MAY BE INDEPENDANT OF THE PRESENCE OF 
DELAMINATION FLAWS. 


IMPLANTED FLAW CONSTRUCTION 
ON COMPRESSIVE PERFORMANCE. 
OF SELECTING AN APPROPRIATE 
MODELING A DELAM I NAT I ON . 


CAN HAVE A SUBSTANTIAL INFLUENCE 
THIS HIGHLIGHTS THE NECESSITY 
DEFECT CONSTRUCTION WHEN 


THE OCCURRENCE OF BUCKLING PLAYS AN IMPORTANT ROLE IN THE 
FAILURE PROCESS BECAUSE DELAMINATION GROWTH OCCURS WITH 
THE ONSET OF DEFECT INSTABILITY. 


THE BUCKLING ANALYSIS IS 
BUCKLING LOAD FOR THE TH 
OF THE APPROXIMATION IS 


USEFUL IN PREDI 
ICKER LAMINATES 
LOST FOR THE THI 


CTING THE DEFECT 
BUT THE ACCURACY 
N 6 PLY RESULTS. 
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COMMENT 


The fact that It Is difficult to measure compressive strength In composites 
does not mean that there Is no such property. Even If It were Impossible to 
measure It, we would have to define It In order to deal with structures subjected 
to combined states of stress. It can be defined macroscoplcally as a structural 
failure of the particular lamina or laminate. 

I. M. Daniel 

Illinois Institute of Technology 
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THE IN SITU TRANSVERSE LAMINA STRENGTH 


OF COMPOSITE LAMINATES 


Donald L. Flaggs 
Applied Mechanics Laboratory 
Lockheed Palo Alto Research Laboratory 
Palo Alto, California 


The objective of the work reported In this presentation is to determine the 
in situ transverse strength of a lamina within a composite laminate. From a 
fracture mechanics standpoint, in situ strength may be viewed as constrained cracking 
that has been shown to be a function of both lamina thickness and the stiffness of 
adjacent plies that serve to constrain the cracking process. From an engineering 
point of view, however, constrained cracking can be perceived as an apparent increase 
in lamina strength. With the growing need to design more highly loaded composite 
structures, the concept of in situ strength may prove to be a viable means of 
increasing the design allowables of current and future composite material systems. 

As a means of Introducing the current research, some earlier results for cross- 
plied graphite/epoxy and glass/epoxy laminates and [±25/90j^]g graphite/epoxy 
laminates are presented. In this work the strain at onset of transverse cracking 
is seen to be a strong function of 90°-laminae thickness. In the figures that 
follow, E represents the total or measured strain and E represents the elastic 
strain (e = £ - aAT) . Data are also presented indicating that for an off-axis 
loading angle greater than about 5°, matrix cracking is the initial lamina failure 
mode. Experimental results from this research are given, and the 02 stress (in situ 
transverse strength) in the 90° laminae of the [±0/9On]s laminates is plotted, based 
on laminated plate theory predictions, at the onset of transverse cracking. The 
separate contributions of applied and residual thermal stresses on the predicted 
in situ strength are Illustrated, and the effect of the different constraining 
±9 laminae on in situ strength (which has been termed an interaction effect) is 
examined . 

A simplified one-dimensional analytical model is presented that is used to 
predict the strain at onset of transverse cracking. While it is accurate only for 
the most constrained cases, the model is important in that the predicted failure 
strain is seen to be a function of a lamina's thickness d and of the extensional 
stiffness bEo of the adjacent laminae that constrain crack propagation in the 
90° laminae. 

The possibility is investigated of using a Weibull brittle failure model, in 
which the probability of failure is related to the volume of stressed material, to 
predict the in situ strength. Using a shape parameter determined for a material 
system (T300/5208) similar to that used in the current investigation (T300/934), it 
was seen that the Weibull model underestimated the observed thickness effects. 
Additionally, there is no mechanism in the model to account for the interaction effect 
on in situ strength. However, it was interesting to note that by assuming that the 
shape parameter was not a material parameter, we could quite accurately fit the 
experimental behavior. 

In closing, perhaps the most important conclusion that can be drawn from these 
results is that strength predictions for failure modes associated with matrix cracking 
in composite laminates should not be based upon the concept that strength is an 
Intrinsic material property. Work is currently under way to investigate the possi- 
bility of integrating these results into existing strength theories. 
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Fracture modes in off-axis tensile loading 
of ^ 03 ^ Mod mo r- 1 /epoxy laminates 


Mode 

Load angle 

Longitudinal tension 

0 < 0 < 5° 

Intralamina shear 

5 < e < 20° 

Mixed mode 

20 < 0 < 45° 

Transverse tension 

45 < e < 90 ° 


Laminate geometry of [to/ 90n sj specimens 
z 
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laminates; 


Ea 

.798 

.627 

.495 

.461 


Inner ply thickneas (m) 


Calculated transverse stress in 90° plies of [±25/90n]s 
T300 / 934 laminates at transverse cracking onset 
{From ref. 2 \ 


N 

YT 

e^(%) 

.5 

106 (15,370) 

## 

.650 

1 

111 (16,095) 

.665 

2 

83.9(12,165) 

.415 

3 

76.5(11,092) 

.360 


71.2(10,324) 

.325 

6 

69.4(10,063) 

.330 

8 

69.2(10,034) 

.345 


*MPa (Psi) 

Ultimate failure strain 



Mechanical strain at onset of transverse cracking in |0/90/0]. 
glass fIber/epoxy laminates; tg - 03mm {From ref. 3 } ‘ 





Transverse ply thickness 2d (mm) 


Experimental results 


Specimen 
Fami ly 

n 

Number of 
Specimens 

Ex (MSI) 

(02/90n)s 

1 

4 

13.73 


2 

4 

10.92 


4 

4 

7.51 


8 

5 

5.29 

(+30/90n)s 

1 

4 

7.73 

2 

4 

6.41 


4 

3 

4.82 


8 

4 

3.43 

(+60/90n)s 

1 

2 

2.08 

2 

4 

1.95 


4 

4 

1.86 


8 

4 

1.73 



Ex(X) 


Nx(’b/1n) 

.555 

(.549-. 562 

) 2428.9 (2390.4-2491.0) 

.319 

(.287-. 378 

) 1456.3 (1381.5-1676.5) 

.236 

(.221-. 242' 

) 1158.4 

(1083.7-1184.6) 

.188 

(.173-.225; 

) 1034.7 

(894.8-1289.9) 

.498 

(.489-. 506] 

1 1219.9 

(1206.3-1259.8) 

.386 

(.371-. 431] 

1 1036.3 

(995.5-1150) 

.293 

.267-. 338 

1 865.4 

(799.2- 997.5) 

.242 

(.231-. 254] 

1 860.5 

(796.9- 898.8) 

.592 

(.568-. 615) 

1 387.1 

(374.4- 399.7) 

.553 

(.475-. 596 

449.5 

(399.8- 474.5 

.453 

(.405-. 503) 

511.1 

(471.7- 524.3) 

.405 

(.387-. 416) 

717.2 

(701.8- 742.2) 
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n 



n 
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Analytical model for predicting onset of transverse cracking 


Griffith energy balance used In conjunction with a one-dimensional shear lag 
model to account for load transfer resulting from a Mode I crack spanning 
the 90 laminae 


Griffith energy balance; AW - AU ■= 2Ygp (per unit area) 


One-dimensional shear lag model: Ao = ^ E 2 EXP(-(ji^ X) 


Resulting expression for in situ transverse strength; 

Ik 


TC 


b Egm 


(b+d)d“Ejj 


°23 “le 


Weibull material (7) 


Pj(v) = exp 




where a = shape parameter 
0 ^ = strength cut-off 
Ojj = scale parameter 




1/a 


®2-°u 


I/O 


For T300/5208, a ■ 7.68; ^ ~ 

Now. consider the fainily (assuming o^ = 0): 


A - - 1 89 ? 

Og “ 8,350 
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Summary of in situ transverse strengths 


Specimen 


exp. 

_pred. 


Fami ly 

n 

^2(5!) 

E2{%) 

YT(psi) 

(Oz/90n)s 

1 

.953 

.950 

15,847 


2 

i700 

.662 

11,751 


4 

.590 

.465 

9,944 


8 

.497 

.316 

8,350 

(+30/90n)s 

1 

.879 

.839 

14,154 


2 

.739 

.579 

12,140 


4 

.604 

.397 

9,869 


8 

.496 

.269 

8,063 

(+60/90n)s 

1 

.721 

.815 

12,503 


2 

.666 

.545 

10,909 


4 

.540 

.352 

8,540 


8 

.463 

.223 

7,234 


Experimental versus predicted in situ strength ratios 

YTi ” \vj 

Specimen Family a Vol . Ratio, 'V?' Experimental , 7T^ Predicted, 


{02/90n)s 

3.245 

2:1 

1.349 

1.238 



4:1 

1.594 

1.533 



8:1 

1.898 

1.898 


(+ 30 / 90 n)s 

3.697 

2:1 

1.166 

1.206 



4:1 

1.434 

1.455 



8:1 

1.755 

1.755 


(+ 60 / 90 n): 

3.802 

2:1 

1.146 

1.200 


4:1 

1.464 

1.440 



8:1 

1.728 

1.728 
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Summary 


Design Implications: 

0 Individual plies should be dispersed rather than stacked when practical 

0 Laminate failure analyses should include these in situ strength effects 
AS soon as a methodology has been established for quantifying them 
FOR GENERAL LOADING CONDITIONS 

Work in Progress: 

0 Experimental determination of effect of on onset of transverse 
CRACKING 

0 Analytical simulation of transverse cracking using a generalized 

plane strain FE model 

1) Better analytical representation of cracking phenomenon 

2) Can account for mixed mode effects 

3) Analytical derive design curves 



EFFECTS OF DEFECTS IN COMPOSITE 


STRUCTURES 


G. P. Sendeckyj 
Structural Integrity Branch 
Flight Dynamics Laboratory 
Air Force Wright Aeronautical Laboratories 
Wrlght-Patterson Air Force Base, Ohio 



DESIGN CRITERIA 


FIRST PLY FAILURE 

• DESIGN ULTIMATE STRAIN (DUS) IS STRAIN FOR FIRST PLY FAILURE 
STATE OF THE ART 

• DUS IS FAILURE STRAIN OF LAMINATE WITH 0.25-in. HOLE 
POSSIBLE 

• DUS IS FAILURE STRAIN OF LAMINATE WITH LOW-ENERGY IMPACT 
DAMAGE 

• BUFFER AND SOFTENING STRIPS USED 
UPPER BOUND 

• DUS IS FAILURE STRAIN OF LAMINATE 

• STRESS CONCENTRATIONS DESIGNED AROUND LEAST CONSERVATIVE 


PREPREG DEFECTS 

HOLLOW FI BERS 

EXCESSIVE VARIABILITY IN FIBER PROPERTIES 
RESIN-STARVED OR FIBER-STARVED AREAS 
WRINKLES, WAVINESS, MISCOLLIMATION 
FOREIGN PARTICLES, CONTAMINATION 
PILLS AND FUZZ BALLS 
NONUNIFORM AGGLOMERATION OF HARDENER 
PREPREG OUT OF SPECS. 
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DEFECTS IN LAMINATES 

HOLLOW FIBERS DELAMINATIONS 

FIBER BREAKS PLY GAPS 

EXCESSIVE POROSITY, VOIDS 

RESIN-RICH AND RESIN-STARVED AREAS 

FIBER WAV I NESS, WRINKLES, MISCOLLIMATION 

FOREIGN PARTICLES, CONTAMINATION, INCLUSIONS 

INCOMPLETE AND/OR VARIABLE CURE 

WRONG STACKING SEQUENCE 

DENTS, TOOL IMPRESSIONS, SCRATCHES 

LAMINATE POROSITY 

STUDIED EXTENSIVELY 

MATRIX DOMINATED PROPERTIES DEGRADED (DELAMINATION 
NOT INCLUDED) 

5% STRENGTH REDUCTION FOR 1% POROSITY 
50% LIFE REDUCTION FOR 1% POROSITY 

FIBER-DOMINATED PROPERTIES NOT AFFECTED 

DELAMINATION GROWTH AFFECTED - NOT WELL DOCUMENTED 

MOISTURE ABSORPTION 

EQUILIBRIUM MOISTURE LEVELS INCREASED 
AGGRAVATES THERMAL SPIKE PHENOMENOM 



EFFECT OF PLY GAP DEFECT 


(REF. 1) 

[(0/45/90/ -45)g]2 LAMINATE 

- 16.9% STRENGTH REDUCTION FOR GAP(S) IN 90 PLIES 

- 8.7% REDUCTION FOR GAP(S) IN 0 PLIES 
[ (0/45/0/ -45/0)g]2 LAMINATE 

- 6.5% STRENGTH REDUCTION FOR GAPS IN 0 PLIES 


(REF. 2 ) 

[(0/±45/90)^]2 laminate 

- 12.8% STRENGTH REDUCTION FOR GAPS IN OUTER 45 PLIES 
[(0/-/45/0)532 LAMINATE 

- 6.2% STRENGTH REDUCTION FOR GAPS IN OUTER 45 PLIES 
DEFECT CRITICALITY - BENIGN FOR DESIGN ULTIMATE STRAIN 0.7% 

EFFECT OF PLY WAVINESS DEFECT 


(REF. 3 ) 

SURFACE 0 PLY WAVINESS IN [(0/45/90/ -45 ) 512 LAMINATE 

STATIC TENSILE STRENGTH REDUCTION 

• 10% FOR SLIGHT WAVINESS 

• 25% FOR EXTREME WAV I NESS 

FATIGUE LIFE REDUCTION 

• AT LEAST A FACTOR OF 10 

• CONSISTENT WITH STATIC STRENGTH REDUCTION 

DEECT CANNOT BE FOUND BY STANDARD NDE 

STRENGTH LOSS CAN BE PREDICTED BY ASSUMING LOSS OF 
LOAD CARRYING CAPACITY DUE TO THE WAVINESS 

DEECT CRITICALITY - INSUFFICIENT DATA FOR ACCURATE ASSESSMENT 

• SHOULD BE BENIGN FOR DESIGN ULTIMATE STRAINS 0.7% 



MACHINING DEFECTS 


EDGE DELAMINATIONS 
OVERSIZE HOLES 
UNDERSIZE HOLES 
TILTED HOLES 


EDGE NOTCHES AND SURFACE NOTCHES 
HEAT-DAMAGED MACHINED EDGES 
FIBER BREAK-OUT ON HOLE EXIT SIDE 
OUT-OF-ROUND HOLES 


TILTED COUNTERSINKS IMPROPER DEPTH OF COUNTERSINKS 

DENTS, FIBER BREAKING FROM IMPACT 
TEAROUT OR PULL-THROUGH IN COUNTERSINKS 


EFFECT OF SURFACE NOTCHES 

EXPERIMENTAL DATA 

STATIC STRENGTH REDUCED UP TO 50% 

LOCAL DELAMINATION AT NOTCH 
FATIGUE LOADING REDUCES STRESS CONCENTRATION 
RESIDUAL STRENGTH HIGHER THAN STATIC STRENGTH 
DATA AVAILABLE FOR VARIOUS STACKING SEQUENCES 

ANALYSIS 

ECCENTRIC BEAM MODEL PREDICTS STRENGTH REDUCTION 

STRENGTH REDUCTION IS SMALL FOR SIZES EXPECTED 
IN SERVICE 

DEFECT CRITICALITY 

NOT CRITICAL FOR DUS <0.7 



BOLTED ASSEMBLY DEFECTS 


OVERTORQUED FASTENERS IMPROPER FASTENER SEATING 

MISSING FASTENERS FASTENER INSTALLATION DAMAGE 

OVERSIZED AND UNDERSIZED FASTENER 

BONDING DEFECTS 

ADHESIVE POROSITY 
MISCURE 

ADHESIVE-STARVED AREAS 


IMPROPER SURFACE PREPARATION 
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EXPERIMENTAL METHODS FOR 
IDENTIFYING FAILURE MECHANISMS 


I. M. Daniel 

Illinois Institute of Technology 
Chicago , Illinois 




I 


Fiber composite 



Single ply 


Microscopic 

Matrix failure (tensile^ compressive^ shear) 
Bond failure 
Fiber failure 

Miniscopic 

First-ply failure 
Lamina failure criteria 

Macroscopic 

Laminate failure criteria 

Scales of observation of failure 
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1. Photoelastic (2-D^ 3-D/ microphotoelastiC/ dynamic/ 

BIREFRINGENT COMPOSITES/ BIREFRINGENT 

coatings) 

2. Moire 

3. Strain Gages 

4. Interferometric and Holographic Methods 

5. Nondestructive Evaluation (Ultrasonics/ acoustic 

EMISSION/ X-RAY^ thermography) 

6. Fractography 


Experimental methods 


STRESS-OPTIC LAW 


where 


a, - O2 “ difference of "secondary" 
principal stresses 


n = fringe order 


f = material fringe value 
(constant for material) 


t = specimen thickness 


Photoelastic method 


^ R 


where 


S = gage factor (function of alloy and backing 
S of gage) 



(6-a) S AT + yAT 

O 


a = thermal coefficient of expansion of gage material 


B = thermal coefficient of expansion of base material 
Y = coefficient of resistivity of gage material 


Electrical resistance strain gages 



59 lines per inch 


5 I lines 
per inch 


Fringes due to 
rotation alone 


Fringes due to 
difference in 
pitch alone 


Fringes due to 
combination of 
rotation and 
difference in pitch 


Mechanism of formation of Moire fringes 
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Strain-optic law; 



where 

“ strain fringe value 
N « fringe order 
h » coating thickness 

c,s «* refer to coating and specimen, respectively. 
Conditions at boundary: 

At interface between coating and specimen, 

_ s _ s s 

C22 ^22 ■''12 ^11 

At top surface of coating, 

c _ c c _ c s 
^22 ^11 ^ ^11 

Principal strain along boundary, 

Nf 1 

= g . 1 

^11 ~2Tr 1 j. c 
1 + V 

Photoelastic coating method (refs. 1 to 4) 




Hologram recording 


Reconstruction 


Holographic processes 


' OBJECT 
/| (viriuol image) 


hologram 
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T“ TRANSDUCER 
S- SPECIMEN 


Ultrasonic pulse echo method 


t t t t 


INEFFECTIVE 

LENGTH 



Failure model of unidirectional composite under longitudinal tension (ref. 5) 






Sequence of photographs showing distribution of fiber breaks 
in unidirectional composite under longitudinal tension 









Fixture for dynamic tensile loading of composite models 



Transient Isochromatlc fringe patterns In a glass-plastic 
composite model under dynamic tension (Camera speed : 
200,000 frames per second) 






[+ 45 ^/ 0 ^/ 0 ^] 


[0^/+A5^/0"] 


Characteristic failure patterns of three graphite/S-glass/high 
modulus epoxy specimens under uniaxial tensile loading 


Isochromatlc fringe patterns around hole in [0/±45/0/90] boron/ 
epoxy specimen for applied uniaxial stresses of 166 MPa (24.0 
ksi), 225 MPa (32.6 ksi) , and 293 MPa (42.4 ksl) 





Typical failure pattern around hole in [0/±45/0/90]^ boron/epoxy 
specimen under uniaxial tensile loading 



198 MPa (28.7 ksi) 



206 MPa (29.9 ksi) 


206 MPa (29.9 ksi) 


210 MPa (30.4 ksi) 




Sequence of Moire fringe patterns corresponding to vertical 

displacements in [0/±45/0/90]g glass/epoxy specimen at various 
applied uniaxial stresses 



STRENOTH REDUCTION 


ExperlmenCal 


Average Stress Criterion 


HOLE RADIUS, a 


Strength reduction as a function of hole radius for [02/±45]2 
graphite/epoxy plates with circular holes under uniaxial 
tensile loading 



ApproxiinaCe Stress Distribution 


<Jy (x,0) - U + 7 p'^ + + I (k^- 3) (5p’^ - 7p-®)] 

o ■ far field stress 
o 


*■ anisotronic stress concentration factor 


Strength Reduction Ratio 


_ ? 

S„ (1+i;) [2 + (k - 3) 


- characteristic length dimension 

Syy, - strengths of notched and unnotched laminates, respectively. 


Strength reduction of unlaxlally loaded composite plate with hole 
according to average stress criterion 





76 MPa (11.0 ksl) Oy = 102 MPa (14.8 ksi) Oy = 128 MPa (18.5 ksi) 


Moire fringe patterns around crack in glass/epoxy composites 10/90/0/90], 


at three levels of applied stress 




121 MPa (18 ksi) 


162 MPa (23 ksi) 


202 MPa (29 ksi) 


222 MPa (32 ksi) 


243 MPa (35 ksi) 


273 MPa (40 ksi) 


Isochromatic fringe patterns in photoelastic coating around 1.27-cm 
(0.50 in.) crack of [0/±45/90]g graphlte/epoxy specimen at various 
levels of applied stress 






166 MPa (24 ksi) 


202 MPa (29 ksi) 




237 MPa (3'; ksi) 




261 MPa (38 ksi) 


285 MPa (41 ksi) 


293 MPa (42 ksi) 


Isochromatic fringe patterns in photoelastic coating of [0/±45/90]g 
graphite/epoxy specimen with 2 , 54-cm-diameter (1 in.) hole under 
equal biaxial tensile loading (Far-field biaxial stress marked) 


0* 90* 

-e - 22.5*. 67.5* 



_ 

’xx <3 
- 100 ^ 
in 


FRINGE ORDER, n 


STRAIN, 10 e 

Fringe order and circumferential strain at two locations on the hole 
boundary for [0/±45/90jg graphite/epoxy specimen with 2. 54-cm-diameter 
(1 in.) hole under equal biaxial loading 




STRENGTH REDUCTION, 


Failure pattern in [0/145/90] g graphite/epoxy specimen with 1.91-cm- 
diameter (0.75 in.) hole under equal biaxial tensile loading 



Strength reduction as a function of hole radius for [0/±45/90Js graphite/ 
epoxy plates with circular holes under 1:1 biaxial tensile loading 



CTxx = 152 MPa 
(22 ksi) 


202 MPa 


(29 ksl) 


Oxx = 253 MPa 
(37 ksi) 


Moire fringe patterns around crack in uniaxially loaded [0/±45/90]g graphite/ 
epoxy specimen for three levels of applied stress 




Crack opening displacement and far-field strain for [0/±45/90]g 

graphite/epoxy specimen with a 1.27-cm (0.50 in.) horizontal crack 


■f 


t 


\ 


*1 


4 I 


>1 I 




a 


Failure patterns in [02/145] g graphite/epoxy specimens with holes of 
various sizes under uniaxial tension (Hole diameters are 2.54 cm 
(1 in.), 1.91 cm (0.75 in.), 1.27 cm (0.50 in.), and 0.64 cm (0.25 in.) 


audible failure 
audible failure 

visible delar>ination on vertical axis 



79 nPa (11. A X lO” psi) 

0.67 

3.31 

AlA MPa (60 ksiV 


STRAIN. lO^f 


Vertical strains along horizontal axis of [02/±45]g graphite/epoxy specimen 
with 1 . 91-cm-diameter (0.75 in.) hole under uniaxial tensile loading 



Specimen no . 43 
[±45/02/0] g 


Specimen no. 45; Bo/E 
[02/±45/0]g 


Failure patterns of boron-epoxy tensile panels with holes 


[O2/+45/O] 


[+45/0. /O] 


Isochromatic fringe patterns in photoelastic coating around hole in boron/ 
epoxy specimens of two different stacking sequences (a = 392 MPa (56.8 
ksi)) ^ 



[0/90/0/90]^ 


, [45/90/0/-45]g 


r±45/0/±45]g 


[±45/±45]g 


Failure patterns of boron-epoxy panels with holes of various laminate constructions 




[±45/±45]g; Oy = 77 MPa (11.1 ksi) 




[0/90/0/90]g; (Jy = 170 MPa (24.6 ksi) 

Isochromatic fringe patterns in photoelastic coating around hole in boron/epoxy 
specimens 
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Failure patterns in uniaxially loaded [0/±45/90]g graphite/epoxy plates with 
cracks of various lengths (Crack lengths are 0,64 cm (0.35 in.), 1.27 cm 
(0,50 in.), 1.91 cm; (0.75 in,), and 2.54 cm (1.00 in.) 



Strength reduction as a function of notch size for [0/±45/90]g graphite/epoxy 
plates with circular holes and horizontal cracks under uniaxial tensile loading 
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50 


30-^ 

IB L3O 


S /it (•+■„} 


20-1 

^cr X20 


15 (n) 


HALF CRACK LQUmi. a 


Critical stress intensity factor as a function of crack length for [0/±45/90] 
graphite/ epoxy plates with horizontal cracks under uniaxial tensile loading 


C2 ” IwJ 


”v 1 + k - (I-Wcooze 
" (l-k)sin5s 


Stress transformations of the far— field biaxial state of stress around a crack 
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Biaxial loading of [0/±45/90]g graphite/epoxy specimens with cracks 




Isochromatic fringe patterns in photoelastic coating around 1.27-cm (0.5 in.) 
crack in (0/±45/90)g graphite/epoxy specimen under biaxial loading - 
CTyy = 2.0302^^ at 30 deg with crack direction 


335 








Front surface 



Back surface as viewed from 
the front through the specimen 
(ref. 6) 


Thermally induced holographic fringe patterns in fatigue-loaded [ 0/145/90] g 
graphite/epoxy specimen with circular hole 



TBE enhanced X-ray photographs showing fatigue-induced damage in [ (0/±45/90)g]2 
graphite/epoxy specimens 
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[0z/±45]2s [(0/±45/90)s] 


Pen-lift scan Analog scan (normal) Analog scan (perspective) 

Ultrasonic C-scans of [0/±45/90)g]2 graphite/epoxy specimen with a film patch 


Q INITIAL FLAW 
P O LIFETIMES 
DB 1/2 LIFETIME 
® 1 LIFETIME 
^ I >4 LIFETIMES 

0 2 LIFETIMES 
H 2 I 4 LIFETIMES 
DH 3 LIFETIMES 

1 3 '.4 LIFETIMES 
■ 4 LIFETIMES 


Flaw growth under spectrum fatigue loading in [ (0/±45/90)g]2 and [02/±45]2s 

graphite/epoxy specimens with four types of initial flaws (Ambient environment) 
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GENERAL 

DISCUSSION 



Ij 


DISCUSSION OF MATRIX CRACKING IN 
GRAPHITE/ EPOXY COMPOSITE MATERIALS 

Walter L. Bradley 
Texas A&M University 


Recent work on the epoxy 3502 with and without carbon fiber reinforcement has 
demonstrated that reinforcement of a weak brittle material with a strong brittle 
material enhances the fracture toughness, even when cracking is parallel to the 
fibers. Neat material compact tension coupons indicated that = 75 J/m^. 

Delamlnatlon cracking of ASl/3502 gave a toughness of 150 J/m^, whereas transverse 
cracking parallel to the fibers gave a toughness of 225 J/m^. In delamination 
cracking, occasional crack tip Interaction with the fibers results in either fiber 
breakage or crack arrest, requiring renucleation in the adjacent matrix material. 
Coalescence of the existing crack and the new crack gives a serrated or scalloped 
appearance to the fracture surface. For transverse cracking in a 90° compact tension 
specimen, the resin path through which the cracking occurs is quite "lean" compared 
to the resin-"rich" region between plies which is cracked during delamination. Thus, 
the incidence of crack tip interaction with the fiber is enhanced and the energy 
required to propagate the crack increases accordingly. 

Finally, mixed-mode delamination fracture appears to occur microscopically on 
the plane of principal normal stress, thus increasing the angle between the delami- 
nation plane and the plane of principal normal stress. The result is that crack 
propagation is continually being redirected into the fibers, resulting in either 
fiber breakage or crack arrest, with the need for renucleation in adjacent matrix 
material. As was discussed previously, this has the effect of increasing the energy 
required to propagate the crack. For G j j/ (G j + Gjj) increasing from 0 to 0.35, 
Gj^jj = Gj + Gj-j Was observed to increase from 150 J/m^ to 550 J/m^. 
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DISCUSSION OF COIIPOSITE STRUCTURES 


John N. Dickson 
Lockheed -Georgia Company 


Composite Structures , Are We Ready ? 


SECONDARY STRUCTURE 

OH SURE! 

CONTROL SURFACES 

0. K. I 

HORIZONTAL OR VERTICAL TAIL 

MAYBE! 

WING OR PRESSURIZED FUSELAGE 

OH NO! 


Composite Primary Structure 

PRINCIPAL DESIGN DRIVERS 

• DAMAGE TOLERANCE 

• OUT-OF-PLANE DEFORMATIONS 


• JOINTS AND ATTACHMENTS 



Damage Tolerance Requirements 
Metallic Airplanes 


MILITARY 

REQUIREMENTS PER MIL-A-83444 

(AIRPLANE DAMAGE TOLERANd REQUIREMENTS. JULY 1974) 

COMMERCIAL 

MANUFACTURER/FAA 


Damage Tolerance Requirements 
Composite Structures 


• NO SPECIFIC REQUIREMENTS ESTABLISHED YET 


• WILL NOT ACCEPT COMPOS ITE STRUCTURES THAT 
ARE LESS DAMAGE TOLERANT THAN METALS 



Composite Structures 

DAMAGE CAUSES 


• MANUFACTURING 

• MAINTENANCE AND REPAIR 

• IN-SERVICE IMPACT 

• STRUCTURAL OVERLOAD 


Damage Tolerance 
Metallic vs Composites 

METALLIC 

FLAWS, CRACKS; TENSION CRITICAL 
METHODOLOGY FAIRLY WELL ESTABLISHED 


COMPOSITES 


MULTIPLE DELAMINATIONS, DISBONDS; 
COMPRESSION CRITICAL 
METHODOLOGY VIRTUALLY NON EXISTENT 



Composite Structures 


DO WE DESIGN STRUCTURE FOR UNLIMITED LIFE 
IN THE PRESENCE OF NONVIS IDLE DAMAGE? 


Composite Structures 
In-Service Impact 


DEGREE OF DAMAGE / INSPECTION 


• NONVISIBLE / ? 

• VISIBLE, BACK SURFACE ONLY / DEPOT 

• VISIBLE, BOTH SURFACES / FLIGHT EVIDENT 

GROUND EVIDENT 
SPECIAL VISUAL, ETC. 
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Composite Structure 
Out-of-Plane Deformations 

PRESSURIZATION , POST-BUCKLING 

SKIN/STRINGERS J 

SKIN/S UBSTRUCTURE ( ATTACHMENT? 

SPAR WEB/STIFFENERS ) 

CO-CURE; ADHESIVE BONDING; MECHANICAL FASTENERS; 
STITCHING; STAPLING; COMBINATIONS 



ANALYSIS EXPLAINS VARIABLE RATE OF INSTABILITY-RELATED 

DELAMINATION GROWTH 

John D. Whitcomb 
NASA Langley Research Center 


Under compression fatigue, delaminations in composites sometimes induce 
localized buckling, causing high interlaminar stresses at the ends of the delami- 
nation. Rapid delamination growth and loss of structural stability often ensue. 
Because delamination growth can lead to structural instability, the growth process 
must be understood. 

To Improve our understanding, through-width delamlnatlons (see figure) were 
studied experimentally and analytically. The figure shows a comparison of measured 
growth rates and calculated strain energy release rates, which are a measure of the 
intensity of stresses at the crack tip. In the figure, Gj and Gji are the 
energy release rates related to peel and shear stresses, respectively. Note that 
both Gx and the growth rate first increase then decrease rapidly with crack 
extension, whereas Gxi Increases monotonlcally and does not reflect the change in 
growth rate. Apparently the rate of delamination growth is governed by the intensity 
of the peel stress field. Hence, prediction of growth depends on an accurate 
assessment of the peel stress field around the crack tip. 


COMPARISON OF STRAIN-ENERGY RELEASE RATES 
AND DELAMINATION GROWTH RATES 
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THE APPROACH TO FAILURE ANALYSIS OF COMPOSITES 


Kenneth Reif snider 
Materials Response Group 

Virginia Polytechnic Institute and State University 


Much of the philosophy of failure analysis expressed at this workshop took the 
general direction of damage tolerance concepts. This is certainly a reasonable 
approach for structures, since, for example, a structure such as an airplane must be 
periodically inspected and the service of that structure continued (or discontinued) 
based on an interpretation of the residual strength and serviceability of the 
component. Such an approach requires nondestructive examination methods that can 
identify and quantitatively measure observables that are precisely and reliably 
related to engineering response. It is also essential that failure Induced by 
long-term loading be sufficiently well understood and represented so that accurate 
anticipations of residual strength, stiffness, and life (the materials response) can 
be made based on the observables. 

The experience of our Materials Response Group at VPI&SU and our perception of 
the results of other investigators suggest that the failure of composite laminates 
after long periods of loading at engineering levels can be typified by three rather 
distinct aspects. One of these aspects, which was highlighted at this workshop, and 
is the primary component of most current efforts to discuss damage tolerance. That 
aspect is stress redistribution. Currently it is typical, for example, to try to 
use nondestructive methods to find a major delamlnatlon, introduce the geometry of 
that defect into an approximate stress analysis of the component, and attempt to 
predict the residual strength of the component by typifying the failure resistance 
of the material using, say, a critical strain energy release rate or toughness 
determined from quasi-static laboratory tests. 

However, such an approach neglects the other two major aspects of long-term 
failure which we believe must be addressed. In addition to stress redistribution, 
we believe that material degradation and localization of damage are critical aspects 
of long-term damage development and consequent failure. Even under uniaxial uniform 
stress states, unnotched coupons (not dominated by edge effects) eventually fail 
under load levels which may be only 60 to 70 percent of the static ultimate (tensile) 
strength (i.e., the residual strength is reduced by 30 to 40 percent). As a result 
of recent work in our laboratory, we know that a very complex, widely distributed, 
and dense array of microdamage develops under such circumstances, and that there is 
a distinct localization process that precedes the final failure event. We also know 
that strength reductions of as much as 30 to 40 percent can result from these damage 
aspects, as we have mentioned. While we have made some preliminary efforts (and 
other investigators are undoubtedly making similar attempts), a philosophy and 
(appropriate) analysis have not yet been established such that these aspects can be 
described and anticipated by a rational and rigorous model. Our efforts continue in 
this direction. We believe that the ultimate success of failure analysis efforts, 
especially those based on a damage tolerance approach, will depend on the degree to 
which all three of the basic aspects of long-term damage development (stress 
redistribution, material degradation, and localization) are successfully investigated 
and characterized. 
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Finally, we also believe that there is a need for investigative and develop- 
mental efforts directed towards the mechanics of nondestructive testing and 
evaluation. We illustrate this opinion by mentioning two examples which emphasize 
an important contrast. The contrast is between methods which present analog or 
strictly geometric information (the most typical schemes) and methods which function 
based on the actual mechanical response of the material component. One example of 
the latter is a technique developed in our laboratory at Virginia Tech called 
vibrothermography. That technique, simply put, uses variable-frequency mechanical 
excitation to generate hysteresis energy heat patterns which are directly related to 
the stress distributions and the disturbances in mechanical response produced by 
Internal damage. The sensitivity and. Indeed, the physical basis of this technique 
are primarily mechanical, and therein lie the contrast and strength of this technique 
for the prediction of subsequent strength and life of a damaged laminate. The second 
example is the use of stiffness change as a method of detecting and interpreting 
damage. As micro-damage develops in composite laminates, stiffness components change 
in precise and significant amounts which reflect the nature and severity of damage. 
For long-term load histories, the engineering Young's modulus, for example, typically 
changes by the same order of magnitude as the residual strength (i.e., by as much as 
30 to 40 percent) even for fiber-dominated laminates. These stiffness changes enter 
the constitutive equations and, therefore, the analytical mechanics and modeling 
quite naturally. The measurement and interpretation of stiffness and stiffness 
changes are familiar exercises. 

We feel strongly that there is a great need for a research and development 
effort that addresses the need for better understanding and exploitation of the 
mechanics of nondestructive testing and evaluation. In summary, we believe that a 
proper failure analysis must consider the distribution of damage, the consequent 
local state of stress, and local state of the material. The analysis should be 
based on a sound rational mechanics development, which in turn is based on physical 
observables directly related to the mechanical response that is to be anticipated. 

If such a body of philosophy is to be successfully and completely developed, it will 
require that specialists in mechanics, materials, and structures work closely and 
effectively together. 
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DISCUSSION OF INTERLAMINAR STRENGTH EFFECTS ON 
AIRCRAFT DESIGN 




D. J. Watts 
Douglas Aircraft Co. 


A primary issue in the design of high— load-intensity primary aircraft structures 
is the identification and determination of small out-of-plane forces due to exter- 
nally applied or deflection— induced loads which produce interlaminar forces within 
the laminate. Interlaminar shear forces from in-plane loadings can also be 
significant in high-load-intensity structural designs in the vicinity of dis- 
continuities. These interlaminar forces can initiate and propagate delaminations 
which can drastically reduce the residual compressive strength of structural 
elements. Methodology and data are required to establish interlaminar static 
strength, delamination initiation and growth characteristics, the effect of flaws 
on these characteristics, and the residual compressive strength of delaminated 
structures. Resolution of these Issues not only affects the design but is important 
in guiding manufacturing technology development to avoid developing methods and 
processes which cannot be applied to designs that strive to minimize interlaminar 
forces. 
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